THE  BASIS  FOR  RAINFOREST  DIVERSITY  AND  BIOSPHERE  2 


By 
LINDA  SUSAN  LEIGH 


A  DISSERTATION  PRESENTED  TO  THE  GRADUATE  SCHOOL 

OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 

OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 

DOCTOR  OF  PHILOSOPHY 


UNIVERSITY  OF  FLORIDA 
1999 


ACKNOWLEDGMENTS 


Dr.  H.  T.  Odum  has  been  the  quintessential  mentor,  finding  no  idea  too 
trivial  to  explore  and  no  emergy  analysis  or  comparison  too  radical  to 
consider.     I  thank  him  for  his  encouragement  of  both  creative  imagination 
and  scholarship.   I  would  like  to  thank  my  committee  members  for  the 
inspiration  they  provided  throughout  my  doctoral  program:  Dr.  M.  T.  Brown, 
Dr.  C.  L.  Montague,  Dr.  C.  S.  Holling,  and  Dr.  F.  N.  Scatena. 

I  am  indebted  to  the  original  Biosphere  2  rainforest  design  team  with 
whom  I  worked  in  my  role  as  Biome  Design  Coordinator  from  1985  to  1993: 
Ghillean  Prance  (rainforest  biome  design  'captain'),  Tony  Burgess  (desert 
biome  design  'captain'),  Peter  Warshall  (food  web  design),  Scott  Miller 
(entomology),  Harry  Scott  (biome  design  and  species  selections),  John  Druitt 
(biome  design  and  species  selections),  Walter  Adey  (aquatic  systems),  Julia 
Bennett  (plant  propagation  and  management),  Bob  Scarborough  (soils), 
Stephen  Storm  (soils)  Phil  Hawes  (architectural  design),  Michael  Balick  (plant 
selections),  Andrew  Henderson  (plant  selection,  collections),  and  others. 
I  completed  the  third  plant  survey  with  the  assistance  of  Jessica  Bolson,  and 
measurements  in  1998  with  assistance  of  the  students  of  the  Spring  1998 
Columbia  University  Earth  Semester  at  Biosphere  2. 

Some  of  the  data  for  this  dissertation  were  assembled  from  the 
following  reports,  which  are  on  file  at  Columbia  University's  Biosphere  2 
Center  in  Oracle,  Arizona:  biomass  estimates  for  1990  and  1991  (Alan 
Haberstock,  Dec.  1991);  biomass  estimates  for  1993  (Mark  Bierner,  Oct.  1993); 


soil  chemical  and  physical  properties  in  1993  (Harry  J.  Scott,  Dec.  1994);  maps 
and  survey  data  (Jeremiah  Teague  and  Co.,  Sept.  1991);  soils  design,  assembly, 
and  placement  (Robert  Scarborough,  Mar.  1993).  Access  to  the  Biosphere  and 
logistic  support  was  facilitated  by  Adrian  Southern  and  John  Adams.  Eda 
Melendez  provided  data  from  EI  Verde. 

An  educational  grant  was  received  from  Edward  Bass,  a  scholarship 
fund  from  the  Drylands  Institute,  and  a  Graduate  Assistantship  in  the 
contract  between  the  Institute  of  Tropical  Forestry,  U.S.  Department  of 
Agriculture  and  the  University  of  Florida,  H.  T.  Odum,  Principal  Investigator. 
Unconditional  thanks  are  extended  to  Richard  Felger,  Silke  Schneider, 
Roy  Walford,  Kathleen  Dyhr,  Bernd  Zabel,  Dan  Levinson,  Randall  Gibson, 
and  Bonnie  Knickerbocker;  Carlos  Nagel,  the  late  Alfredo  Rivera,  Sherry  and 
Casey  Brandt-Williams,  Joanie  Breeze,  Tony  Burgess,  Jerry  Leigh;  and  a  long 
list  of  loyal  friends. 

"Never  say  die"  and  "Keep  the  voyage  going,  going,  but  never  gone" 
reflect  the  attitudes  of  me  and  my  fellow  crew  members  inside  Biosphere  2  - 
Gaie  Ailing,  Taber  MacCallum,  Jane  Poynter,  Roy  Walford,  Mark  Nelson, 
Sally  Silverstone,  and  Mark  Van  Thillo  -  and  they  continue  to  inspire. 

I  extend  a  special  acknowledgment  and  heartfelt  thanks  to  Ren  Hinks 
who  encouraged  me  on  a  daily  basis. 


TABLE  OF  CONTENTS 

page. 

ACKNOWLEDGMENTS ii 

LIST  OF  TABLES viii 

LIST  OF  FIGURES xii 

ABSTRACT xvii 

1    INTRODUCTION ! 

Review  of  Diversity  Influencing  Processes  in  Other  Studies ?. 

Species  and  Area 2 

Species  and  Individuals 6 

Succession 7 

Theories  of  Diversity 9 

Diversity  Decline 14 

A  Conceptual  Model  of  Ecosystem  and  Diversity 20 

Carrying  Capacity  of  Systems  for  Diversity 30 

Spatial  Organization 30 

Hierarchical  Organization 30 

Rainforest  in  Biosphere  2 31 

Design  Elements  and  Description 34 

A  Conceptual  Model  of  Diversity  in  Rainforest  of 

Biosphere  2 38 

Rainforest  in  Puerto  Rico 41 


2    METHODS.. 


..43 


Construction  and  Operation  of  the  Rainforest  in  Biosphere  2     43 

Procedures  Used  in  Starting  the  Ecosystem ZZ45 

Soil  Components  and  Placement 45 

Collection  and  Initial  Placement  of  Plants .....47 

Habitat  Assembly  ... '49 

Climate  Maintenance  Systems 53 

Human  Intervention c^ 

Plant  Mapping  and  Identifications ..,   57 

Initial  Mapping  and  Identifications:   1990-1991 ...     .     ."       S7 

Species  Additions  and  Removals 57 


Field  Measurements eg 

Second  Survey:  Transition  period,  1993-1994  58 

Third  Survey    June- August  1996.           _. 59 

Surveys  of  Species  Found  per  1000  Individuals  Counte        ...       61 

Leaf  Area  Index '^2 

Number  of  Seedlings  per  mJ 6? 

Fercent  of  Holes  in  Leaves     „.,  .             ..  63 

Number  of  Green  or  Yellow  Fallen  Leaves  per  m2 ........63 

Calculations ^ 

Diversity  Index ^0 

Biomass  Estimates ^ 

Growth  Form  Spectra 6c 

Poisson  Distribution 65 

Simulation  Method 66 

Emergy  Evaluation 67 

5    RESULTS 72 

Characterization  of  and  Changes  in  Biosphere  2  Rainforest  72 

Sci1 72 

Plants £ 

Spatial  Distribution  of  Plants .ur ^ 90 

Cumulative  Species  per  Cumulative  Individuals  Ccunted     90 

Leaf  Area  Index "g. 

Number  of  Seedlings  and  Green  and  Yellow 

Fallen  Leaves  per  m2 94 

Percent  of  Hcles  in  Leaves       . .......    ,..„ 97 

Litterfall  and  Decomposition  Rates  ,.    ,         .     ,  97 

Carbon  Uptake  and  Respiration '         '     99 

Cutting  and  Consumers '"" 

Diversity  Index 

Biomass  Estimates 

Comparisons  Between  Biosphere  2  and  Tabonuco  Forest  in 

Puerto  Rico . 

Physical  Environment 1no 

Soil JJ3 

plants :::::::::::::::::; ^ 

Overview  Simulation  Models 122 

Production  and  Biomass  Minimodel  Description  126 

Results  of  Simulating  the  Production  and  Biomass  Minimodei- 130 

1  roduction  and  Diversity  Minimodei  Description  130 

Results  of  Simulating  the  Production  and  Diversity  Minimodei 143 

Emergy  Evaluation  of  Rainforest  in  Biosphere  2  '  " 14g 


..101 
..101 


4   DISCUSSION.. 


..154 


Succession  with  Declining  Diversity        . 154 

Effect  of  First  Arrivals 155 

Limited  Access  to  Genes 156 

Declining  Diversity  and  a  Species  Plateau 156 

Extrapolation  of  Species  Composition  According  to 

Reproduction 158 

Simulation  of  Species  Decline 160 

Mechanisms  Affecting  Diversity  in  Biosphere  ?. = „.,  163 

Excess  Resources  and  Diversity 163 

Effect  of  Trimming  Weedy  Growth 164 

Pulsing  Disruption I65 

Consumers 166 

Graphical  Representation  of  Cumulative  Species-Individual 

Counts 166 

Effect  of  Spatial  Characteristics  of  Plants  in  Biosphere  2 169 

Comparisons  with  El  Verde I69 

Comparisons  Between  Rainforest  Structure  in  Biosphere  2  and 


El  Verde.., 


170 


Soil  Structure 170 

Absence  of  Most  animals  and  Simplification  of  Food  Webs .      171 

Plant  Reproduction 171 

Species  and  Individuals 172 

General  Implications 174 

Comparison  of  Biosphere  2  with  Role  and  Trends  of 

Global  Biodiversity 174 

Declining  Diversity  and  Downsizing 175 

Comparing  Succession  of  Diversity  in  Biosphere  2  with 

Global  Cultural  Change 175 

Carrying  Capacity  of  Systems  for  Diversity "..176 

APPENDICES 

A       BIOSPHERE  2  RAINFOREST  PLANT  MAPS 177 

B       BIOSPHERE  2  RAINFOREST  PLANT  LISTS 199 

C        PRODUCTION  AND  BIOMASS  MINIMODEL 306 

D        PRODUCTION  AND  DIVERSITY  MINIMODEL 323 

REFERENCES 332 

BIOGRAPHICAL  SKETCH M2 


LIST  OF  TABLES 
Table 


page 


1-1.      Definitions  and  concepts  used  in  this  dissertation 

(after  Odum  1996) 3 

1-2.      Symbols  for  the  energy  systems  language  used  in  this 

dissertation 22 

1-3.     Parts  and  processes  of  ecosystem  model  that  affect  species 

numbers  and  account  for  observed  declines 24 

2-1.      Chronology  of  Biosphere  2  rainforest  construction 

and  operation 44 

2-2.      Rainforest  soil  components  and  their  origins.  From 

Scarborough  (1994) 46 

2-3.      Percent  of  components  in  topsoil  mixture  specifications  and  their 
corresponding  habitats.  From  Scarborough  (1994).  Totals 
are  approximate 48 

2-4.  Spreadsheet  used  to  calculate  coefficients  for  single  tank 
biodiversity  model  in  Figure  2-3,  calibrated  at  steady 
state 69 

3-1.      Characteristics  of  Biosphere  2  rainforest  soils  from 

samples  taken  on  different  dates 73 

3-2.      Soil  bulk  density  and  percentage  coarse  fragments  in  Biosphere  2 
rainforest  from  samples  made  in  November  1993. 
From  Scott  (1999) 76 

3-3.      Total  number  of  species  and  individual  plants  seeded  in  the 

Biosphere  2  rainforest 73 

3-4.      Number  of  individual  plants  and  species  recorded  in  1991, 1993, 
and  1996  surveys  of  Biosphere  2  rainforest.  The  first  number' 
in  each  entry  includes  plants  from  the  1991  planting,  the  second 
number  from  the  1993  planting  and  third  is  plants  that  have 
self-propagated.  Species  reported  are  for  species  new  to 
the  rainforest gj 


3-5.      Leaf  and  branch  interceptions  per  observation  through  vertical 
points  in  Biosphere  2  rainforest.  Data  for  1993  from 
Odum  et  al.  (1993) 96 

3-6.      Distribution  of  seedlings  found  in  30,  0.54  m2  circular  plots 

in  Biosphere  2  rainforest  understory  in  1998 96 

3-7.      Distribution  of  green  and  yellow  leaves  in  0.54  m2  circular  plots 

on  lowland  rainforest  ground,  April  1998 98 

3-8.      Percent  of  holes  in  leaves  in  Biosphere  2  rainforest  trees. 

Data  for  1993  from  Odum  et  al.  (1993) 98 

3-9.      System  metabolism  of  Biosphere  2  rainforest  for  summer,  1996. 

NEE=Net  ecosystem  exchange,  R^Soil  respiration,  Ac=Canopy 

assimilation,  RUE=Radiation  use  efficiency,  Rp=Plant 

respiration.  From  Lin  et  al.  (1999) 100 

3-10.    Shannon-Wiener  index  of  diversity  in  Biosphere  2  rainforest  for 

1991, 1993,  and  1996 102 

3-11.    Aboveground  plant  biomass  estimates  in  Biosphere  2  rainforest 
for  1990, 1991,  and  1993.  From  Haberstock  (1991)  and 
Bierner  (1993) 102 

3-12.  Environmental  variables  at  El  Verde  and  Biosphere  2  forests. 
From  Odum  (1970),  Romer  (1985),  Ahrain  et  al.  (1998),  and 
Cuevas  et  al.  (1991) 104 

3-13.    Comparison  of  soil  properties  of  Puerto  Rico  and  Biosphere  2 
forests.  From  Edmisten  (1970),  Silver  et  al.  (1994), 
Scott  (1999),  and  Lin  et  al.  (1998  and  1999) 106 

3-14.  Soil  bulk  density  in  Biosphere  2  rainforest  and  Puerto  Rico 
tabonuco  forest.  From  Scott  (1999),  Edmisten  (1970),  and 
Silver  et  al.  (1994) l0g 

3-15.    Plant  species  occurring  in  El  Verde  that  were  originally 
planted  in  Biosphere  2  forest,  with  an  inventory  of 
individuals  counted  per  species  in  1991  and  1996 109 

3-16.    Sampling  sites  for  species  /individual  counts,  and  number 

of  species/ 1000  individuals 113 

3-17.    Leaf-area  indices  of  Puerto  Rico  tabonuco  forest  and 

Biosphere  2  rainforest.  From  Odum(1970),  Jordan  (1969), 

Odum  and  Pigeon  (1970),  and  Odum  et  al.  (1993) 120 


3-18.    Seedling  density  in  Biosphere  2  and  tabonuco  rainforest. 

Tabonuco  forest  data  from  Odum  (1970b) 121 

3-19.     Percent  of  holes  in  leaves  in  Biosphere  2  and  tabonuco 

forest  trees.  Tabonuco  forest  data  from  Odum  (1970b) 121 

3-20.    Litterfall  rates  in  Biosphere  2  rainforest  habitats  and  tabonuco 
forests.  From  Nelson  (1999),  Wiegert  (1970), 
Cuevas  et  al.  (1991),  Lugo  (1992),  and  Lodge  (1991) 123 

3-21.  Decomposition  rates  of  litter  in  Biosphere  2  rainforest  and 
tabonuco  forest  in  Puerto  Rico.  From  Nelson  (1999)  and 
Wiegert  (1970) 123 

3-22.    Nutrient  content  of  leaves  in  percent  of  dry  weight  in 

Biosphere  2  and  El  Verde  forests.  From  Odum  (1970a), 

Ovington  (1970),  Medina  et  al.  (1981),  and  Lin  et  al.  (1999) 124 

3-23.    Carbon  exchange  in  Biosphere  2  (summer)  and  tabonuco 
rainforests.  R^Soil  respiration,  Ac=Canopy  assimilation 
From  Lin  et  al.  (1999) 125 

3-24.    Estimates  of  aboveground  biomass  density  in  Biosphere  2  and 
tabonuco  forests.  From  Haberstock  (1991),  Bierner  (1993), 
Ovington  and  Odum  (1970),  Scatena  et  al.  (1993),  and 
Lugo  (1992) 125 

3-25.    Accumulated  emergy  inputs  to  Biosphere  2  rainforest  for 

start-up  of  the  system  prior  to  material  closure  in  1991 150 

3-26.    Comparison  of  emergy  in  developing  rainforest  in  Biosphere  2 

in  1991  and  25-year-old  natural  tabonuco  forest 153 

4-1.      Extrapolation  of  species  composition  of  Biosphere  2  rainforest 
to  100  years  with  30  species.  H=Herb,  T=Tree,  P=Palm, 
G=Giantherb,  A=Bamboo,  C=Climber,  S=Shrub 161 

4-2.      Species  / 1000  individuals  and  slope  of  log-log  graphs  for  all 

study  sites.  K=slope  of  log-log  graph 173 

B-l.      Individual  plants  in  the  Biosphere  2  rainforest 199 

B-2.      Species  from  the  first  planting  of  the  Biosphere  2  rainforest, 
with  inventories  from  1991, 1993  and  1996.  T=Tree, 
S=Shrub,  P=Arboreal  palm,  R=Graminoid,  C=Climber, 
A=Woody  graminoid,  such  as  bamboo,  H=Herb, 
G=Giantherb,  E=Epiphyte 283 


B-3.      Species  from  the  second  planting  of  the  Biosphere  2 
rainforest,  with  inventories  from  1993  and  1996,  or 
self-propagated.  T=Tree,  S-Shrub,  P=Arboreal  palm, 
R=Graminoid,  C=Climber,  A=Woody  graminoid,  such  as 
bamboo,  H=Herb,  G=Giantherb,  E=Epiphyte 299 

C-l.     EXCEL  spreadsheet  used  to  calculate  coefficients  for 

Biosphere  2  rainforest  production  and  biomass  minimodel 

under  predicted  steady  state  conditions 308 

C-2.     Program  in  BASIC  for  the  simulation  of  metabolism  of  the 

rainforest  in  Biosphere  2,  Model  B2METAB  in  Figure  3-19  313 

D-l.     EXCEL  spreadsheet  used  to  calculate  coefficients  for 

Biosphere  2  rainforest  diversity  and  production  minimodel 

SPDIV 324 

D-2.  Program  in  BASIC  for  the  simulation  of  production  and 
diversity  in  the  Biosphere  2  rainforest,  Model  SPDIV 
in  Figure  3-28 327 


LIST  OF  FIGURES 

Figure 
**^  page 

1-1.      Systems  diagram  representing  theories  of  ecosystem  diversity 21 

1-2.      Areas  within  Biosphere  2 32 

1-3.  Continuous  photosynthetically  active  radiation,  C02, 
and  temperature  in  terrestrial  wilderness  biomes  of 
Biosphere  2  during  the  first  material  closure 33 

1-4.      Side  view  of  rainforest  biome        .     . 35 

1-5.      Habitats  delineated  within  the  Biosphere  2  rainforest  37 

1-6.      Systems  diagram  representing  ecosystem  diversity  within 

Biosphere  2  rainforest 39 

2-1.      Hydrologic  connection  of  Biosphere  2  rainforest  to  the  rest  of 

the  Biosphere  2  system ^ 

2-2.  Diagram  and  equations  for  species  diversity  on  an  island  as 
described  by  MacArthur  and  Wilson  (1967)  and  drawn  bv 
Beyers  and  Odum  (1993) 6g 

2-3.      Single  tank  model  of  biodiversity  as  a  balance  of  steady  inflow 

pathway  coefficients 6g 

2-4.  Single  tank  model  of  biodiversity  with  calibration  values  for 
flows  and  storages 

2-5.      Simulation  of  number  of  species  Q  for  different  starting  values 

using  the  single  storage  model  in  Figure  2-3 „ 70 

3-1.      Map  section  #3  showing  location,  growth  form,  and  canopv 
size  of  individual  plants  in  the  Biosphere  2  rainforest  in 
September  1991.   The  twenty  map  sections  are  given  in 
Figure  A-l.  Approximate  scale:   1  cm=0.6  m 79 

3-2.  Change  in  plant  species  abundance  and  distribution  in 
Biosphere  2  rainforest,  1991-1996,  for  plants  from  the 
first  planting 


3-3.      Number  of  individual  plants  in  Biosphere  2  rainforest 

time  including  2  planting  periods.  1991  and  1993  .  84 

3-4.      Number  of  species  in  Biosphere  2  rainforest  through  time  and 

with  two  planting  periods,  1991  and  1993 84 

3-5.      Relationship  of  number  of  species  to  number  of  individuals 
over  time  in  Biosphere  2  rainforest  on  linear  (upper  figure) 
and  semilog  (lower  figure)  scale.  Open  triangles  are  plants 
from  1991  planting,  squares  from  1993  planting,  and  circles 
combine  plants  from  both.   Survey  data  are  from 
1991, 1993  and  1996 86 

3-6.      Number  of  individuals  within  growth  forms  in  Biosphere  2 
rainforest  in  1991,  1993  and  1996.  Data  include  only  plants 
from  first  planting.   H=Herb,  T=Tree,  S=Shrub,  C=CIimber, 
G=Giant  herb,  P=Arborescent  palm,  E=Epiphyte, 
R=Graminoid,  A=Bamboo 87 

3-7.      Number  of  species  within  growth  forms  in  Biosphere  2  rainforest 
in  1991,  1993  and  1996.  Data  include  only  plants  from  first 
planting.  H=Herb,  T=Tree,  S=Shrub,  C=Climber,  G=Giant 
herb,  P=Arborescent  palm,  E=Epiphyte,  R=Graminoid, 
A=Bamboo RR 

3-8.      Change  in  number  of  species  within  growth  forms  from  1991 

planting  in  Biosphere  2  rainforest  through  time     .  .89 

3-9.      Change  in  distribution  of  species  within  growth  forms  from  1991 

planting  in  Biosphere  2  rainforest  through  time 90 

3-10.    Spatial  distribution  of  plants  in  Biosphere  2  rainforest  in 

1991  (solid  squares)  compared  to  Poisson  distribution  (open 
triangles) *  F  gj 

3-11.    Cumulative  plant  species  as  a  function  of  cumulative 

individuals  in  Biosphere  2  rainforest,  1991.  Data  for  all 
habitats  except  cliff  faces  and  surface  aquatic  systems  are 
plotted  on  linear  (a),  semi-log  (b)  and  log-log  (c)  scales 92 

3-12.    Cumulative  plant  species  as  a  function  of  individuals  in 
Biosphere  2  rainforest,  April  1998.   Open  squares  have 


Scindapsus  and  Syngonium  subtracted.. 


..93 


3-13.    Cumulative  species  as  a  function  of  cumulative  individuals 
in  Biosphere  2  rainforest  in  1991  (solid  circles),  1998a 
(all  plants,  open  squares)  and  1998b  {Syngonium  and 
Scindapsus  deleted,  crosses) 95 

3-14.  Growth  form  spectra  of  Biosphere  2  rainforest  and  El  Verde 
tabonuco  forest.  B2=Biosphere  2.  Data  for  EI  Verde  from 
Smith  (1970) m 

3-15.    Cumulative  plant  species  as  a  function  of  cumulative 
individuals  >0.5  m  tall  counted  on  an  11-year-old 
landslide  at  El  Verde  site H4 

3-16.    Cumulative  plant  species  as  a  function  of  cumulative 

individuals  counted  in  Biosphere  2  rainforest  (solid  circles) 
and  on  an  11-year-old  landslide  at  El  Verde  site 
(open  squares) 215 

3-17.    Cumulative  plant  species  as  a  function  of  cumulative 

individuals  counted  on  El  Verde  landslide  (open  circles), 

Radiation  Center  (open  triangles),  and  Bisley  grid  (crosses) 116 

3-18.    Cumulative  plant  species  as  a  function  of  cumulative 

individuals  counted  on  EI  Verde  landslide,  Radiation 

Center,  Bisley  grid  and  Biosphere  2  rainforest  for 

1991  and  1998 U8 

3-19.    Systems  diagram  for  simulating  production  and  biomass  in 
Biosphere  2  rainforest  showing  storages,  pathway 
coefficients,  and  equations 127 

3-20.    Systems  diagram  for  simulating  production  and  biomass  in 

Biosphere  2  rainforest  showing  calibration  values  projected 

for  steady  state  after  20  years  of  operation 131 

3-21.    Ten-year  simulation  of  production  and  biomass  minimodel 
B2METAB  in  Figure  3-19.   Final  values  for  plant  biomass ' 
are  on  the  right  side  of  their  graphs,  g/m2  -  grams  per 
square  meter,  dry  weight  for  biomass 132 

3-22.    One-hundred-year  baseline  simulation  of  production  and 
biomass  in  the  minimodel,  B2METAB  in  Figure  3-19. 
Final  values  for  plant  biomass  are  on  the  right  side  of 
their  graphs,  g/m2  =  grams  per  square  meter,  dry  weight 
for  biomass ° 


3-23.    Comparison  of  simulations  of  biomass  by  the  model  in 

Figure  3-19  showing  effects  of  light,  elimination  of  weedy 

biomass,  limited  pruning,  and  altered  airflow  on  mature 

and  weedy  plant  biomass.  Values  for  biomass  at  100  years 

are  given  on  the  right  side  of  the  graphs  in  grams  of  dry 

biomass  per  square  meter.  Scales  are  the  same  for  each 

box.   M=Marure  plants,  W=Weedy  plants 135 

3-24.    Comparison  of  effects  of  light,  elimination  of  weedy  biomass, 
limited  pruning,  and  altered  airflow  on  system  production 
ratio,  Pg/Rp.  Scales  are  the  same  for  each  box,  with  the 
center  line  =  1 136 

3-25.    Systems  diagram  for  simulating  diversity  and  production  in 
Biosphere  2  rainforest  showing  storages,  pathway 
coefficients  and  equations 140 

3-26.    Systems  diagram  for  simulating  production  and  diversity  in 
Biosphere  2  rainforest  showing  calibration  values  for 
storages  and  flows 144 

3-27.    Baseline  simulation  of  diversity  and  production  minimodel 
SPDIV,  in  Figure  3-25  for  Biosphere  2  rainforest  for 
10  years.  Final  values  for  plant  biomass  and  species  are  on 
the  right  side  of  their  graphs,  g/m2  -  grams  of  dry  biomass 
per  square  meter 145 

3-28.    Baseline  simulation  of  diversity  and  production  minimodel 
SPDIV,  in  Figure  3-25  for  Biosphere  2  rainforest  for  100 
years.  Final  values  for  plant  biomass  and  number  of  species 
are  on  the  right  side  of  their  graphs,  g/m2  =  grams  of  dry 
biomass  per  square  meter 147 

3-29.     Emergy  analysis  summary  diagram  of  Biosphere  2  rainforest 149 

4-1.      Extrapolation  of  number  of  plants  in  El  Verde  rainforest  to  a 
possible  diversity  plateau  (solid  triangle)  for  the 
Biosphere  2  rainforest.   Measured  data  from  Biosphere  2 
are  shown  with  solid  circles  (1991)  and  open  squares  (1998)  159 


A-l. 


Map  cells  showing  the  location  of  individual  plants  in  the 
Biosphere  2  rainforest 


..177 


C-l.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with 
light  increased  20%  over  baseline.   Final  values  for  plant 
biomass  are  on  the  right  side  of  their  graphs.  g/m2=grams 
per  square  meter,  dry  weight  for  biomass 316 

C-2.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with 
no  weedy  biomass.   Final  values  for  plant  biomass  are 
on  the  right  side  of  their  graphs.  g/m2=grams  per  square 
meter,  dry  weight  for  biomass 317 

C-3.      One-hundred-year  simulation  of  Biosphere  2  rainforest 

production  and  biomass  minimodel,  B2METAB,  where 

all  pruned  biomass  is  put  onto  the  soil.   Final  values  for 

plant  biomass  are  on  the  right  side  of  their  graphs. 

g/m2=grams  per  square  meter,  dry  weight  for  biomass 318 

C-4.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with 
no  pruning  by  humans.  Final  values  for  plant  biomass 
are  on  the  right  side  of  their  graphs.  g/m2=grams  per 
square  meter,  dry  weight  for  biomass 319 

C-5.       One-hundred-year  simulation  of  Biosphere  2  rainforest 

production  and  biomass  minimodel,  B2METAB,  human 

effort  is  reduced  to  .75  of  the  baseline.  Final  values  for 

plant  biomass  are  on  the  right  side  of  their  graphs. 

g/m2=grams  per  square  meter,  dry  weight  for  biomass 320 

C-6.     One-hundred-year  simulation  of  Biosphere  2  rainforest 

production  and  biomass  minimodel,  B2METAB,  human 

effort  is  reduced  to  .85  of  the  baseline.  Final  values  for 

plant  biomass  are  on  the  right  side  of  their  graphs. 

g/m2=grams  per  square  meter,  dry  weight  for  biomass 321 

C-7.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with 
airflow  cut  off  from  the  rest  of  the  Biosphere.   g/m2=grams 
per  square  meter,  dry  weight  for  biomass 322 

D-l.     Simulation  of  Biosphere  2  rainforest  with  continuous 

additions  of  5  and  20  species  per  year  over  10  years.  Final 
values  for  plant  biomass  and  species  are  on  the  right  side 
of  their  graphs.  g/m2=grams  of  dry  biomass  per  square  meter         329 


D-2.    Simulation  of  Biosphere  2  rainforest  with  continuous 

additions  of  5  and  20  species  per  year  over  100  years.  Final 
values  for  plant  biomass  and  species  are  on  the  right  side 
of  their  graphs,  g/  m2=grams  of  dry  biomass  per  square  meter 330 

D-3.     Simulation  of  Biosphere  2  rainforest  when  50  species  are 

added  after  the  second  year.  Final  values  for  plant  biomass 

and  species  are  on  the  right  side  of  their  graphs,  g/  m2=grams 

of  dry  biomass  per  square  meter 331 


Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 


THE  BASIS  FOR  RAINFOREST  DIVERSITY  AND  BIOSPHERE  2 

By 

Linda  Susan  Leigh 

August  1999 

Chairman:  H.  T.  Odum 

Major  Department:  Environmental  Engineering  Sciences 

A  miniature  rainforest  created  in  the  glass-enclosed  Biosphere  2 
mesocosm  in  Arizona  was  compared  with  montane  rainforest  in  Puerto  Rico  to 
study  the  basis  for  biodiversity  and  succession.  Initial  seeding  of  high  plant 
diversity  was  a  technique  for  determining  by  extrapolation  the  carrying  capacity 
in  Biosphere  2  and  in  disturbed  areas  in  Puerto  Rico.  In  Biosphere  2  the 
rainforest  area  was  started  in  1991  with  1890  plants  of  316  species  in  0.19 
hectares.  Three  hundred  and  thirty-nine  plants  in  92  species  were  added  in  1993. 
A  diversity  index,  species  per  1000  individuals  counted,  decreased  from  250  in 
1991  to  96  in  1998  compared  to  60  in  the  comparison  forest  while  the  Shannon- 
Wiener  index  decreased  from  an  initial  5.39  to  4.64  in  1996  compared  to  4.62  in 
the  comparison  forest.  Although  normal  populations  of  insects  and  pollinators 
were  absent,  vegetative  reproduction  established  productivity,  diversity, 


hierarchy  and  developed  soil  profiles  approaching  those  in  the  comparison 
forest. 

Graphs  of  cumulative  species  versus  individuals  counted  were  analyzed 
to  identify  mechanisms  affecting  diversity  including  seeding,  competition  for 
excess  nutrients  and  carbon  dioxide,  extinction  rate,  seasonal  change  in 
insolation,  and  pruning  and  management  by  humans. 

The  productive  basis  for  diversity  was  simulated  with  a  model  of 
production,  consumption,  and  recycle  that  included  management  alternatives. 
Simulation  runs  showed  total  biomass  production  highest  when  the  weedy  vine 
biomass  was  regularly  trimmed  and  removed. 

Another  model  of  main  factors  affecting  diversity  was  simulated  with  a 
minimum  species  limit  to  help  explain  the  observed  patterns  and  anticipate 
species  carrying  capacity.  With  an  addition  of  50  new  species  after  2  years,  the 
simulation  showed  an  increase  in  diversity  and  decrease  in  biomass  after  10 
years  when  compared  to  the  simulation  without  further  additions.  Emergy  and 
emdollar  evaluation  showed  resources  required  for  developing  the  rainforest  in 
Biosphere  2  were  2300  times  larger  than  those  for  the  natural  succession. 


CHAPTER  1 
INTRODUCTION 


The  basis  for  biodiversity  is  a  principal  question  in  ecology  and  other 
fields.   Whereas  most  studies  of  succession  and  development  of  ecosystems 
have  concerned  periods  of  increase  and  diversification,  less  studied  are  the 
limits  to  diversity.    Are  there  general  systems  principles  controlling  the  self- 
organization  of  diversity  in  all  systems?   What  happens  when  there  is  initial 
seeding  of  high  diversity  and  limited  area?  In  this  dissertation  the  self- 
organization  of  the  rainforest  in  Biosphere  2,  the  glass-enclosed  mesocosm  in 
Arizona,  is  compared  with  the  tabonuco  rainforest  in  Puerto  Rico.   Floristic 
studies,  measurements,  and  models  are  used  to  account  for  patterns  of  growth 
and  diversity. 

Suppose  a  system  is  defined  to  include  a  newly  added  set  of  individual 
organisms  and  species.   Immediately,  relationships  develop  between  these 
and  the  flows  of  energy  and  matter  and  some  between  each  other.  Parts  and 
relationships  can  be  represented  with  systems  diagrams  that  give  an 
overview  of  the  whole  and  the  parts.   Diversity  and  structure  are  the 
information  of  the  system  that,  like  other  storages,  requires  energy  for  its 
maintenance.   The  amount  and  quality  of  energy  available  to  a  system  in  part 
determines  which  of  the  energy  requirements  for  diversity,  structure,  and 
their  connections  can  be  supported,  thus  shaping  its  biodiversity.   Science  of 
the  parts  of  ecosystems,  such  as  physiology,  and  science  of  the  whole,  such  as 
systems  ecology,  are  converging  with  respect  to  understanding  biodiversity. 


To  improve  the  synthesis  when  science  of  the  parts  is  applied,  ecosystems  and 
energy  are  included  in  systems  models. 

This  dissertation  considers  the  self-organization  of  the  diversity  of 
systems  with  limited  resources.  Do  systems  have  a  carrying  capacity  for 
diversity?   How  do  resources  and  mechanisms  of  system  interaction  account 
for  observed  patterns?  The  next  section  is  a  literature  review  of  the  theories 
and  processes  that  have  been  offered  to  account  for  decreases  in  diversity  and 
structure.  In  order  to  clearly  relate  diversity  to  its  resource  basis,  Table  1-1 
provides  definitions  of  the  measures  used. 

Review  of  Diversity-Influencing  Processes  in  Other  Studies 

Observations  of  diversity  patterns  over  scales  of  time,  space,  or  energy 
inputs  suggest  causal  mechanisms  affecting  the  processes  that  contribute  to 
diversity.    A  comparative  approach  is  often  used  in  studying  diversity  in 
ecological  systems. 

Species  and  Area 

Diversity  is  the  number  of  kinds  of  units  in  an  area.   In  an  ecosystem  it 
may  be  the  number  of  species  found.  Patterns  of  diversity  have  been 
observed  and  compared  for  decades  over  spatial  scales  ranging  from  the  entire 
biosphere  (latitudinally)  to  continents,  regions  and  localities 
(Rosensweig  1995). 

A  general  pattern  in  the  relationship  of  species  and  area  show  that 
within  a  taxonomic  group,  large  areas  support  more  species  than  small  areas. 
A  clear  description  of  the  pattern  was  first  published  by  Jaccard  (1912),  and  in 
1921  Arrhenius  published  a  mathematical  description  of  the  curve,  which 


Table  1-1.   Definitions  and  concepts  used  in  this  dissertation  (after  Odum 
1996). 


TERM 


DEFINITION 


Emergy 


Power 

Empower 

Transformity 

Species  richness 
Diversity 


Diversity  indices: 
Shannon  index 


Margalef  index 


Species/ 1000 
individuals 


The  amount  of  available  energy  of  one  particular  type 
that  was  used  up,  directly  and  indirectly,  to  make 
something  (product  or  service).   The  unit  with  the 
prefix  "em"  is  used  with  the  units  summed  to  indicate 
that  it  is  an  emergy  rather  than  an  energy  scale.   This 
study  uses  solar  joules  as  energy  units  and  solar 
emjoules  as  emergy  units. 

Flow  of  useful  energy  per  unit  time. 

Flow  of  emergy  per  unit  time. 

A  measure  of  energy  quality  defined  as  the  emergy  per 
unit  available  energy  of  a  product  or  service.  Units  are 
emjoules  per  joule. 

The  number  of  species  found  in  an  area.  This  is  often 
used  interchangeably  with  diversity,  though  richness 
only  accounts  for  the  number  of  species  and  not  their 
distribution  in  an  area. 

The  accounting  of  different  types  within  a  system 
using  taxonomic,  functional,  size,  trophic  levels, 
information  content,  etc.,  as  the  categories.  The' 
number  of  species  and  their  relative  abundances 
within  an  area  are  often  used  to  calculate  an  index  of 
diversity  using  various  formulas. 

H'  =  -I  p  In  Pi 

Pi   is  the  proportion  of  individuals  found  in  the  z'th 

species. 

(S-l)/ln  N 

S  =  number  of  species,  N  =  number  of  individuals 

A  form  of  the  Margalef  index 


Table  1-1  --  continued,. 


TERM 


Logarithmic 
index 


DEFINITION 


logS/logN 


Information 


Maximum 

empower 

theory 


System 


Energy  hierarchy 


The  units,  connections,  and  configurations  of  a  system 
requiring  some  form  of  energy  as  their  carrier. 
Information  is  useful  if  it  can  make  its  system  operate. 
Genetic  codes  are  an  example  of  useful  information. 

Self-organizing  systems  prevail  by  developing  designs 
that  maximize  emergy  inputs  and  effective  use  (Odum 
1996). 

A  part  of  nature  with  components  defined  by  artificial 
or  real  boundaries,  where  'windows'  of  time  and  space 
delimit  the  system  of  study.   Single  organisms  or 
single  populations  are  systems  since  they  have 
relationships  to  their  external  factors  and  among  their 
parts.  A  set  of  organisms  is  a  system  even  if  the  units 
are  new  and  not  interacting  much  with  each  other. 
They  are  interacting  with  the  light,  material  and  heat 
flows.  A  system  can  be  any  size. 

The  natural  arrangement  of  energy  transformation 
processes  in  a  chain  in  which  each  process  uses  energy 
of  one  kind  to  generate  less  energy  of  a  different  kind. 
Available  energy  decreases  along  the  chain,  but 
transformity  increases. 


was  quickly  refuted  by  Gleason  (1922).    Extensive  dialogue  regarding  the 
significance  of  the  mathematical  relationship  has  followed. 

A  plot  of  the  increase  of  species  as  the  sample  size  increases  within  an 
ecosystem  is  called  the  species /area  curve.  The  first  species /area  curve  is 
attributed  to  H.  C.  Watson  in  1835  (Connor  and  McCoy  1979).  The  curve  is 
used  to  indicate  the  number  of  species  within  an  area  of  study  larger  than  the 
sampled  area;  to  determine  the  rate  at  which  species  accumulate  within  an 
area  for  comparison  with  other  areas;  and  to  determine  the  minimum 
sampling  size  needed  to  adequately  characterize  diversity  of  an  area.  Vestal 
(1949)  attempted  to  calculate  equivalent  reference  areas  for  vegetation  types  in 
about  240  areas  worldwide.  The  reference  areas  were  used  to  compare  plant 
diversity  within  the  area  as  well  as  to  calculate  sampling  size  required  to 
characterize  the  vegetation  type.  The  following  site  differences  were  cited  as 
reasons  for  the  necessity  of  varying  the  area  sampled  to  determine  diversity: 
size  of  plants  or  clumps  of  plants;  degree  of  species  dominance;  low  density  of 
plants  overall;  productivity  of  the  site;  size  and  shape  of  the  site;  degree  of 
heterogeneity  of  the  environment;  complexity  of  the  community;  floristic 
richness;  abundance  among  species;  successional  stage;  stand  condition. 
Studies  of  diversity  have  to  be  carefully  related  to  area.   Though  these  reasons 
are  plausible  explanations  for  the  need  of  varying  sample  sizes,  his  results 
have  been  considered  "dubious",  "purely  hypothetical  and  probably  wrong" 
due  to  methods  he  used  in  data  conversions  (Goodall  1952). 

Though  the  pattern  of  diversity  increasing  with  area  may  seem 
obvious  or  trivial,  Huston  (1994,  p.  35)  claims  that,   "...the  underlying 
mechanisms  include  most  of  those  that  are  potentially  important  in 
regulating  diversity."    Three  hypotheses  were  suggested  by  Connor  and 
McCoy  (1979)  in  examining  the  causal  mechanisms  behind  the  pattern.   The 


first  is  considered  the  null  hypothesis,  and  states  that  the  pattern  is  a  statistical 
sampling  phenomenon  (Gotelli  and  Graves  1996).    The   alternate  hypotheses 
are  the  result  of  biological  processes.   The  first  suggests  that  diversification  of 
species  is  greater  in  heterogeneous  habitats,  which  occur  with  greater 
frequency  as  total  area  increases.  The  second  suggests  that  the  extinction  rate 
is  higher  in  smaller  areas  where  populations  are  smaller,  and  that 
immigration  rates  don't  change  between  areas. 

As  the  area  over  which  species  occur  increases,  the  available  resources 
for  species  support  increase  proportionally.   Species-area  curves  comparing 
local  and  regional  scales  show  the  species  increasing  as  the  total  energy 
available  to  support  diversity  increases.  In  a  study  of  82  different  variably 
sized  terrestrial  ecosystems,  Orrell  (1997)  found  correlation  between  species 
richness  and  the  energy  flow  of  each  ecosystem. 

Due  to  difficulties  in  comparing  species /area  curves  where  different 
areas  have  been  sampled,  other  researchers  have  recommended  using  a  plot 
of  the  number  of  species  that  accumulate  for  the  number  of  individuals 
counted.  Research  by  Condit  et  al.  (1996b)  in  tropical  forests  demonstrate  that 
diversity  estimates  can  be  compared  over  sites  where  an  identical  number  of 
stems  have  been  counted. 

Species  and  Individuals 

A  common  pattern  measured  in  nature  is  the  relationship  between  the 
cumulative  number  of  species  found  and  the  individuals  counted  within  a 
system.  The  typical  pattern  is  a  curve  of  decreasing  slope  (example  Figure  3- 
11),  and  is  the  same  type  of  pattern  as  the  species /area  curve  (Huston  1994). 
The  pattern  holds  for  total  enumeration  of  individuals  within  a  site  as  well  as 


for  subsets,  such  as  trees  within  a  particular  stem  diameter  range  (Condit  et  al. 
1996b).   Documentation  of  this  pattern  has  been  repeated  with  various  taxa 
such  as  aufwuchs  in  natural  springs  and  a  laboratory  microcosm  (Odum  and 
Hoskin  1957;  Yount  1956),  fishes  (Angermeier  and  Schlosser  1989),  and  insects 
(Fisher  et  al.  1943). 

Succession 

Succession  is  the  process  of  self-organization  of  an  ecosystem  after  a 
disturbance.    It  is  linked  to  diversity  in  that  species  composition  and 
dominance  change  in  an  area  over  time  due  to  processes  within  the  system  or 
processes  at  a  larger  scale.    Clements  (1936)  initially  suggested  succession  as  a 
phenomenon  whereby  early  successional  species  facilitate  later  successional 
species.  According  to  his  conceptual  model,  each  stage  of  succession  paves  the 
way  for  the  next  stage  by  altering  soils,  humidity,  nutrients  and  water 
availability.  Connell  and  Slatyer  (1977)  suggest  3  categories  of  causal 
mechanisms  for  successional  changes,  adding  to  Clements'  facilitation 
category:     facilitation,  inhibition,  and  tolerance.     With  their  model,  the  net 
effect  of  an  earlier  successional  species  on  the  establishment  of  a  later  species 
can  be  either  positive  (facilitation),  negative  (inhibition),  or  neutral 
(tolerance).     A  neutral  effect  allows  equal  probability  of  all  organisms  to 
become  established  in  an  area,  and  "leads  to  a  community  composed  of  those 
species  most  efficient  in  exploiting  resources,  presumably  each  specialized  on 
different  kinds  or  proportions  of  resources."  (Connell  and  Slatyer  1977).   All 
three  of  these  mechanisms  along  with  life  history  characteristics  of  organisms 
are  important  in  any  given  ecosystem  (Ricklefs  1990). 


Gutierrez  and  Fey  (975)  developed  a  model  simulating  the  general 
dynamic  patterns  of  succession,  relating  ecosystem  structure  to  successional 
patterns.  Simulations  of  their  model  suggest  that  the  internal  feedback 
structure  of  an  ecosystem  drive  secondary  succession.  Ulanowicz  (1980)  uses 
community  flow  networks  through  ecosystems  to  describe  the  development 
and  information  content  of  the  systems,  suggesting  that  flows  by  themselves 
adequately  describe  ecosystem  succession. 

Four  stages  of  succession  are  suggested  by  Holling  (1986):  exploitation, 
conservation,  creative  destruction,  and  renewal  (recycle).     During  the 
exploitation  and  conservation  phases,  an  ecosystem  develops  structure, 
functions,  and  diversity  from  available  energy  and  matter.   Disturbances  -  the 
creative  destruction  phase  -  occur  as  pulses  on  widely  varying  spatial  and 
time  scales  which  are  generally  characteristic  of  a  given  ecosystem. 
Disturbances  may  be  caused  by  abiotic  elements  or  by  biotic  factors.  The  pulse 
of  nutrients  made  available  allows  their  recycle  into  a  newly  developing 
ecosystem,  which  may  have  very  different  configurations  than  the  previous 
ecosystem  (Holling  1986). 

When  a  disturbance  is  widespread,  the  successional  process  may  be  in 
the  same  phase  over  the  entire  area  of  disturbance.   In  the  case  of  treefalls,  the 
successional  process  starts  over  in  small  local  areas,  creating  a  discontinuous, 
mosaic  landscape.  Remmert  (1991)  calls  this  a  'mosaic-cycle',  where  various- 
aged  patches  of  different  sizes  occur  throughout  an  ecosystem.  The 
importance  of  spatial  and  temporal  landscape  patterns  to  vertebrate  species 
diversity  has  been  documented  for  avian  species  in  Panama  by  Karr  and 
Freemark  (1983). 

Diversity  of  an  ecosystem  is  generally  reduced  after  the  creative 
destruction  phase.    During  the  establishment  and  conservation  phases, 


diversity  generally  increases,  though  there  are  cases  where  diversity  declines 
during  these  2  phases.    Diversity  decline  has  been  observed  just  after  the 
point  where  all  species  (both  opportunist  and  mature  life-histories)  occur  at 
the  same  time.   Eventually,  the  opportunist  species  are  out-competed  for 
resources,  and  a  species  decline  is  seen  (Remmert  1991)  where  later 
successional  species  dominate. 

Theories  of  Diversity 

Since  species  diversity  reflects  the  sum  of  additions  and  removals  of 
species  from  an  area,  local  and  regional  processes  that  add  or  subtract  species 
are  important  to  its  understanding.    Species  are  added  to  an  area  by  evolution 
or  by  immigration.  They  are  subtracted  from  an  area  by  total  extinction  or  by 
one  of  following  processes:   change  in  physical  conditions  that  no  longer  meet 
the  needs  of  the  species,  such  as  catastrophes  or  climate  change;  random 
events  that  affect  rates  of  propagation  and  mortality;  or  exclusion  by  local 
processes  due  to  interactions  between  individuals,  such  as  competition, 
predation,  and  pollination   (Hawksworth  and  Kalin-Arroyo  1990). 

Many  theories  have  been  offered  to  explain  differences  in  diversity. 
They  include  both  equilibrium  and  non-equilibrium  processes.  Theories 
invoking  equilibrium  dynamics  suggest  that  the  rate  of  addition  and  removal 
of  species  from  a  region  is  equivalent,  and  that  areas  have  their  own 
equilibrium  value  of  species  diversity.    If  species  are  lost,  the  processes  of 
adding  or  subtracting  species  will  adjust  the  value  so  that  it  remains  constant. 
On  the  local  scale,   equilibrium  theory  suggests  that  species  interactions 
account  for  the  level  of  diversity  saturation,  and  that  local  diversity  is  fed  by 
the  regional  species  pool  (Ricklefs  1990).    Non-equilibrium  or  dynamic 


equilibrium  theory  suggests  that  communities  have  fluctuating  species 
diversity  resulting  from  an  approximate  balance  between  local  and  regional 
processes  (Ricklefs  1990).   Competitive  exclusion  and  its  prevention  are  the 
basis  of  this  theory,  according  to  Huston  (1994). 

The  following  broad  theories  suggest  underlying  processes  of  observed 
patterns  of  diversity:   1)  Diversity  through  fluctuations;  2)  Time  theory;  3) 
Diversity  proportional  to  resources;  4)  Energy  or  empower  diverted  for 
physiological  adaptation;   5)  Population  mechanisms  affecting  diversity;  and 
6)  Habitat  heterogeneity  theory. 
Diversity  through  fluctuations 

More  diversity  is  predicted  when  fluctuations  of  conditions  affecting 
species  is  large  enough  so  that  no  one  species  can  gain  competitive 
dominance,  but  not  so  large  as  to  cause  extinctions.  This  is  the  basis  of  the 
"intermediate  disturbance  hypothesis"  proposed  by  Connell  (1978),  who 
applied  it  to  rainforests.   Rainforests  have  numerous  gaps,  evidence  that 
there  is  almost  perpetual  disturbance  of  small  areas.    Gaps  are  thought  to 
occur  before  competitive  exclusion  within  the  forest  can  reduce  its  diversity. 
Hence,  a  high  diversity  is  maintained  in  rainforests. 

Previously,  Hutchinson  (1961)  questioned  the  mechanism  of  co- 
existence of  many  species  of  phytoplankton  in  relatively  large  lakes  as  well  as 
the  co-existence  of  large  numbers  of  species  in  other  ecosystems.    He  framed 
his  argument  using  the  competitive  exclusion  theory,  and  maintained  that  if 
disturbance  in  the  ecosystem  recurred  more  frequently  than  the  time  over 
which  species  would  have  been  excluded  due  to  interspecific  competition, 
species  richness  should  remain  high. 


Fluctuations  may  cause  extinction  of  rare  species,  leaving  a  lower 
diversity  at  least  temporarily.  Pimm  (1991)  discusses  2  ways  in  which  species 
extinction  rates  increase  for  small  populations:   one  due  to  random 
demographic  causes,  such  as  all  individuals  growing  in  one  year  are  the  same 
sex,  and  the  second  due  to  fluctuations  that  add  to  the  risk  of  random 
extinction.  Supporting  this,  Tilman  and  EI  Haddi  (1992)  reported  a  40% 
extinction  rate  of  species  in  grassland  ecosystems  due  to  drought.  Most  of  the 
species  lost  were  rare  before  the  drought.    Hall  et  al.  (1992)  suggest  that  the 
geographic  distribution  of  individuals  of  a  species  reflects  the  net  energy 
balance  of  the  species,  occurring  in  a  gradient  that  can  be  divided  over  3  parts 
of  the  distribution:   the  range  over  which  the  stored  energy  is  sufficient  for 
reproduction  and  long  term  existence  of  the  species;  the  range  over  which 
stored  energy  is  adequate  for  survival  but  not  reproduction;  and  the  range 
wherein  energy  reserves  must  be  consumed  for  survival.    The  extremes  of 
the  gradient  may  be  the  areas  over  which  a  species  becomes  rare,  and  has  a 
higher  probability  of  extinction  with  fluctuations  that  will  take  energy  for 
adaptation. 

Seasonal  fluctuations  can  increase  diversity  of  an  ecosystem.   For 
example,  in  arid  lands  a  large  annual  flora  may  flourish  during  seasonal 
rainfall  events.   Since  the  perennial  plants  of  arid  lands  are  widely  spaced, 
probably  due  to  competition  for  water,  seasonal  invasion  by  annual  plants  is 
possible  when  water  is  available  (Colin vaux  1986). 
Time  theory 

The  amount  of  time  over  which  a  region  has  been  undisturbed  by 
major  climatic  changes,  such  as  glaciations,  may  influence  the  number  of 
species  that  are  in  the  area.   Tropical  areas  may  have  had  more  time  for 


evolution  to  occur  than  temperate  regions.   Periods  of  drying  may  also  have 
created  refugia  for  tropical  moist  species,  providing  the  isolation  believed 
important  to  speciation.   Temperate  and  arctic  regions  also  have  ancient  floral 
elements,  however,  so  the  time  that  floras  have  had  to  evolve  may  not  be  as 
different  as  originally  proposed  (Ricklefs  1990),  though  temperate  and  arctic 
regions  have  been  periodically  nearly  eliminated  by  expansion  of  glaciers. 

Sanders  (1968)  proposed  the  stability-time  hypothesis,  which  addressed 
the  high  diversity  of  the  deep  sea  benthos  compared  to  continental  shelf 
benthos.  He  suggested  that  the  difference  was  due  to  stability  of  the 
environment  of  the  deep  sea,  allowing  evolutionary  specialization  and  low 
extinction  rates  due  to  environmental  fluctuations.    In  addition,  respiratory 
costs  of  living  in  the  environment  would  be  very  low. 
Diversity  proportional  to  resourrps 

More  diversity  is  predicted  when  there  is  more  energy  to  support  the 
greater  complexity  of  the  system  of  more  species.  Odum  (1960,  1970  and  1971) 
related  diversity  to  the  energy  requirements  using  permutations  and  square 
functions  to  relate  number  of  species  to  energy  requirement.    Where  the 
levels  of  energy  supporting  an  ecosystem  decrease,  species  diversity  has  been 
observed  to  decrease.  For  example:  Energy  theory  was  suggested  by  Connell 
and  Orias  (1964)  but  they  later  retracted  their  view.  Odum  (1970a,  1970b) 
offered  energy  theory  of  diversity  for  a  rainforest  in  Puerto  Rico. 

Wright  (1983)  proposes  a  species-energy  theory,  predicting  diversity  of 
species  from  either  actual  evapotranspiration  (for  plants)   or  total  net  primary 
production  (for  breeding  land  or  freshwater  bird  species)  of  island  biota. 
Wright  showed  that  on  islands  of  varying  sizes  and  locations,  as  the  basic 
energy  resources  decrease,  as  with  increased  latitude  or  seasonal  changes, 


species  diversity  decreases.  Since  lower  populations  of  species  are  supported 
with  less  energy,  there  is  a  higher  probability  of  species  extinctions.  More 
recently  Rosenzweig  and  Abramsky  (1993)  suggested  a  decrease  of  diversity 
with  a  decrease  of  primary  productivity,  as  well  as  a  diversity  decrease  with  an 
increase  of  primary  productivity  above  a  certain  level.   In  other  words,  the 
proposed  pattern  on  a  regional  scale  is  hump-shaped. 

The  relationship  between  energy  and  species  richness  is  scale- 
dependent.   Wright  et  al.  (1993)  show  the  relationship  on  regional  and  local 
scales.   He  emphasizes  the  importance  of  using  extensive  measures  of 
incoming  energy  to  a  system  (total  amount  per  area)  for  the  correlation  rather 
than  intensive  (amount  per  m2). 

The  maximum  empower  theory  predicts  this.    Empower  (definition 
Table  1-1)  is  a  measure  of  the  rate  of  resource  use  by  an  ecosystem  which 
combines  different  resources  in  one  energy-evaluated  unit.   Emergy  is  used  to 
combine  inputs  and  processes  on  a  common  basis  (insolation,  precipitation, 
transpiration,  evaporation,  and  primary  production).  Energy  concentration 
has  been  proposed  as  a  third  measure  of  scale  (with  time  and  space)  for 
relation  to  diversity. 
Population  mprhanisms  affpct  divprsity 

Population  approaches  to  diversity  describe  processes  occurring  at  the 
scale  of  individual  organisms  and  their  interactions  with  and  effects  on  each 
other.    These  include  predation,  competition,  and  mutualism.    The 
implication  is  that  population  densities  change  by  each  individual 
responding  to  its  local  conditions  rather  than  in  response  to  average 
conditions  across  an  area  (Tilman  et  al.  1997). 


Three  approaches  to  understanding  population  dynamics  are  discussed 
by  Tilman  et  al.  (1997):    metapopulation-like  models,  cellular-automaton-like 
models,  and  reaction-diffusion  models.   Huston  (1994)  concludes  that 
competition  has  a  direct  effect  on  species  diversity  only  under  a  very  restricted 
set  of  condition. 
Spatial  heterogpneity  theory 

Large-scale  equilibrium  dynamics  are  implied  with  the  spatial 
heterogeneity  theory  (Huston  1994).   Higher  diversity  is  maintained  by 
having  enough  different,  local  patches  within  a  region  so  that  there  is  a 
balance  of  losses  and  additions.  Patches  that  are  caused  by  disturbance  will 
support  species  of  different  successional  stages.  Patches  that  are  caused  by 
environmental  heterogeneity  will  support  different  species  than  the 
surrounding  environment  due  to  availability  of  different  resources  or  the 
rate  at  which  the  resources  can  be  exploited. 

Diversity  DerlinP 

Insight  into  the  way  diversity  may  be  limited  comes  from  the  study  of 
situations  with  declining  diversity.   In  studies  of  contained  microcosms 
started  with  high  diversity,  an  initial  decline  in  diversity  is  generally 
observed  (Beyers  and  Odum  1993).  For  any  system,  diversity  decline  may 
result  from  processes  within  the  system  or  from  those  outside  of  the  system. 

A  survey  of  common  knowledge  and  published  examples  suggested 
the  following  circumstances  where  species  diversity  has  been  observed  to 
decrease  in  natural  and  experimental  systems. 


Catastrophic  destruction 

Whether  by  human  hands  or  by  large  scale  planetary  events, 
destruction  of  part  of  an  ecosystem  reduces  the  number  of  species. 
Externally  driven  diurnal  and  seasonal    changes 

Declines  in  species  are  observed  as  part  of  the  ups  and  downs  of  abiotic 
influences  such  as  climatic  cycles.   For  example:   The  Sonoran  Desert  in 
Arizona  has  a  rich  winter  ephemeral  flora  appearing  only  in  the  years  when 
the  rainfall  regime  meets  certain  criteria.   When  the  dry  season  begins, 
ephemeral  species  die  and  the  local  diversity  decreases.   Animal  migrations 
also  follow  seasonal  patterns  of  resource  abundance,  with  a  decline  in  local 
diversity  occurring  on  a  yearly  cycle. 
Oscillations 

Declines  in  species  are  observed  during  parts  of  the  oscillatory  cycles 
due  to  internal  rhythms.  Diversity  has  been  observed  to  rise  and  fall.   For 
example:  Ecosystems  pulse  with  recurring  cycles  of  production  and 
consumption  (Odum  1994),  often  due  to  recurring  natural  perturbations  such 
as  forest  fires  or  floods.  In  general,  systems  grow  and  exhibit  successional 
stages  towards  a  climax  stage  which  is  followed  by  a  descent;  this  oscillation  is 
repeated  through  time. 
Concentrated  seeding 

Wherever   species  are  assembled  in  higher  concentrations  than  can  be 
supported,  declines  follow.   For  example:   Experiments  using  microcosms 
showed  that  an  initial  high  seeding  of  species  in  a  microcosm  were  followed 
by  a  decrease  of  species  to  a  lower  level  (Dickerson  and  Robinson  1986).  This 
may  be  the  result  of  competitive  exclusion  due  to  limited  resources  in  a  small 
volume,  as  well  as  due  to  the  increased  maintenance  costs  of  species  as  they 
aged  (Beyers  and  Odum  1993). 
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Competitive  overgrowth  by  pnriched  populations 

Where  conditions  for  growth  become  enriched,  rapid  growth  by  some 
species  can  displace  others.  The  result  is  generally  an  increase  in  the  density 
or  biomass  of  one  or  a  few  species,  reducing  the  available  resources  (such  as 
light)  to  remaining  species  which  may  then  become  locally  extinct. 

Cultural  eutrophy  is  an  example  of  available  resources  increasing 
abruptly.     In  a  five  month  field  plot  experiment,  Kent  (1996)  documented 
species  declines  due  to  the  surge  of  competition  with  the  abundant  growth  of 
cattails,  water  hyacinths,  and  other  species  that  could  most  effectively  use  the 
excess  resource. 

In  a  study  of  old  field  succession  in  Michigan,  Tilman  (1993)  recorded  a 
decline  of  species  in  plots  that  were  fertilized  with  nitrogen  compared  to 
those  that  were  not  fertilized.  An  increase  in  productivity  of  the  plots  was 
correlated  both  with  an  increased  loss  of  species  and  a  decreased  rate  of 
establishment  of  species,  both  contributing  to  a  reduction  in  the  number  of 
species.  Increased  loss  of  species  was  attributed  to  competitive  interactions, 
and  decreased  establishment  of  new  species. 
Changed  conditions 

Where  environmental  conditions  change  so  that  existing  species  are 
not  adapted,  species  decreases  are  observed  at  least  temporarily.  For  example: 
In  a  Barro  Colorado  Island  study,  a  long  term  drying  trend  showed  a  decline 
and  predicted  extinction  of  16  species,  while  colonizer  species  were  increasing 
(Condit  et  al.  1996a).  A  similar  decline  in  species  was  documented  in  a  prairie 
ecosystem,  where  it  was  attributed  to  drought  (Tilman  and  El  Haddi  1992). 
Decreased  resources 

Also,  diversity  decreases  with  decreasing  production.   Production  is 
strongly  correlated  with  precipitation.   Though  diversity  in  general  decreases 


with  less  rain,  correlations  made  of  diversity  and  precipitation  is  not  as 
compelling  as  diversity  and  evaporation.   Increased  competition  due  to  a 
change  in  resources  may  cause  the  diversity  decline. 
Decreasine  diversity  along  gradients 

The  number  of  woody  plant  species  in  0.1  ha  plots  within  69 
neotropical  forests  was  compared,  and  richness  was  correlated  with  soils  and 
climatic  data  (Clinebell  et  al.  1995).   The  most  important  variables  correlated 
with  number  of  species  were  annual  rainfall  and  rainfall  seasonality. 

Gentry  and  Dodson  (1987)  showed  richness  of  epiphytes  in  western 
Ecuador  and  southern  Central  America  to  decrease  with  decreasing  absolute 
precipitation,  using  data  from  local  floras. 

Wright  (1992)  reports  a  consistent  decrease  in  number  of  species  with 
decreased  rainfall  in  tropical  forests,  emphasizing  the  inverse  relationship  of 
rainfall  and  seasonality  throughout  the  wet  tropics.  The  standard  forest  plots 
that  he  examined  showed  a  fivefold  decrease  in  plant  species  densities  in  the 
Neotropics  in  rainfall  gradient  from  4000  to  1000  mm;  and  sixfold  decrease 
over  a  Ghanan  rainfall  gradient  of  1750  to  750  mm. 

The  number  of  species  declines  polewards  from  the  equator  and 
coincides  with  decreasing  rainfall,  temperature,  evapotranspiration, 
insolation,  and  primary  production.     Whether  these  relationships  are  merely 
correlative  or  causal  has  been  addressed  widely,  but  the  patterns  have  been 
observed  for  many  life  forms.    For  example:  Based  on  74  samples  of  0.1  ha 
lowland  sites,  Gentry  (1988)  reports  an  order  of  magnitude  downsizing  in 
number  vascular  plant  species  between  rich  tropical  forests  and  temperate 


forest. 


Gentry  (1982)  used  total  plant  species  numbers  in  0.1  ha  samples  across 
the  upland  Neotropics,  and  showed  a  strong  correlation  between  diversity 
and  precipitation  in  the  Neotropics.   He  later  discussed  that  this  relationship 
may  be  a  special  case  relevant  to  the  Neotropics  where  there  is  a  strong 
correlation  of  total  annual  rainfall  and  strength  of  dry  season  (Gentry  1988). 
Such  a  relationship  was  not  apparent  in  the  samples  from  tropics  of  Africa. 
In  addition,  he  found  that  the  relationship  at  the  high  end  of  actual 
precipitation  was  nonlinear,  becoming  asymptotic  above  4,000  mm.   This  may 
show  saturation  level  for  species  diversity.  He  also  warns  that  the  asymptote 
may  be  an  artifact  of  sampling  methods  -  species-area  curves  do  not  level  off 
for  the  highest  diversity  sites;  thus  the  area  sampled  at  the  lower  sites  may 
not  reflect  diversity  at  the  highest-diversity  sites.   He  further  strengthens  the 
relationship  between  rain  and  number  of  species  with  data  from  1  ha  plots 
measuring  trees  and  lianas  >  10  cm  diameter.    In  upper  Amazonia,  adjacent 
forest  types  were  shown  to  have  different  species  and  were  growing  on 
different  substrates,  with  little  change  in  species  diversity. 

Overlaying  the  latitudinal  pattern  of  diversity  are  regional  or  local 
declines  in  diversity  with  reduced  rainfall.    This  often  is  associated  with 
effects  of  local  or  regional  topography  influencing  temperature  and  rainfall 
patterns.   In  arid  northern  Chilean  Andes  above  the  Atacama  Desert,  Arroyo 
et  al.  (1988)  report  a  steep  gradient  of  plant  species  diversity  corresponding 
with  the  rainfall  gradient.     Surveys  in  1620  minimum  area  quadrats  were 
made  on  6  transects  1/4  degree  latitude  wide  and  between  18-28  S,  running 
from  the  edge  of  the  Atacama  Desert  (1500-3000m)  to  the  elevational  limit  for 
vascular  plants  (4500-5000  m  depending  on  latitude).   Data  showed  species 
diversity  decline  with  increasing  aridity  on  a  latitudinal  gradient  along  the 


western  side  of  the  Andes  as  well  as  from  east  to  west  across  the  Andes  along 
the  rainfall  gradient. 

Based  on  surveys  in  48  sites  in  the  Neotropics  (lowland  tropics  of 
Middle  and  South  America,  montane  cloud  forests,  and  supra-treeline  of 
Andes),  Duellman  (1988)  found  a  decline  in  diversity  from  wet  to  dry  regions 
and  from  low  to  high  elevations  for  the  entire  anuran  fauna. 

Both  a  decrease  in  structure  and  diversity  occurs  with  the  decreased 
input  of  rainfall  and  temperature.     Evapotranspiration  is  a  function  of  both 
rainfall  and  temperature  directly  related  to  plant  production. 
Decreased  information  exrhange  due  to  fragmentation 

Diversity  has  been  observed  to  decrease  where  ecosystem  areas  are 
fragmented  by  insertion  of  other  uses  of  space,  and  are  isolated  from 
exchanges  of  species  and  seeding  with  adjacent  areas.  Bierregard  et  al.  (1992) 
documented  loss  of  species  in  lowland  tropical  rainforest  where  habitat 
islands  were  created  by  the  removal  of  forest  around  various  size  fragments. 
Diversity  decline  due  to  fragmentation  of  habitat  is  also  reported  on  the 
islands  that  were  created  by  the  formation  of  Gatun  Lake  for  the  opening  of 
the  Panama  Canal  (Leigh  et  al.  1993). 
Onset  of  stable  conditions 

Some  researchers  have  related  species  decline  to  the  onset  of  stable 
conditions  after  times  of  greater  fluctuation.    Whereas  competitive  exclusion 
may  be  pre-empted  due  to  disturbance,  the  lack  of  disturbance  will  allow 
interspecies  competition  to  cause  a  higher  rate  of  extinction. 
Growth  in  si?P 

Where  space  is  limited  during  a  time  of  growth,  increasing  size  of 
dominant  species  by  occupying  more  space,  may  exclude  species.  This  is 


demonstrated  in  a  successional  sequence  of  a  very  limited  area  of  forest  where 
a  dominance  hierarchy  develops,  shading  out  species  in  the  understory. 
Loss  of  controllers 

Where  species  are  dependent  on  control  actions  by  other  species,  loss  of 
the  control  species  may  cause  loss  of  species.  For  example,  plants  dependent 
on  insects  for  pollination  were  lost  when  insects  were  not  available. 
Physiological  demand 

Where  special  conditions  require  energy  to  be  diverted  to  sustain  life, 
diversity  has  been  observed  to  decrease.  For  example:  Increased  levels  of 
toxic  stress  may  result  in  reduced  diversity  accompanied  by  an  increase  in 
biomass.  Pratt  et  al.  (1987)  report  a  decrease  of  microorganism  species  present 
in  microcosms  with  increasing  zinc  concentration. 

A  Conceptual  Model  of  Ecosystem  and  Diversity 

In  order  to  overview  the  principal  parts  and  processes  of  an  ecosystem 
and  the  role  of  diversity,  a  conceptual  model  was  developed  in  Figure  1  -1.  By 
using  the  energy  systems  language,  Table  1-2,  the  main  components  and 
processes  that  appear  to  be  involved  in  determining  the  diversity  are  shown 
not  separately  but  as  part  of  a  connected  whole.  Enough  parts  and  pathways 
of  interaction  are  included  to  consider  the  observations  and  alternative 
theories  discussed  in  the  previous  sections.   The  model  diagram  is  a  way  to 
combine  and  synthesize  the  factors  affecting  diversity  such  as  energy, 
competitors,  nutrient  materials,  outside  pulses,  and  internal  oscillations.    The 
conceptual  model  may  be  used  for  all  scales  over  which  diversity  is  studied. 

A  systems  model  is  a  complex  hypothesis.   To  account  for  diversity 
declines,  a  systems  model  must  have  the  main  parts  and  processes  that  affect 


Table  1-2.  Symbols  for  the  energy  systems  language  used  in  this  dissertation. 
From  Odum  (1983). 

Symbol  Explanation 


^.        Energy  circuit:  A  pathway  whose  flow  is  proportional  to  the 

quantity  in  the  storage  or  source  upstream. 

(      )>        Source:  Outside  source  of  energy  delivering  forces  according  to 
v — '  a  program  controlled  from  outside;  a  forcing  function. 

— C  y=T"        Tank:  A  compartment  of  energy  storage  within  the  system 

v— '    }        storing  a  quantity  as  the  balance  of  inflows  and  outflows;  a  state 
variable. 

r  Heat  sink:  Dispersion  of  potential  energy  into  heat  that 

accompanies  all  real  transformation  processes  and  storages;  loss 
of  potential  energy  from  further  use  by  the  system. 

Interaction:  Interactive  intersection  of  two  pathways  coupled  to 
produce  an  outflow  in  proportion  to  a  function  of  both;control 
action  of  one  flow  on  another;  limiting  factor  action;work  gate. 

Consumer:  Unit  that  transforms  energy  quality,  stores  it,  and 
feeds  it  back  autocatalytically  to  improve  inflow. 

Switching  action:  A  symbol  that  indicates  one  or  more  switching 
actions.  ° 

Producer:  Unit  that  collects  and  transforms  low-qualitv  enerev 
under  control  interactions  of  high-quality  flows. 

Box:  Miscellaneous  symbol  to  use  for  whatever  unit  or  function 
is  labeled. 

Constant-gain  amplifier:  A  unit  that  delivers  an  output  in 
proportion  to  the  input  /  but  is  changed  by  a  constant  factor  as 
long  as  the  energy  source  S  is  sufficient. 
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Jg)  Transaction:  A  unit  that  indicates  a  sale  of  goods  or  services 

<-y~y  S2hdJE&  H1  exchange  for  payment  of  money  (dashed  line). 

|  1  nee  (P)  is  shown  as  an  external  source 
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the  species  numbers  and  are  capable  of  accounting  for  observed  declines. 
Figure  1-1  is  used  here  to  suggest  what  parts  and  processes  are  important  and 
how  they  may  account  for  diversity  decline  in  an  ecosystem.  Of  course, 
whether  the  model  and  its  inherent  hypotheses  are  appropriate  can  only  be 
verified  by  study  of  the  behavior  of  the  real  systems.  In  Table  1-3, 
explanations  for  ways  that  diversity  may  decrease  are  suggested  as 
consequences  of  the  way  the  model  was  formulated. 
External  sources  entering  the  system 

The  box  drawn  in  Figure  1-1  defines  the  boundary  of  a  system. 
Everything  outside  of  the  box  is  not  considered  part  of  the  system.  Inflows 
from  sources  outside  the  system  are  represented  by  circles  outside  of  the  box, 
with  pathways  crossing  the  boundary  into  the  system.    The  external  sources 
operate  on  larger  spatial  scales  than  the  system. 

The  circle  on  the  left  labeled  "energy"  represents  the  energy  sources 
sun,  wind  and  rain  that  enter  the  system.     The  inflow  may  vary  over  time 
scales,  such  as  diurnal,  seasonal,  and  geological  depending  upon  the  scale  of 
the  system  being  studied.  The  total  amount  of  energy  entering  is  in 
proportion  to  the  area  of  the  system  being  studied,  as  represented  by  the 
interaction  between  the  storage  tank  labeled  "area"  and  the  energy  source. 
Thus,  the  total  energy  entering  a  larger  area  will  be  greater  than  that  entering 
a  smaller  area  in  the  same  region.    Some  of  the  energy  entering  the  system  is 
used  by  the  producers,  represented  with  bullet-shaped  characters.  Some  of  the 
energy  entering  the  system  flows  through  the  system  without  being  used,  and 
is  shown  as  a  pathway  with  an  arrow  exiting  the  system.  An  example  is  sun 
that  is  reflected  from  surfaces  in  the  system. 
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Over  a  latitudinal  gradient  from  the  equator  to  the  poles,  the  energy 
from  outside  sources  decreases.    The  reduction  of  energy  to  the  system  causes 
a  decrease  in  all  the  storages  and  flows  in  the  system,  including  the  number  of 
species  per  unit  or  area.  If  the  area  over  which  energy  is  flowing  into  the 
system  is  reduced,  such  as  the  fragmentation  of  a  forest  surrounded  by  clear- 
cutting  for  agriculture,  primary  production  and  number  of  species  decrease 
proportionally. 

External  nutrient  sources  can  enter  the  system  with  rainfall,  run-in, 
wandering  wildlife,  air  pollution  or  additions  such  as  fertilizer  or  sewage.    A 
sudden  increase  of  nutrients  may  result  in  a  rapid  rise  in  productivity  of  the 
system,  in  which  case  one  or  several  plant  species  may  respond  with  rapid 
growth.     Some  species  may  become  locally  extinct  from  competitive 
interactions  such  as  shading,  resulting  in  a  system  with  fewer  species. 

Examples  of  stress  from  external  sources  are  changes  in  water  relations 
such  as  drought,  temperature  extremes,  or  input  of  substances  resulting  in 
toxic  buildup  in  the  soil.  Plant  response  to  adapt  to  stress  may  require  use  of 
energy,  taking  energy  away  from  other  uses.  This  may  result  in  the  diversion 
of  energy  from  biomass  production  and  diversity  support  in  order  to  support 
adaptations  for  survival  of  the  stressful  condition. 

Seeding  from  sources  external  to  the  system  depends  on  dispersal 
methods  and  distance  from  the  seed  source.  As  distance  (d)  from  the  source 
increases,  dispersal  of  propagules  into  the  system  decreases  (1/d2).  Where 
ecosystems  are  separated  in  space  from  similar  ecosystems  that  supply 
propagules,  species  with  seeds  that  do  not  have  a  vector  for  farther  dispersal 
may  be  eliminated  from  seeding  an  area. 

Catastrophes  cause  a  pulse  of  nutrients  to  a  system  and  a  resetting  of 
successional  pathways  to  an  earlier  stage.    Systems  are  often  adapted  to 


certain  frequencies,  amplitudes  and  sizes  of  catastrophes  such  as  hurricanes  or 
fires.    If  one  of  these  3  changes,  the  results  may  be  systems  with  lower 
diversity  than  the  initial  system. 
Internal  processes 

Competition  for  light  is  shown  between  the  opportunistic  plant  species 
and  other  species.  Opportunistic  species  are  those  that  establish  as  early 
successional  species,  or  species  that  are  'weedy'.     The  opportunistic  species 
draw  light  from  the  energy  pathway  first,  before  those  species  that  are  later 
successional  species  indicative  of  a  more  mature  system  (labeled  "other 
species  structure").  As  the  later  successional  species  grow,  they  are  able  to  pull 
more  light  for  production  and  overshade  the  early  successional  species.    Local 
diversity  may  show  an  increase  when  species  from  both  early  and  later 
successional  stages  are  present,  and  a  decrease  when  the  earlier  species 
become  locally  extinct  from  overshading.   As  external  energy  sources  increase, 
primary  production  increases.     The  symbol  showing  the  mature  species  also 
shows  physiology,  which  reflects  the  respiratory  costs  of  physiological 
adjustment  to  changing  conditions  within  the  system.   As  stress  from  the 
outside  increases,  the  energetic  cost  of  plant  adaptation  increases.  This  is 
shown  on  the  diagram  as  an  increase  in  the  constant-gain  amplifier  (triangle), 
which  in  turn  increases  the  extinction  rate  from  the  species  diversity  tank. 

Plant  biomass  is  shown  with  4  outflowing  pathways.  The  first  pathway 
flows  into  the  nutrient  storage.  This  occurs  through  decomposition  of  plant 
material  from  litterfall,  plant  mortalities  including  those  due  to  catastrophe, 
and  exudates.  Another  pathway  shows  biomass  feeding  consumers,  where 
all  consumers  in  the  system  are  aggregated  into  one  consumer  storage.  The 
other  2  pathways  from  biomass  lead  to  the  species  diversity  tank.      Biomass 
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(representing  suitable  habitat)  interacts  with  seeding  on  one  pathway, 
allowing  colonization  of  new  species  when  the  energy  base  can  support  it.  As 
the  rate  of  seeding  decreases,  the  rate  of  colonization  also  decreases.  The 
second  flow  to  diversity  represents  the  habitat  base  supporting  extant  species. 

Pathways  from  the  diversity  tank  are  either  extinction  pathways  or 
pathways  that  feed  back  to  the  primary  production  of  the  system.  The 
quadratic  extinction  pathway  represents  extinction  due  to  competition  with 
other  species.   The  rate  of  extinction  increases  when  opportunist  species 
increase,  such  as  with  increase  due  to  a  pulse  of  nutrients.  It  also  increases 
when  the  stress  from  outside  the  system  increases,  requiring  more  energy  for 
plant  adaptation  to  new  conditions.   The  second  extinction  pathway  is  a  linear 
pathway,  representing  extinction  due  to  causes  other  than  competition. 
Feedback  pathways  show  an  increased  efficiency  of  the  system  due  to  new 
species  or  additional  species  added  to  the  system. 

Nutrients  are  added  to  the  system  either  from  outside  sources  or  from 
cycling  of  biomass  or  through  consumers.  Available  nutrients  decrease  when 
they  are  immobilized  by  plant  uptake,  or  when  they  flow  out  of  the  system  in 
processes  such  as  runoff  or  soil  erosion. 

In  any  real  ecosystem  more  than  one  of  these  mechanisms  affecting 
diversity  may  be  operating.   Simulation  is  used  to  relate  the  complex 
hypothesis  to  the  observations. 

Because  of  its  feedback  amplifier  loops  (autocatalytic  input  designs  and 
material  cycles),  the  diversity  model  is  consistent  with  the  maximum 
empower  concept.  This  principle  may  be  stated  as  follows:  Self-organizing 
systems  prevail  by  developing  designs  that  maximize  emergy  inputs  and 
effective  use. 


The  model  in  Figure  1-1  includes  the  general  premise  that  some  species 
variety  is  necessary  to  maximize  a  system's  processes.   To  fit  the  principle, 
increasing  or  decreasing  diversity  needs  to  be  consistent  with  increasing 
and  /or  sustaining  resource  inputs  or  increasing  functional  efficiencies. 

Carrying  Capacity  of  Systems  for  Diversity 

A  general  question  asked  by  those  responsible  for  maintaining 
diversity  in  public  lands  and  parks  is  whether  there  is  a  carrying  capacity  of  an 
ecosystem  for  diversity.   If  the  maximum  empower  principle  applies, 
ecosystems  sustain  that  diversity  that  promotes  total  system  function. 
Priority  in  self-organization  may  go  for  physiological  adaptation  of  fewer 
species  where  this  is  necessary  to  maximize  system  productivity  and 
efficiency. 

Spatial  Organization 

There  are  extensive  published  observations  on  the  spatial  patterns 
developed  by  ecosystems,  including  the  patterns  of  horizontal  plant 
organization  on  the  forest  floor  and  patterns  of  vertical  structure. 
Explanations  for  observed  distributions  have  to  include  the  role  of  light 
energy  and  the  apparently  general  tendency  for  self-organization  to  form 
hierarchies.    Obviously  changes  of  diversity  are  affected  by  the  seeding,  which 
determines  the  locations  and  survival  of  plants.    In  situations  where  humans 
do  the  planting,  how  does  self-organization  develop  afterwards? 


Hierarchical  Organization 

Like  other  ecosystems,  the  rainforest  develops  components  and 
processes  in  an  energy  hierarchy,  which  is  represented  in  model  diagrams  by 
position  from  left  to  right.   Items  on  the  left  have  more  energy  (low  quality 
energy),   turnover  faster,   and  have  smaller  territory.  Position  of  a 
component  or  pathway  flow  in  an  energy  hierarchy  is  measured  by  its 
transformity.  Transformity  is  the  emergy  per  unit  energy  (Table  1-1).    Emergy 
is  the  energy  of  one  kind  previously  used  directly  and  indirectly  to  make  a 
product  or  service.   Calculating  emergy  and  transformities  of  ecosystem 
components  and  flows  is  a  way  of  determining  the  pattern  of  energy 
hierarchy.     Maps  of  transformity  and  changes  of  transformity  over  time  may 
be  estimated. 

Rainforest  in  Biosphere  2 

Biosphere  2  is  a  1.25  hectare  mesocosm  in  Arizona  designed  to  be 
closed  to  material  exchanges  with  the  environment  outside  of  its  glass 
boundary,  but  open  to  exchanges  of  energy  and  some  types  of  information. 
Areas  within  Biosphere  2  were  established  to  represent  different  biomes  of 
earth,  and  included  aquatic  (ocean  and  estuary)  and  terrestrial  (savannah, 
desert,  and  rainforest)  wilderness  areas  and  a  habitat  for  human  residents  and 
their  agricultural  systems  (Figure  1-2).   Whereas  the  atmospheric  chemistry 
and  water  systems  were  global  phenomena  among  all  of  the  areas  of 
Biosphere  2,  temperatures,  rainfall  regimes,  and  humidities  were  maintained 
locally  within  each  area.    Figure  1-3  shows  the  photosynthetically  active 
radiation,  C02  concentration  and  temperature  for  terrestrial  wilderness  areas 
in  the  Biosphere. 
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Figure  1-3.  Continuous  photosynthetically  active  radiation,  C02,  and 
temperature  in  terrestrial  wilderness  biomes  of  Biosphere  2  during  the 
first  material  closure. 


The  footprint  of  the  tropical  rainforest  biome  within  Biosphere  2  is 
1900  m2.  The  highest  point  of  the  pyramidal  glass  structure  enclosing  it  is  22 
m,  measured  from  the  soil  surface  in  the  lowland  rainforest  (Figure  1-4).   The 
total  volume  of  the  rainforest  is  about  35,000  m3  (Dempster  1993),  nearly  17% 
of  the  total  Biosphere  2  volume.   The  contribution  of  the  rainforest  to  the 
atmospheric  chemistry  of  the  larger  scale  mesocosm  has  been  modeled  by 
Engel  and  Odum  (1999).  Their  simulations  suggest  that  the  rainforest  may  be 
responsible  for  up  to  50%  of  the  whole  Biosphere  2  metabolic  rate. 

The  rainforest  biome  was  designed  to  be  functionally  analogous  to  the 
planetary  rainforest  biome.  It  was  built  and  its  climate  managed  to  emulate 
the  general  structure  and  function  of  a  New  World  tropical  rainforest.  It  is  a 
geographical  island,  disjunct  from  the  larger  system  from  which  its 
components  were  obtained. 
Design  Elements  and  Description 

The  climate  of  the  Biosphere  2  rainforest  was  created  by  control  of  rain, 
temperature,  humidity  and  air  flow  so  as  to  support  plant  species  from 
various  humid  tropical  regions.   The  annual  changes  in  day  length  were 
greater  in  Biosphere  2  than  in  a  planetary  rainforest,  and  temperature  and 
atmospheric  C02  ranges  were  greater  on  both  a  daily  and  annual  basis. 

The  rainforest  soils  were  created  from  a  local  desert  grassland  soil 
combined  with  other  organic  materials.   The  design  objectives  were  to 
produce  a  functional  equivalent  of  tropical  rainforest  soils,  to  allow  soils  to 
develop  in  place,  and  to  provide  a  horticulturally  adequate  substrate.    Over 
time,  the  formation  of  humic  and  fulvic  acids  was  expected  to  reduce  the  soil 
pH,  similar  to  a  planetary  rainforest  (Scarborough  1994).    The  soils  were 
designed  to  be  deep  enough  to  allow  expansion  of  roots  needed  to  stabilize 


Figure  1-4.  Side  view  of  rainforest  biome. 


aboveground  canopy  expansion  and  to  contain  a  storage  of  elements 
necessary  plant  growth. 

The  design  challenge  of  the  Biosphere  2  rainforest  was  to  consider 
diversity  over  a  hundred  year  time  scale,  as  well  as  across  the  scales  of 
population,  community,  and  ecosystem,  including  functional  diversity.    The 
rainforest  was  excessively  species-packed  for  self-organization  over  time 
under  minimal  influence  of  human  management.     Since  no  biotic 
introductions  were  to  be  made  during  the  study  period  after  the  initial  closure 
in  1991,  initial  planting  had  to  include  all  of  the  species  for  succession  and 
mature  development.    In  contrast,  an  unconfined  rainforest  experiences 
continuing  immigration  during  successional  changes.    Therefore,  a  much 
larger  number  of  species  was  planted  than  could  survive,  letting  extinction 
occur  as  a  natural  process  in  the  system.   Though  effort  was  made  through 
pruning  to  reduce  competition  for  light  during  the  first  years  of  operation  to 
ensure  the  survival  of  certain  key  species,  decline  in  the  number  of  plant 
species  was  expected.  It  was  hypothesized  that  the  rainforest  structure  would 
change  over  time  from  that  of  a  recently  cleared  ecosystem  to  a  complex 
primary  forest  (Prance  1991).  The  final  mix  of  plant  species  was  eclectic  with  a 
bias  toward  Neotropical  taxa. 

Eight  separate  habitats  were  initially  delineated  in  the  rainforest 
(Prance  1991).   The  habitats  were  named  lowland  rainforest,  ginger  belt, 
varzea,  cloud  bowl,  surface  aquatic  systems,  bamboo  belt,  mountain  terraces, 
and  cliff  faces  (Figure  1-5). 

The  conceptual  design  for  the  rainforest  began  in  1985;  plant  collections 
were  made  from  1986  to  1991;  planting  began  in  1990;  animal  introductions 
began  in  1991;  and  a  total  survey  and  mapping  of  every  plant  was  completed 
prior  to  the  first  material  closure,  which  lasted  from  September  1991  through 
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Figure  1-5.  Habitats  delineated  within  the  Biosphere  2  rainforest. 


September  1993.  A  total  re-survey  of  plants  in  the  rainforest  was  made  after 
the  2-year  closure  and  again  in  1996. 

Measurements  of  energy  flows  and  storages  made  during  and  after  the 
period  of  closure  included  temperature,  relative  humidity,  rainfall  duration, 
plant  sizes,  light  level,  atmospheric  C02  and  02 ,  and  trace  gas  concentrations. 

An  artificial  ecosystem  such  as  the  Biosphere  2  rainforest  or  a  botanic 
garden  might  be  expected  to  behave  according  to  models  of  a  natural 
ecosystem  with  equivalent  energy  inputs  and  available  information.    What 
are  the  different  options  for  ecosystems  to  maximize  power,  and  what  aspects 
of  human  management  of  'artificial'  systems  will  alter  this? 

A   Conceptual  Model  of  Diversity  in  Rainforest  of  Biosphere  7 

Whereas  Figure  1-1  has  the  main  features  of  any  system  believed  to 
affect  the  diversity  of  stored  information,   Figure  1-6  represents  these  concepts 
for  the  rainforest  biome  area  of  Biosphere  2.   The  diagram  is  a  complex 
hypothesis  for  understanding  diversity  changes  in  Biosphere  2.  The 
differences  between  parts  and  processes  of  the  earlier  conceptual  model  of 
diversity  in  ecosystems  (Figure  1-1)  and  the  model  of  the  Biosphere  2 
rainforest  are  as  follows. 
External  sources  entering  thp  system 

The  striped  border  in  Figure  1-6  defines  the  glass  boundary  of  the 
Biosphere  2  structure  that  physically  separates  it  from  the  outside.     The  inner, 
solid  box  in  the  diagram  represents  the  rainforest  biome,  one  of  7  biomes 
inside.   There  was  no  physical  barrier  separating  the  rainforest  biome  from 
the  other  biomes.   Though  many  of  the  parts  and  processes  within  the 


rainforest  model  are  the  same  as  the  earlier,  general  model,  the  external 
sources  are  quite  different. 

The  primary  energy  source  entering  the  rainforest  was  sunlight;  rain 
and  wind  do  not  enter  from  an  external  source.   Since  the  Biosphere  2 
rainforest  was  contained,  its  area  was  constant  and  the  tank  for  land  area 
interacting  with  energy  has  been  deleted.   During  periods  of  material  closure, 
no  nutrients  entered  the  Biosphere  from  outside.    Certain  stress-causing 
external  conditions  could  have  affected  the  producers,  such  as  a  long  period  of 
reduced  sunlight.   Electricity  and  information  were  external  sources  used  to 
operate  and  design  the  mechanical  systems  providing  the  subsidies  normally 
produced  by  nature,  such  as  waves  in  the  ocean,  airflow  throughout  the 
Biosphere,  and  rainfall  operated  by  a  sprinkler  system. 

The  species  that  were  planted  in  the  rainforest  were  selected  from  the 
earth's  species  pool,  with  human  decisions  rather  than  habitat  and  dispersal 
characteristics  determining  which  plant  species  would  colonize  the  area. 
Catastrophes  from  outside  sources  would  be  those  that  would  have  affected 
the  operating  systems  or  caused  massive  collapse  of  the  external  structure.  A 
total  shutdown  of  the  power  systems,  for  example,  may  have  caused 
overheating  during  the  day,  resulting  in  plant  mortalities.   This  would  have 
resulted  in  a  pulse  of  dead  biomass  that  would  have  been  placed  in  the  dry 
biomass  storage,  shown  with  the  interaction  of  catastrophe  and  biomass 
pathways,  flowing  into  the  dry  biomass  storage  tank. 
Internal  processes 

The  cycling  of  nutrients  between  the  rainforest  and  other  areas  inside 
the  Biosphere  included  flows  driven  by  the  operating  equipment.  Airborne 
nutrients  such  as  C02  were  exchanged  by  airflow  through  the  air  handlers. 


Soil  borne  nutrients  exited  the  rainforest  in  water  that  percolated  through  the 
soil  profile  and  was  pumped  to  a  central  storage  reservoir,  where  it  was  mixed 
with  water  from  the  other  areas.  Some  of  these  nutrients  were  then  pumped 
back  to  the  rainforest  in  rain  water. 

Biospherians  lived  in  their  habitat  area  outside  of  the  rainforest,  but 
visited  the  rainforest  regularly.   The  model  shows  an  interaction  of  the 
Biospherian  pathway  with  biomass,  representing  pruning  of  the  weedy 
biomass.  Much  of  the  weedy  biomass  was  dried  and  stored,  shown  with  the 
pathway  to  the  dry  biomass  tank. 

Some  of  the  rainforest  consumers  traveled  between  areas  in  the 
Biosphere.   This  is  shown  with  a  flow  pathway  crossing  over  and  then 
returning  to  the  rainforest  system.     Examples  were  the  galagos  that 
concentrated  their  activities  in  the  rainforest  but  traveled  throughout  the 
wilderness  biomes. 

Though  the  process  of  colonization  of  the  rainforest  by  plant  species 
differed  from  wild  systems,  the  extinction  process  occurred  similarly.    This  is 
shown  in  both  conceptual  models  with  a  quadratic  function  representing 
competition  with  other  plant  species  and  a  linear  function  representing  other 
reasons  for  extinction. 

Rainforest  in  Pnprtp  Rirn 

The  rainforest  compared  with  that  in  Biosphere  2  is  in  the  Luquillo 
Experimental  Forest  (USDA)  of  northeastern  Puerto  Rico.   The  area  has  4 
different  vegetation  zones  -  tabonuco  forest,  Colorado  forest,  sierra  palm 
brake,  and  elfin  forest.   The  tabonuco  rainforest,  named  after  its  dominant 
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tree  species,  has  been  the  most  extensively  studied  (Reagan  and  Waide  1996) 
and  is  the  type  used  for  comparison  in  this  study. 

Long-term  data  are  available  for  the  El  Verde  study  site  and  the  Bisley 
watersheds  and  MAB  biodiversity  plots.  An  extensive  study  of  the  tabonuco 
forest  type  funded  by  the  Atomic  Energy  Commission  was  undertaken  from 
1963-1968,  before  and  after  irradiation  (Odum  and  Pigeon  1970).  Subsequently, 
a  volume  was  published  (Reagan  and  Waide  1996)  describing  the  food  web  of 
EI  Verde. 

Among  the  many  papers  that  deal  with  plant  diversity  and 
successional  patterns  in  the  tabonuco  forest  are  the  following:    Lugo  (1992) 
compared  the  plant  species  in  pairs  of  tree  plantations  and  secondary  forests; 
Taylor  et  al.  (1996)  analyzed  numbers  and  densities  of  species  by  life  form  in 
an  irradiated  area  of  the  tabonuco  forest  over  23  years  after  irradiation  with 
data  from  8  surveys;  landslide  succession  has  been  documented  by  Myster  and 
Walker  (1997),  Walker  et  al.  (1996),  and  Guariguata  (1990);  Crow  (1980)  studied 
species  changes  in  a  0.72  ha  plot  over  30  years,  and  Johnston  (1990)  followed 
Crow's  work  with  a  dissertation  documenting  successional  change  in  the 
same  plots. 

Recent  papers  describing  the  effects  of  Hurricane  Hugo  on  the  tabonuco 
forest  system  and  the  re-development  of  the  system  to  pre-hurricane 
conditions  are  published  in  a  special  issue  of  the  journal  Biotropica,  Volume 
28(4a),  1996. 


CHAPTER  2 
METHODS 


In  this  study,  measurements  of  diversity  were  related  to  measures  of 
the  resources  and  productivity  in  Biosphere  2,  in  the  rainforest  in  Puerto  Rico 
and  other  situations  where  limited  resources  may  have  restricted  diversity. 

Construction  and  Operation  of  the  Rainforest  in  Biosphere  7 

A  chronology  of  the  construction  and  operation  of  the  mechanical, 
structural  and  life  systems  of  the  Biosphere  2  rainforest  from  1987  through 
1999  is  shown  in  Table  2-1.   Three  different  modes  of  operation  of  Biosphere  2 
are  distinguished  by  degree  of  material  closure  from  the  outside  or  the 
separation  between  biomes  on  the  inside  of  the  structure,  and  degree  of 
control  over  composition  of  the  atmosphere,  as  follows:   (1)   Complete 
material  closure,  with  the  entire  atmosphere  inside  Biosphere  2  interacting  as 
a  unit  separate  from  the  outside  atmosphere.   (2)  Forced  airflow  through  the 
entire  Biosphere  from  the  outside.   (3)  Physical  separation  of  the  wilderness 
systems  from  the  habitat  and  agriculture  systems;  separation  of  the 
wilderness  biomes  from  each  other  with  curtains;  and  forced  airflow  from 
outside.   Construction  and  engineering  details  of  Biosphere  2  are  given  in 
Zabel  et  al.  (1999). 
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Procedures  Used  in  Starting  the  Ecosystem 

Soil  Components  and  Placement 

The  soil  profile  was  assembled  with  subsoil  and  topsoil  layers.  The 
subsoil  was  a  mix  of  a  coarse  sandy  loam  (highly  weathered  granitic  grus) 
extracted  from  a  local  quarry  and  5-15  cm  granite  rocks.  It  had  uniform 
composition  throughout  the  rainforest,  and  its  thickness  varied   from  0  to 
about  5  m  depending  upon  specifications  made  for  root  growth  and  the 
architectural  structure  over  which  it  was  placed.   Criteria  for  subsoil  selection 
were  good  downward  percolation  of  water  after  repeated  wettings;  radon 
emissions  not  greater  than  background  crustal  emissions;  low  heavy  metal 
content;  price-competitiveness;  and  proximity  to  the  Biosphere  2  site 
(Scarborough  1994).  The  subsoil  volume  was  estimated  to  be  3340  cubic 
meters,  based  on  measured  volumes  of  the  structure  into  which  it  was  placed. 
Table  2-2  lists  the  quarry  locations  of  the  materials  extracted  for  the  rainforest 
soils.   The  subsoil  was  placed  over  a  shallow  gravel  layer  which  covered  and 
surrounded  drainage  pipes  that  carried  water  percolated  through  the  soils  to  a 
collection  trough,  from  which  it  was  pumped  to  a  global  storage  system  inside 
the  Biosphere. 

Topsoils  were  placed  on  top  of  the  subsoil.  Their  thickness  varied  from 
0.3  -  3.2  m,  averaging  0.9  m.  Four  mixtures  formed  the  major  volume  (1750 
cubic  meters)  of  topsoil  in  the  rainforest,  with  small  amounts  of  3  additional 
mixes  (16  cubic  meters)  used  for  more  specialized  habitats  (Scarborough  1994). 
Each  mixture  consisted  of  a  different  combination  of  the  following  materials: 
a  local  desert  grassland  soil  ('Wilson  Pond  soil'),  organic  matter,  gravelly 


Table  2-2.  Rainforest  soil  components  and  their  origins.  From 
Scarborough  (1994). 


Material 

Coarse  sandy  loam  (subsoil) 

Wilson  Pond  soil 


Gravelly  sand 
Compost,  Coarse  organic 

material 
Fine  peat 
Coarse  peat 

Canadian  sphagnum  moss 
Pumice  chunks 


Origin  ~ 

Kalamazoo  sand  and  gravel  quarry,  just 

southwest  of  San  Manuel  Copper  Mine, 

Mammoth,  AZ 

Along  boundary  between  Sees  2  &  3,  T10S, 

R14E,  1200  feet  south  of  highway  82  and  1.3 

miles  west  of  the  Biosphere  Road 

Quarried  on  location 

AAA  Fertilizer,  Tucson,  AZ;  various  other 

local  suppliers  of  organic  materials 

Purchased  commercially 

Purchased  commercially 

Purchased  commercially 

Purchased  commercially 
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sand,  sandy  loam,  and  pumice  chunks.   Table  2-3  shows  the  different 
combinations  and  volumes  of  materials  that  were  used  and  their  placement 
within  the  ecosystem.   The  'Wilson  Pond  soil'  was  an  organic-rich  (4-5%  OM) 
silt  loam  to  clay  loam  which  was  quarried  from  a  large  earthen  cattle  tank. 
The  organic  amendments  were  compost  (made  from  forest  mulch,  cotton  gin 
trash,  rotted  cattle  silage,  alfalfa  hay,  and  cattle  manure);  coarse  peat;  coarse 
organic  material;  and  Canadian  sphagnum  moss.   The  gravelly  sand  was 
quarried  on  location.    Additional  materials  used   in  much  smaller  amounts 
were  the  coarse  sandy  loam  subsoil,  pumice  chunks,  and  fine  peat. 

The  soil  materials  were  either  mixed  with  heavy  equipment  outside  or 
directly  on  conveyor  belts  feeding  from  the  outside  to  the  inside  of  the 
Biosphere.   Various  types  of  heavy  equipment  (both  track  and  wheel 
vehicles),  conveyor  belts  and  /or  shovels  were  used  for  soil  placement  inside 
the  rainforest.   Leaf  litter  and  humus  from  local  natural  ecosystems  were 
added  to  the  soil  surface  after  placement,  and  thousands  of  purchased 
earthworms  and  a  much  smaller  number  of  earthworm  species  that  were 
field-collected  in  southeastern  Texas  (Scott  1994)  were  placed  in  the  soil. 

Collection  and  Initial  Placement  of  Plant* 

Starting  in  1986,  plants,  seeds,  and  other  propagules  assembled  for  the 
rainforest  were  grown  in  greenhouses  located  near  Biosphere  2.     Plant 
accessions  were  acquired  from  botanic  gardens,  plant  nurseries,  and  private 
collections,  and  from  field  collections  in  Puerto  Rico,  Belize,  Venezuela,  and 
Brazil.   Permits  were  acquired  for  exportation  from  countries  of  origin  and  for 
importation  to  the  United  States  and  Arizona;   state  and  federal  quarantine 
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and  inspection  protocols  were  followed.     Following  soils  placement  in  the 
rainforest  biome,  plant  accessions  were  transplanted  from  greenhouses  into 
Biosphere  2.  The  largest  plants  were  placed  first  due  to  the  logistics  of 
lowering  them  into  the  Biosphere  with  cranes  before  the  glass  was  installed. 
Plants  were  fertilized  and  watered  in  place  with  drip  irrigation,  sprinklers,  or 
by  hand  with  hoses. 

Each  plant  was  mapped  by  standard  survey  methods  (Thompson  1992) 
using  different  symbols  to  represent  plant  growth  forms.   Measurements  such 
as  basal  diameter,  crown  width,  and  stem  length  were  made.  Each  plant  was 
assigned  a  unique  number,  and  data  on  location,  geographic  origin,   size,  and 
phenology  were  recorded.   Thus  a  history  of  every  plant  in  the  rainforest 
biome  could  be  tracked  through  time. 

Habitat  Asspmhly 

The  initial  design  (Prance  1991)  created  8  habitats  within  the  rainforest 
biome,  each  planted  with  a  different  assemblage  of  plants  including  canopy 
trees,  ground  cover,  shrubs  and  intermediate  level  life  forms,  and  epiphytes 
and  vines.    Initial  plantings  of  fast  colonizers  and  secondary  forest  species 
(Clitoria  racemosa,  Carica  spp.,  Leucaena  spp.,  Cecropia  schreberiana) 
provided  shade  and  structure  to  the  emerging  and  more  characteristic 
primary  rainforest  species  (Prance  1991). 
Lowland  rainfnrP^t 

The  concave  topography  of  the  lowland  rainforest  area  centered  in  the 
southeast  quadrant  of  the  rainforest  allowed  trees  to  reach  a  height  sufficient 
for  a  layered  canopy.  The  topsoil  mixed  for  the  lowland  rainforest  was  50% 


loam,  25%  gravelly  sand,  and  25%  coarse  organic  material.  There  were  lianas, 
epiphytes,  and  many  broad-leaved  trees  and  shrubs, 
ginger  belt 

The  ginger  belt  surrounded  the  rainforest  on  all  sides  bounded  by  glass, 
excluding  only  the  south  edge  which  abutted  the  top  of  the  beach  cliff  and  the 
upper  savanna.  Its  purpose  was  to  shield  the  forest  understory  from  excessive 
light.   Studies  on  isolated  patches  of  planetary  rainforests  show  that  the 
typical  understory  vegetation  is  altered  for  some  distance  beyond  the  forest 
edge  in  part  by  light  penetration  (Bierregaard  et  al.  1992).  Thus  this  peripheral 
dense  belt  of  vegetation  was  planted  to  filter  the  light  and  allow  the  interior 
of  the  forest  to  develop  the  shaded  conditions  that  would  foster  an  understory 
more  typical  of  extensive  rainforests. 

Ginger  belt  topsoils  were  a  mix  of  80%  loam  and  20%  compost.   Soil 
depth  was  about  60  -  90  cm  on  the  northwest  and  west  side  where  no  subsoil 
underlies  the  topsoil,  and  1  m  on  the  northeast,  east  and  southeast  sides, 
where  it  was  underlain  by  a  subsoil  of  variable  thickness.     The  ginger  belt 
ranged  from  about  1  m  to  about  4  m  wide.  Plants  in  the  Order  Zingiberales 
dominated  the  ginger  belt.   Plants  in  these  genera  were  most  abundant: 
Musa,  Heliconia,  Alpinia,  Strelitzia,  and  Costus. 


Varzea 


The  varzea  habitat  was  designed  to  resemble  a  forest  that  is  seasonally 
flooded.   It  featured  a  tightly  meandering  stream  that  ran  from  the  pond  to 
the  edge  of  the  savanna  biome.   The  stream  course  was  made  by  first  filling 
the  regular  topsoil  layers  to  the  specified  level;   then  excavating  the  stream 
courses;  and  then  lining  them  with  concrete  and  thick  PVC.     The  varzea 
topsoil  was  60%  loam  and  40%  coarse  peat.    The  varzea  was  not  flooded 


during  the  study  period.    It  was  planted  with  Phytolacca  dioica,  Pachira 
aquatica,    Pterocarpus  indicus,  and  palms. 
Cloud  bowl 

Suggested  by  the  tepui  sandstone  formation  from  eastern  Venezuela, 
the  central  mountain  'cloud  bowl'  was  planned  as  a  cooler  and  more  humid 
microclimate.  The  soil  was  less  than  1  m  deep.  It  was  a  mix  of  40%  Canadian 
sphagnum  moss,  50%  coarse  peat  moss,  and  10%  loam.   Soil  in  the  planter 
pockets  extended  all  the  way  to  the  bottom  of  the  mountain.    Bryophytes, 
carnivorous  plants,  shrubs,  gingers  and  aquatic  plants  were  introduced  to  this 
habitat. 
Surface  aquatic  habitat* 

Surface  aquatic  habitats  in  the  rainforest  biome  were  hydrologically 
connected.    The  moss  seep  and  upper  pond  in  the  cloud  bowl  spilled  over  a 
cut  in  the  edge  of  the  cloud  bowl,  creating  a  waterfall.  The  waterfall  poured 
into  a  splash  pool,  which  overflowed  into  the  larger,  lower  pond  and  then 
into  the  varzea  stream.   At  the  bottom  of  the  stream  water  flowed  over  a  weir 
into  a  sump,  from  which  the  water  was  pumped  back  to  the  lower  pond. 
Another  pump  moved  water  back  to  the  cloud  bowl.  The  water  systems  were 
underlain  with  a  thick  vinyl  liner  beneath  a  concrete  layer.   Taxa  inhabiting 
the  aquatic  habitats  at  the  beginning  of  the  first  material  closure  were  Azolla, 
Eichhornia,    Nymphaea,  and  Typha. 
Bamboo  hplt 

A  bamboo  belt  was  constructed  along  the  south  edge  of  the  rainforest  to 
baffle  any  airborne  salt  particles  from  the  forest  interior.  Soils  of  the  bamboo 
belt  were  50%  loam,  25%  gravelly  sand,  and  25%  coarse  organic  material, 


varying  between  30  -  90  cm  deep  and  laying  over  air  delivery  plenums.   There 
was  no  subsoil  beneath  the  topsoil  in  the  bamboo  belt.    Bambusa  multiplex,  B. 
tuldoides  and  other  species  of  bamboo  initially  formed  the  major  structure  of 
this  habitat. 
Mountain  terraces 

The  mountain  terraces  skirted  the  rainforest  mountain  on  the  east, 
north,  and  west  sides,  extending  to  the  west  side  of  the  lower  pond.   Flagstone 
walls  were  in  place  around  the  periphery  of  the  terraces,  separating  them 
from  the  ginger  belt  habitat  on  the  west  and  north  sides  and  from  the  lowland 
rainforest  habitat  on  the  east.   The  flagstone  was  placed  during  construction  to 
prevent  soil  from  eroding  down  the  somewhat  steep  grade  from  the 
mountain  to  the  ginger  belt.  The  topsoil  mix  for  the  mountain  terraces  was 
50%  loam,  25%  gravelly  sand,  and  25%  coarse  organic  material.  The  soils  to 
the  west  and  north  of  the  mountain  were  about  0.3  -  3  m  deep,  and  those  east 
of  the  mountain  were  about  0.3  -  6  m  deep.    The  plants  on  the  terraces 
initially  included  Carica  papaya,  Clitoria  racemosa,  Cojfea  arabica, 
Carludovica  palmata,  Inga  sp.,  Hibiscus  rosa-sinensis,  and  Manihot  esculenta. 

The  cliff  faces  of  the  mountain  had  pockets  of  soil  to  support  vines, 
bromeliads,  and  other  plants  that  were  to  cloak  its  surface.    Vines  from  the 
Araceae,  Passifloraceae,  and  Vitaceae  family  and  ferns  planted  in  the  cliff  face 
planter  pockets.  Rhyolite  pumice  chunks  were  blended  at  10  to  15%  of  the 
total  soil  volumes.   The  high  porosity  of  the  pumice  provides  water-holding 
capacity  and  releases  trace  nutrients  as  it  weathers. 


Climate  Maintenance  Systems 

Water  system 

In  contrast  to  the  desert  and  savanna  biomes,  the  rainforest  plants  were 
active  continuously;  dormancy  was  not  part  of  the  yearly  cycle.    In  the 
tropical  rainforest  biome,  water  applied  as  rainfall  flowed  through  and  exited 
the  soil  profile  throughout  the  year.   Figure  2-1  shows  the  water  flows  and 
reservoirs  of  the  rainforest  biome  and  their  connection  with  the  global 
system. 

Rain  was  distributed  to  much  of  the  rainforest  biome  through 
overhead  sprinklers  mounted  in  the  space  frame;  other  areas  were  irrigated 
with  ground  sprinklers  or  a  drip  irrigation  system.    Water  vapor  evaporated 
from  soil,  rain,  or  surface  waters,  transpired  from  plants,   or  delivered 
through  a  misting  system  could  subsequently  have  been  condensed  from  the 
air  using  the  air  handlers  located  in  the  rainforest  basement.    Water 
condensate  that  collected  on  the  inside  surface  of  the  windows  during  late 
autumn,  winter  and  early  spring  was  an  additional,  seasonal  source  of 
condensate.   Water  from  both  of  these  sources  was  re-used  in  the  wilderness 
biomes,  and  was  one  of  the  sources  of  rainwater. 

The  major  water  reservoirs  in  the  rainforest  were  atmosphere 
(humidity),  soil,  water  storage  tanks,  the  surface  aquatic  habitats,  and  plants 
and  other  biota.     The  major  water  flows  among  reservoirs  in  the  rainforest 
were  rainfall  and  irrigation,  subsoil  drainage,  reverse  osmosis  system  flow, 
condensation,  mist,  evaporation,  root  uptake,  transpiration,  and  diffusion  as 
water  vapor  to  other  areas. 
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Figure  2-1.  Hydrologic  connection  of  Biosphere  2  rainforest  to  the  rest  of  the 
Biosphere  2  system. 
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Some  of  the  rain  and  irrigation  water  percolated  through  the  entire  soil 
profile.  This  sub-soil  water,  carrying  substances  leached  from  the  soils,  was 
collected  and  stored  in  a  storage  tank  for  future  use.    The  subsoil  water  in  the 
Biosphere  was  reused  either  as-is  or  after  removal  of  dissolved  solids  by  a 
reverse  osmosis  system.    The  remaining  rain  and  irrigation  water  was  either 
held  in  the  soil  pores  by  matric  forces,  taken  up  by  plants  and  retained,  or 
diffused  back  into  the  atmosphere  via  evaporation  and  transpiration.   The 
average  monthly  relative  humidity  was  over  65%    year-round  in  the 
rainforest  biome,  with  daily  minima  above  50%. 
Air  handling,  temperature,  and  humidity  control  systpm 

Seven  air  handlers  in  the  rainforest  basement  were  operated  to 
regulate  temperatures,  extract  condensate  from  the  air,  and  to  create  air 
movement.     Temperature  was  controlled  by  a  heat  exchange  between  a  water 
coil  and  circulating  air  in  the  handler.   Three  temperature  classes  of  water 
originating  externally  to  the  Biosphere  were  circulated  in  the  coils  of  the 
system:  heating  water,  cooling  tower  water,  and  chilled  water  which  reached 
the  lowest  temperatures.  The  rate  of  the  air  flow  through  an  air  handler  was 
controlled  by  opening  or  closing  an  'econodisc'  located  inside  each  of  the  air 
handlers.   Though  normally  controlled  remotely,  they  were  also  controlled 
manually  from  inside  the  air  handler.    Air  from  the  air  handlers  circulated 
through  a  system  of  plenums,  ducts  and  openings  from  the  basement  into  the 
rainforest.     The  air  returned  through  gratings  on  the  west  and  northwest 
periphery  of  the  rainforest,  where  it  entered  the  basement  and  was  again 
pulled  through  the  air  handlers  and  over  the  heating  coils. 


Human   Intervention 

During  the  2  year  closure  of  the  system  (1991-1993)  certain  herbaceous 
plants  in  the  rainforest  biome  were  extensively  pruned  for  two  reasons.   The 
first  reason  was  to  arrest  primary  succession.    Early  successional  species  were 
pruned  so  that  the  later  successional  species  would  survive.   Pruning  of  vines 
occurred  when  they  appeared  to  be  heavily  shading  the  trees  and  understory 
and  when  vines  had  grown  into  the  tree  canopies.   Both  vines  and  ginger  belt 
herbs  were  pruned  when  they  pressed  against  the  glass  causing  algal  growth 
on  the  glass  surface.   These  selective  harvests  decreased  the  competitive 
advantage  of  high  net-producers  (e.g.,  Ipomoea  and  Passiflora  vines)  over 
species  important  for  the  long-term  structure  of  the  rainforest,  particularly 
the  larger  trees.   The  second  reason  for  harvest  was  an  attempt  to  increase 
sequestration  of  C02  and  to  increase  02  production.  High  levels  of  C02  were 
thought  to  be  decreasing  the  ocean  pH,  and  oxygen  concentrations  eventually 
decreased  to  the  point  of  affecting  the  health  of  the  human  inhabitants.   To 
meet  these  goals,  plants  were  propagated  in  an  attempt  to  cover  vertical 
surfaces  with  photosynthetic  biomass.   Growth  of  high  net-producing  species 
(Ipomoea  sp.,  Passiflora  edulis)  and  herbaceous  plant  species  in  the 
circumferential  ginger  belt  habitat  was  encouraged  by  judicious  pruning  in 
areas  where  they  would  not  overwhelm  other  species,  given  the  constraint  of 
diversity  maintenance. 

The  pruned  biomass  was  removed  from  the  rainforest,  dried  to  retard 
respiration,  and  stacked  in  the  basement.  A  small  percentage  of  the  material 
was  used  for  fodder  for  domestic  livestock.  Of  the  three  terrestrial  wilderness 
biomes,  the  rainforest  required  the  most  time  to  manage  for  regulation  of  the 
atmosphere,  while  seeking  to  maintain  species  richness.   Estimates  of  the 
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amount  of  biomass  removed  through  pruning  were  made  from  weekly 
records  of  time  used  for  each  task  and  from  field  journals. 

Plant  Mapping  and  Identifications 

Initial  Mapping  and  Identifications:    1990-1991 

By  September  1991,  over  1800  individual  plants  had  been  planted  in  the 
rainforest  mesocosm.  Every  plant  was  marked  with  a  pink  flag  that  had  a 
unique  number  written  on  it  which  tracked  it  in  the  database,  linking  its 
location  in  the  rainforest  with  size  measurements  and  other  information, 
including  its  origin,  when  available.   Thus  a  history  of  every  plant  in  the 
Biosphere  2  rainforest  could  be  tracked  through  time.   Specimens  of  many 
accessions  are  maintained  in  herbarium  cabinets  on  location  at  the  Biosphere 
2  Center,  and  some  were  placed  in  the  herbarium  of  the  New  York  Botanical 
Gardens.   Some  plants  have  not  been  identified,  and  many  identifications 
that  were  made  have  not  been  verified. 

Species  Additions  and  Removals 

During  the  first  2-year  closure,  Sept.  1991-Sept.  1993,  additional  papayas, 
bananas,  malanga,  and  canna  were  planted  in  the  rainforest  for  food 
production.  These  were  mapped  and  logged  into  the  database  after  the  first 
closure.   Additionally,  some  plants  that  had  been  planted  before  the  first 
closure  had  been  omitted  from  the  original  survey,  and  the  corrections  were 
made. 


During  the  beginning  of  the  period  (Nov.  1993-Dec.  1993)  between  the 
first  and  second  closures  an  additional  339  individuals  were  planted  in  the 
rainforest.   The  Leucaena  trees  that  had  been  planted  for  initial  shade  of  the 
key  rainforest  species  were  removed  and  weighed  for  subsequent  estimates  of 
biomass. 

Field  Measurements 

The  first  survey  was  the  record  made  at  the  time  of  planting,  shown  in 
Figure  A-l. 

Second  Survey:   Transition  Period.  1993-1994 

Base  maps  used  for  the  second  survey  were  made  from  the  original 
CAD-generated  survey  maps  completed  in  1991.  The  complete  map  of  the 
rainforest  was  divided  into  cells,  each  covering  approximately  10.6  m  x  7  m  at 
a  scale  of  1:24.  Thirty-five  complete  or  partial  cells  provided  total  coverage  of 
the  rainforest. 

Each  plant  found  from  the  original  survey  was  marked  with  a  black 
plastic  tag  (approx.  5  by  8  cm)  with  its  unique  number  etched  in  white.  The 
tags  were  fixed  around  the  base  of  plants  with  plastic  cable  ties.    Data  forms 
that  were  computer-generated  for  each  plant  specified  the  measurements  that 
were  to  be  taken  for  the  plant.   Data  recorded  on  the  forms  were  entered  into  a 
data  base,  and  the  paper  forms  stored  at  the  Biosphere  2  Center.  Plants  that 
had  not  been  mapped  before  the  1991  closure  were  tagged  and  added  to  the 
original  maps.   Additionally,  plants  that  were  introduced  during  the  period  of 
the  survey  were  tagged  and  noted  on  the  map. 


Third  Survey:    Tune-August  1996 

The  third  plant  survey  was  made  in  1996  using  blueline  base  maps 
copied  from  the  original  maps  that  had  been  modified  during  the  1993-1994 
resurvey.   We  searched  for  every  plant  indicated  on  the  base  maps  and  for  all 
plants  that  were  listed  in  the  database  provided  by  Biosphere  2  Center.  Each  of 
the  35  cells  had  a  form  listing  all  of  the  species  within  the  cell.  We  assigned 
one  of  three  status  categories  to  each  plant:   alive,  dead  or  not  found.  When  a 
tag  was  located  with  no  plant  next  to  it  at  the  mapped  location,  the  plant  was 
classified  as  'dead'  and  the  tag  was  removed.  When  neither  tag  nor  plant  was 
found  at  the  mapped  location,  the  plant  was  considered  'not  found'.   In  most 
cases,  the  'not  found'  category  plants  were  dead,  though  since  we  had  neither 
evidence  of  correct  location  by  finding  the  tag  nor  a  live  plant  we  created  this 
category  to  indicate  uncertainty  (plant  tags  are  easily  lost  in  or  were 
previously  removed  from  this  site).     When  a  living  plant  of  the  same  species 
listed  in  the  database  was  found  in  the  location  of  the  survey  point,  it  was 
classified  'alives'. 

When  the  identification  listed  on  the  forms  was  incorrect  or  suspect,  it 
was  corrected  or  questioned  on  the  field  forms.  On  less  than  10  occasions,  a 
plant  classified  as  'dead'  or  'not  found'  in  the  previous  survey  was  found 
alive  in  the  1996  survey.   When  we  were  unable  to  locate  the  tag  on  a  plant, 
we  tied  survey  flagging  tape  to  it  with  the  original  survey  number  and  Latin 
binomial  written  in  waterproof  ink  so  that  replacement  tags  could  be  made 
and  attached  at  a  later  time.   Collections  of  plants  in  flower  were  made  and 
pressed  for  later  identification. 


New  numbers  were  assigned  to  plants  that  were  in  one  of  the 
following  categories:  plants  (excluding  tiny  seedlings  with  only  cotyledons  or 
the  first  few  leaves)  recruited  from  seeds  produced  inside  the  Biosphere  (e.g. 
Pachira  aquatica);  plants  that  had  wandered  to  a  new  location  but  were  clearly 
once  located  at  a  numbered  survey  point  (e.g.,  Dieffenbachia  sp.);  plants  that 
did  not  appear  on  the  maps  and  did  not  have  a  tag,  but  clearly  were 
intentionally  planted  at  one  point  (e.g.,  Eucharis  grandiflora). 
Problems  with  thp  si,irvev  method 

Clonal  or  creeping  understory  species  and  climbing  plants  and  canopy 
vines  were  not  mapped  and  total  number  of  individuals  for  viney  species 
that  root  at  the  nodes  was  not  quantified  during  this  period  due  to  the 
uncertainty  of  delimiting  an  individual.    This  included  Scindapsus  aureus, 
Syngonium  podophyllum,  Passiflora  coriacea,  Ipomoea  sp.,   Cyperus 
alternifolius,  Tradescantia  sp.,  Calissia  fragrans.    Likewise,  the  spread  of  clonal 
ginger  belt  species  was  not  mapped.  These  plants  were  simply  identified  as 
alive,  dead,  or  not  found. 

Not  all  of  the  planter  pockets  on  the  mountain  could  be  accessed,  and 
assumptions  were  made,  in  some  cases  with  only  limited  visibility,  that  since 
there  was  no  longer  a  functioning  irrigation  system  to  the  pockets  that  the 
originally  placed  plants  were  dead.   This  included  fewer  than  20  individual 
plants. 

Species  lists  and  assembly  of  data 

For  the  total  species  count,  individuals  within  a  genus  that  were 
unidentified  to  species  were  counted  as  a  separate  species  only  if  there  were 
no  other  species  within  the  genus.   If  there  were  other  species  in  the  genus, 


the  unidentified  individuals  were  considered  unknown,  since  they  could 
have  overlapped  with  species  already  counted.   Unknown  species  were  not 
included  in  species  abundance  figures  or  growth  form  spectra.  Species  that 
were  classified  as  "not  found"  in  a  survey  were  considered  dead  for  purposes 
of  this  study. 

Surveys  of  Species  Found  per  1000  Individuals  Counted 

Biosphere  2  rainforest 

A  count  of  the  cumulative  number  of  species  found  per  cumulative 
number  of  individuals  was  made  in  the  Biosphere  2  lowland  rainforest 
habitat  in  April,  1998.  Plants  were  identified  and  recorded  in  1.5  m  belt 
transects  covering  the  entire  lowland  rainforest  and  the  east-facing  ginger 
belt.   As  each  individual  was  identified,  it  was  recorded  in  a  column  with  the 
proper  species  label  until  there  were  a  certain  number  of  individuals  per 
column,  starting  with  1  in  the  first  column,  10  in  the  next,  and  about  100  in 
the  subsequent  columns.  All  of  the  plants  in  the  lowland  rainforest  area  were 
counted.   A  diagram  of  the  transects  was  drawn  and  transferred  onto  a  map  of 
the  area  so  that  similar  counts  could  be  made  using  maps  from  previous 
years. 

The  1991  count  of  cumulative  species  per  cumulative  individuals  was 
completed  from  the  original  map  and  species  list  from  1991.   The  1.5  m  wide 
belt  transects  drawn  on  the  maps  were  sampled  in  the  same  order  as  they 
were  sampled  on  the  ground  in  1998.   The  center  of  each  plant  was  included 
only  if  the  surveyed  center  of  the  plant  (as  recorded  on  the  map)  was 
encountered  within  the  transect;  overlaps  were  not  counted.    In  addition, 
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species /individual  counts  were  made  for  the  ginger  belt,  bamboo  belt,  cloud 
bowl,  varzea,  and  mountain  terraces  on  the  1991  map. 
Puerto  Rico  rainforpst 

The  same  field  method  was  used  in  Puerto  Rico  on  a  landslide  near  the 
EI  Verde  Field  Station,  but  the  belt  transect  zig-zagged  across  the  landslide 
area.      All  plants  greater  than  0.5  m  tall  within  a  1  meter  distance  on  either 
side  of  the  observer  were  counted  and  identified.  Where  there  was 
uncertainty  as  to  whether  or  not  2  plants  were  the  same,  collections  were 
made  and  checked  later  with  local  botanical  technicians.   Considerable  care 
was  taken  to  stay  within  the  boundary  of  the  11 -year-old  landslide  area.  In 
addition,  databases  were  sampled  for  number  of  species  per  1000. 

Leaf  Area  Index 

Leaf  area  index  in  the  Biosphere  2  rainforest  lowland  habitat  was 
estimated  from  counts  of  vertical  leaf  overlap  using  a  50-foot  extendable  rod 
expanded  upwards  from  a  ground  point.  The  interception  of  every  leaf  above 
the  ground  point  was  counted.   A  record  of  the  number  of  leaf  interceptions 
per  ground  point  was  made  for  30  points  in  the  lowland  forest  habitat  of  the 
Biosphere  2  rainforest  biome  over  a  ground  area  of  approximately  300  m2. 

Number  of  Seedlings  per  0.54  m2 

Thirty  0.54  m2  circular  plots  were  located  10  paces  apart  in  the  rainforest 
lowland  habitat  in  April  1998.   The  plots  were  defined  with  a  plastic  hoop  and 
were  searched  for  seedlings.  The  counts  were  made  over  a  ground  area  of 
approximately  370  m2. 
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Percent  of  Holes  in  Leaves 

A  visual  estimate  was  made  of  the  percent  of  holes  in  leaf  blades. 
Holes  were  counted  if  they  appeared  to  have  been  caused  by  consumers,  and 
were  not  counted  if  they  appeared  to  be  the  result  of  mechanical  injury  from 
wind  or  human  activities.   The  leaves  that  were  selected  were  within  the 
range  of  view  of  the  ground-based  observer.   Ten  leaves  on  each  of  twenty 
trees  were  counted  using  the  following  percentage  categories:  0,  <1,  1-5,  5-10. 
There  were  no  leaves  that  had  more  than  10%  of  the  blade  consumed. 

Number  of  Green  or  Yellow  Fallen  Leaves  per  0.54  m2 

The  0.54  m2  circular  plots  used  for  seedling  counts  were  also  used  for 
counts  of  the  number  of  green  or  yellow  fallen  leaves  per  0.54  m2.  Each  plot 
was  searched  for  yellow  or  green  leaves.  A  leaf  that  was  100%  green  was 
recorded  as  green.  A  leaf  that  had  any  amount  of  yellowing  was  recorded  as 
yellow.  If  over  50%  of  the  leaf  was  brown,  it  was  considered  brown  and  not 
recorded. 

Calculations 

Diversity  Index 

The  Shannon-Wiener  index  of  diversity,  H'  =  -I  p(  In  p(    where  p,    is 
the  proportion  of  individuals  found  in  the  rth  species,  was  calculated  using 
the  rainforest  survey  data  for  1991,  1993,  and  1996.  Calculations  were  made  to 
follow  the  plants  from  the  first  planting  through  all  three  sampling  dates. 
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Biomass  Estimates 

Three  estimates  of  biomass  have  been  made:   one  for  Nov.  1990,  the 
second  for  July  1991  just  prior  to  the  first  closure,  and  one  in  1993  reflecting 
changes  during  the  first  2  years  of  closure.  The  1990  and  1991  estimates  were 
made  by  personnel  from  Yale  University  School  of  Forestry  and 
Environmental  Studies  (Haberstock  1991).   Estimates  were  made  from  size 
attributes  measured  for  individual  plants  inside  the  Biosphere  such  as  height 
and  diameter  at  breast  height,  using  equations  for  trees  from  the  literature 
(woody  trees  except  iovleucaena  from  Scatena  et  al.  1993)  and  an  equation  for 
the  'Musa'  type  developed  by  Haberstock  based  on  small  destructive  sampling 
at  Biosphere  2. 

For  the  1990  estimate,  every  plant  in  the  Biosphere  2  rainforest  was 
measured.  For  the  1991  estimate,  every  large  woody  tree,  about  1/2  of  the 
'Musa'  category,  and  10%  of  most  other  categories  were  measured. 
Measurements  were  made  on  randomly  selected  subsets  of  plant  categories 
other  than  trees. 

The  July  1993  estimate  (Bierner  1993)  was  extrapolated  from  a  subset  of 
tree  measurements,  using  the  equations  cited  above.   The  subset  included  81 
big  trees,  which  were  compared  with  their  previous  measurements  to 
calculate  a  percent  increase  in  biomass.   That  increase  was  assumed  to  be  the 
same  for  all  growth  forms  in  the  rainforest,  and  was  applied  to  original 
biomass  estimates  to  arrive  at  the  1993  estimate. 


Growth  Form  Spectra 

The  symbols  used  to  represent  the  rainforest  plants  on  the  original 
maps  were  changed  for  this  study.  Nine  growth  forms  are  used  to  describe 
the  plants,  as  follows:  tree  (T),  arborescent  palm  and  palmlike  plant  (P),  shrub 
(S),  giant-leaved  herb  (G),  herb  (H),  graminoid  (R),  woody  graminoid,  such  as 
bamboo  (A),  climber  (C),  and  epiphyte  (E).  Discrimination  between  giant- 
leaved  herbs  (Musa,  Strelitzia)  and  herbs  (Calathea,  Maranta)  was  somewhat 
arbitrary,  but  generally  giant-leaved  herbs  have  the  form  of  herbaceous  trees. 
All  herbaceous  and  woody  vines  and  lianas  were  lumped  into  the  climber 
category.  All  bromeliads  and  orchids  were  considered  epiphytes  if  a 
description  of  the  species  was  not  found  to  the  contrary.  Since  these 
categories  do  not  always  agree  with  those  of  other  studies,  certain  categories 
were  lumped  when  comparisons  were  made  with  other  published  data.   In 
particular,  the  following  growth  forms  were  lumped:  Herb  =  Herb  +  Giant- 
leaved  herb  +  Graminoid;  Tree  =  Tree  +  Arborescent  palm  +  Woody 
graminoid. 

Poisson  Distribution 

A  grid  of  squares  was  printed  onto  transparencies,  with  each  square 
representing  a  ground  area  of  2  m  by  2  m.  The  transparency  was  placed  over 
the  map  of  plants  in  the  rainforest  biome,  and  the  number  of  plants  in  each 
square  was  counted  and  recorded  on  the  transparency.  The  entire  rainforest 
biome  was  counted  with  the  exception  of  squares  containing  concrete  slabs, 
water  bodies,  or  other  features  that  would  have  prevented  planting.  A  total 
of  389  squares  were  counted,  or  about  80%  of  the  rainforest  surface  area. 


The  total  number  of  squares  containing  the  same  number  of  plants  was 
tabulated,  and  the  mean  number  of  plants  per  square  was  calculated.   From 
this,  a  Poisson  distribution  was  generated  for  comparison  with  the  actual 
distribution.   A  chi-square  test  was  used  to  make  the  determination. 

The  ratio  of  the  variance  to  the  mean  was  calculated  to  compare  with 
the  same  ratio  of  the  Poisson  distribution,  which  is  equal  to  one.   A  ratio 
greater  than  one  would  imply  contagious  distribution;  less  than  one, 
regularity  (Whittaker  1975). 

Simulation   MpfhnH 

The  simulation  method  used  in  this  dissertation  is  described  in  detail 
by  Odum  and  Odum  (in  press).  The  first  step  of  the  method  is  to  draw  a 
diagram  of  the  system  being  studied  using  symbols  of  the  energy  systems 
language  (Table  1-2).  The  study  system  is  delineated  inside  the  border  of  a 
window  to  include  material  and  information  storages  in  their  hierarchical 
order,  energy  flows,  feedbacks,  and  energy  drains.    External  energy  sources  are 
shown  flowing  in  from  outside  across  the  system  border.   Next,  rate  equations 
are  written  from  the  diagrams.   The  method  is  explained  with  a  simple 
example. 

The  diagram  and  equations  in  Figure  2-2  represent  the  species  diversity 
on  islands  as  described  by  MacArthur  and  Wilson  (1967)  and  interpreted  and 
drawn  by  Beyers  and  Odum  (1993)  using  energy  systems  symbols.    Here 
diversity  (Q)  is  the  result  of  the  constant  inflow  of  seeds,  spores,  and  other 
propagules  immigrating  from  outside  the  system  and  the  outflow  from  the 
system  due  to  linear  extinction.   The  number  of  species  already  established  on 
an  island  system  as  shown  creates  a  backforce  against  colonization  by  new 


species.  The  backforce  is  represented  diagrammatically  by  the  pathway 
between  P  and  S,  which  lacks  an  directional  flow  arrow;  and  mathematically 
in  the  equations  which  show  the  rate  of  flow  as  dependent  on  the  difference 
between  P  and  S. 

Figure  2-3  shows  a  simple,  one-tank  model  of  biodiversity,  where  Q  is 
the  species  diversity  of  the  system;  J  is  the  constant  flow  of  species  from 
outside  the  system  (seeds,  spores,  other  propagules)  immigrating  into  the 
system;  Kl  is  the  pathway  coefficient  for  the  rate  of  extinction,  which  is 
described  with  a  quadratic  drain  representing  extinctions  in  proportion  to  self 
interactions  such  as  interspecies  competition.     K2  is  the  coefficient  for  linear 
extinction  in  proportion  to  species  present.   The  inflow  of  species  in  this 
model  is  independent  of  the  number  of  species  already  present,  as  shown  by 
the  constant  flow  of  species  P  to  the  system  storage  Q. 

For  calibration,  values  of  storages  and  flows  are  placed  on  the  diagram 
to  help  visualize  consistency.    Then  a  calibration  table  is  made  to  calculate 
the  pathway  coefficients.  Table  2-4  calculates  pathway  coefficients  for  the 
model  shown  in  Figure  2-4  when  the  biodiversity  at  steady  state  is  60  species. 

An  EXCEL  spreadsheet  was  used  to  run  the  simple  biodiversity 
simulation.  Results  of  the  simulation  when  started  with  different  values  of 
diversity  are  shown  in  Figure  2-5. 

Emerffv  Fvalyatinn 

Emergy  evaluations  of  Biosphere  2  construction  and  rainforest 
development  were  made  according  to  the  method  of  Odum  (1996),  as 
summarized  below: 
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Figure  2-2.  Diagram  and  equations  for  species  diversity  on 
an  island  as  described  by  MacArthur  and  Wilson  (1967) 
and  drawn  by  Beyers  and  Odum  (1993). 


dQ/dt  =  J-Kl*Q*Q-K2*Q 


Figure  2-3.  Single  tank  model  of  biodiversity  as  a 
balance  of  steady  inflow  and  linear  and  quadratic 
extinction  showing  storages  and  pathway  coefficients. 


Figure  2-4.  Single  tank  model  of  biodiversity  with 
calibration  values  for  flows  and  storages. 


Table  2-4.  Spreadsheet  used  to  calculate  coefficients  for  single  tank 
biodiversity  model  in  Figure  2-3,  calibrated  at  steady  state. 


Storage  or  flow Steady  state  value    Coefficient  calculation 

Storage 

Number  of  species  Q=  60 

Flow 

Yearly  species  additions       j_  \ 
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Linear  species  K2*Q  =0.1  K2  =  0.001667 
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Figure  2-5.  Simulation  of  number  of  species  Q  for  different 
starting  values  using  the  single  storage  model  in  Figure  2-3. 


First,  a  diagram  of  the  system  was  drawn  showing  the  elements  that 
were  to  be  included  in  the  analysis.  Second,  the  emergy  evaluation  table  was 
constructed,  including  columns  for  the  attribute  of  the  system  being 
measured  and  their  energy  values,  transformity  values,  and  emergy  values. 
Others,  such  as  emvalue,  may  be  added  if  relevant  to  the  analysis  being  made. 
The  first  column  of  each  table  is  a  note  number,  referring  the  reader  to  the 
notes  following  the  table  which  explain  the  origin  of  the  values  and /or  the 
calculations  used  to  derive  the  values.   The  table  was  then  filled  in  with 
available  or  calculated  data. 

The  emergy  analysis  of  the  developing  rainforest  of  Biosphere  2 
included  a  summation  of  everything  that  went  into  creating  the  structure  and 
ecosystems  prior  to  the  first  material  closure.   A  percentage  of  the  total  value 
in  proportion  to  the  relative  size  of  the  rainforest  within  the  Biosphere  was 
calculated,  giving  the  total  cost  for  the  rainforest  only.    The  final  analysis  was 
then  used  to  compare  to  the  tabonuco  rainforest  system  of  Puerto  Rico. 


CHAPTER  3 
RESULTS 


Data  and  analyses  are  presented  for  Biosphere  2  and  the  rainforest  at  El 
Verde  in  Puerto  Rico. 

Characterize  Hon  of  the  Rainforest  in  BiosphprP  7 

Soil 

The  rainforest  soils  were  relatively  homogeneous  in  vertical  profile 
after  placement.    Some  textural  discontinuities  resulted  from  uneven  mixing 
of  the  soil  materials  during  their  placement,  creating  sandier  and  rockier  soils 
in  parts  of  the  lowland  rainforest  where  the  soil  is  a  coarse  sandy  loam  or 
coarse  very  sandy  pebbly  loam  throughout  the  vertical  profile   (Scarborough 
1994).   Various  soils  used  for  potting  mixes  of  individual  plants  became  part 
of  the  rainforest  topsoil  after  planting.   Several  of  the  larger  trees  had  nearly  a 
cubic  meter  of  soil  and  root  volume.  In  addition,  clay  aggregates  in  the  soil 
mixture  tended  to  maintain  their  structure  as  'peds'  from  about  1-12  cm 
diameter  (Scarborough  1994). 

Measurements  made  by  Scott  (1999)  for  December  1993  and  those 
reported  by  Lin  et  al.  (1998)  for  later  dates  are  shown  in  Table  3-1.  Vertical 
development  of  the  rainforest  soils  in  the  lowland  and  ginger  belt  habitats 
was  apparent  in  December  1993,  after  3  years  of  emplacement.  Scott  (1999) 
reports  an  accumulation  of  some  elements  (C,  N,  K,  Ca,  Mg)  in  the  upper 
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stratum  (0-10  cm)  of  the  rainforest  soils  during  this  period.    Similar  trends 
were  reported  by  Lin  et  al.  (1998)  in  Jan.  1996  for  percent  organic  matter  and 
N03-N,  though  in  the  same  study  the  accumulation  was  not  apparent  just 
one  month  earlier.    Over  a  short  time  (Dec.  95-July  96)  the  same  study  reports 
a  clear  decrease  in  concentrations  of  N03-N,  P04-P,  and  Fe  in  the  0-20  cm 
stratum  (the  only  layer  for  which  data  are  available  for  both  dates,  n=5). 

The  pH  values  reported  from  Dec.  93  through  Jan.  96  along  a  vertical 
profile  are  slightly  alkaline,  ranging  from  7.34+0.10  to  7.98+0.10.  Scarborough 
(1994)  noted  a  "tendency  towards  precipitation  of  soil  carbonate  minerals  in 
the  rainforest  soils,  based  upon  several  water  chemistry  signatures",  and  that 
the  water  routinely  used  for  irrigation  during  the  first  2  years  was  enriched 
with  bicarbonate  ions. 

Soil  bulk  density  in  Dec.  1993  of  the  0-10  cm  stratum  was  1.10  g  cm"3  for 
the  lowland  and  1.05  g  cm"3  for  the  gingerbelt  (Scott  1999),  Table  3-2.  Bulk 
density  was  higher  in  all  of  the  other  strata  (to  1.32  in  lowland  and  1.11  in 
ginger  belt),  but  no  further  patterns  were  evident.  Per  cent  volume  of  coarse 
fragments  increased  with  depth. 

Anoxic  areas  were  reported  by  Scarborough  (1994)  after  three  years  of 
soil  development.  He  also  reported  "a  great  [deal]  of  soil  homogenization  of 
the  upper  12-16  inches  due  to  worm  activity,"  and  a  gray  surficial  layer  0-7.5 
cm  deep  forming  in  the  NE  quadrant  of  the  lowland  habitat  that  he  attributed 
to  worm  activity.  Soil  fauna  as  described  by  Scott  (1999)  included  worms  to  a 
depth  of  40  cm,  but  not  in  the  surface  litter  which  he  attributes  to  predation  by 
and  avoidance  of  ants.  Other  macrofauna  counted  by  Scott  were  isopods, 
cockroaches,  and  millipedes. 


Table  3-2.  Soil  bulk  density  and  percentage  coarse  fragments 
in  Biosphere  2  rainforest  from  samples  made  in  November  1993. 
From  Scott  (1999). 


Depth 

Bulk  density  Coarse  fragment 

Location 

cm 

g  cm"3                    vol.,  % 

Biosphere  T 

Lowland  rainforest 

0-10 

1.10                        8.23 

10-20 

1.26                       12.88 

20-40 

1.32                       15.22 

40-80 

1.18                       15.54 

Ginger  belt 

0-10 

0.94                        10.94 

10-20 

1.11                        14.60 

20-40 

1.03                       16.75 

40-80 

1.10                       21.13 

'  Data  from  Scott  (1999),  n=3. 


Plants 

There  were  two  periods  of  plant  introductions  to  the  Biosphere  2 
rainforest.  The  first  was  from  April  1990  to  September  1991;  the  second  from 
October  1993  to  March  1994.  The  number  of  individual  plants  and  species  are 
summarized  in  Table  3-3. 

At  the  start  of  the  1991  closure  of  the  Biosphere,  1890  individual  plants 
were  recorded  in  the  rainforest.  Of  the  original  plants  recorded,  316  species  (or 
monospecific  genera)  were  recognized  in  99  families,  with  315  of  the  original 
1890  plants  not  identified  to  species. 

All  of  the  plants  growing  in  the  rainforest  by  September,  1991  were 
mapped.  Every  plant  was  drawn  on  the  map  with  a  symbol  representing  its 
growth  form  and  with  a  unique  survey  number.   Canopy  widths  were 
mapped  to  scale.  Figure  3-1  shows  a  section  of  the  map.  All  of  the  map  is 
shown  in  20  sections  in  Figure  A-l.     Every  survey  point  with  its  individual 
plant  identification,  growth  form,  and  inventory  status  for  1991,  1993,  and 
1996  is  listed  in  Table  B-l. 

During  the  second  planting  period  at  the  end  of  1993,  339  individuals 
were  added  in  at  least  92  species  (or  monospecific  genera),  50  species  of  which 
were  new  to  the  rainforest. 

In  1996  an  additional  48  individual  plants  in  15  species  were  recorded. 
All  of  these  plants  appeared  to  have  been  self-propagated,  some  from  seed 
(such  asPachira  aquatica  and  Coffea  arabica)  and  others  clonally  (such 
asColocasia). 
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Table  3-3.  Total  number  of  species  and  individual 
plants  seeded  in  the  Biosphere  2  rainforest. 

Approximate  Number  of  Number  of 

dates  of  planting  individuals  species 

4/90-9/91                        1890  316 

11/93-2/94                       339  92/50a 

Total           2229  366 


aOf  the  92  species,  50  were  new  to  the  Biosphere. 


Figure  3-1.  Map  section  #3  showing  location,  growth  form,  and  canopy  size 
of  individual  plants  in  the  Biosphere  2  rainforest  in  September  1991.  The 
twenty  map  sections  are  given  in  Figure  A-l.  Approximate  scale:  1  cm  =  0.6  m. 
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Plant  survey  results,  first  planting 

At  the  end  of  1993,  a  total  of  872  (46%)  of  the  original  plants  (first 
planting)  were  still  alive.    In  1996,  529  (28%)  of  the  original  plants  were 
found.    Table  B-2  shows  the  number  of  individuals  within  each  species  at  the 
time  of  closure  in  1991,  along  with  the  subsequent  inventories  for  1993  and 
1996. 

In  the  1993  inventory,  194  (61%)  of  the  original  species  persisted  and  by 
1996,  137  (43%)  of  the  original  species  remained.   The  distribution  of  plants 
within  species  for  both  inventory  years  shows  a  decline  in  both  numbers  of 
species  as  well  as  number  of  individuals  per  species,  as  summarized  in 
Table  3-4. 

All  or  most  of  the  plants  in  the  following  groups  from  the  first 
planting  died  during  the  first  2  year  study  period:  Adiantaceae  (ferns), 
Aspleniaceae  (ferns),  Blechnaceae  (ferns),  Cyatheaceae  (tree  ferns), 
Orchidaceae,  Polypodiaceae  (ferns),  Selaginellaceae  (fern  allies),  and  most 
Bromeliaceae  (mostly  epiphytes);    Sagittaria,  Pontederia,  Typha,  Nymphaea, 
and  Lycopodium.    This  included  almost  all  of  the  herbaceous  and  tree  ferns, 
fern  allies,  epiphytes  and  aquatic  plants.  The  large  herbaceous  species  planted 
largely  in  the  peripheral  ginger  belt  thrived,  spreading  clonally.   The  most 
successful  taxa  in  terms  of  maintaining  large  populations  were  all  clonal 
species  including  the  following  taxa:   Musaceae,  Marantaceae,  Zingiberaceae, 
and  Strelitziaceae. 


Table  3-4.  Number  of  individual  plants  and  species  recorded  in 
1991, 1993,  and  1996  surveys  of  Biosphere  2  rainforest.  The  first 
number  in  each  entry  includes  plants  from  the  1991  planting,  the 
second  number  from  the  1993  planting,  and  third  is  plants  that 
have  self-propagated.  Species  reported  are  for  species  new 
to  the  rainforest. 


Survey 
date 


Number  of        Number  of 
individuals        species 


1991         1890  316 

1993         872/339  194/50 

1996         529/86/48         137/20/0 


Plant  survey  results   second  planting 

In  1996,  41%  of  the  second  planting  species  and  25%  of  the  individuals 
remained.  Table  A-3  lists  the  plants  that  were  added  to  the  rainforest  after 
September  1993,  and  their  subsequent  inventory  for  1996.    In  contrast  to  the 
first  planting,  many  of  the  Bromeliaceae,  largely  epiphytes,  survived.   Most  of 
the  plants  in  the  following  families  died:  Orchidaceae,  Arecaceae  (palms), 
Aspleniaceae  (ferns),  and  Polypodiaceae  (ferns). 
Abundance  and  distribution  nf  plants  within  speries 

The  change  in  the  number  of  plants  within  species  over  5  years  is 
shown  for  the  plants  from  the  1991  planting  in  Figure  3-2.  For  all  years,  there 
were  relatively  few  common  species  and  many  rare  species.  In  1991,  90 
species  were  represented  by  only  one  individual,  whereas  only  about  15 
species  were  each  represented  by  more  than  15  individuals.   About  45%  of  the 
species  in  1991  had  only  1  or  2  individuals,  which  represented  11%  of  the 
plants.   By  1993,  both  number  of  species  and  numbers  of  individuals  within 
species  had  declined,  and  56%  of  the  species  had  only  1  or  2  individuals, 
representing  only  15%  of  the  individuals,  a  pattern  repeated  again  in  1996. 
Individual  and  sperieS  decline 

The  mortality  rate  of  individuals  and  extinction  of  species  was  higher 
for  the  second  planting  than  the  first  over  the  period  1993-1996  as  shown  in 
Figures  3-3  and  3-4.  Seventy-five  percent  of  the  second  planting  individuals 
died,  whereas  only  39%  of  the  first  planting  individuals  died;  and  about  60% 
of  the  species  from  the  second  planting  and  30%  of  the  species  from  the  first 
planting  went  extinct  over  the  same  period. 
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Figure  3-2.  Change  in  plant  species  abundance  and 
distribution  in  Biosphere  2  rainforest,  1991-1996,  for 
plants  from  the  first  planting. 
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Figure  3-3.  Number  of  individual  plants  in  Biosphere  2 
rainforest  through  time  including  2  planting  periods 
1991  and  1993. 
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Figure  3-4.  Number  of  species  in  Biosphere  2  rainforest 
through  time  including  2  planting  periods,  1991  and  1993. 


The  decline  of  number  of  species  was  proportional  to  the  decline  of 
number  of  individuals,  Figure  3-5. 
Change  in  growth  form  sperrrnm 

The  decline  of  species  and  individuals  within  each  of  9  growth  forms  is 
shown  in  Figure  3-6  and  3-7,  respectively,  for  the  plants  from  the  initial  1991 
planting.   The  largest  percentage  of  survival  in  individual  plants  for  both 
1993  and  1996  was  for  bamboos,  graminoids  and  giant  herbs.  The  decline  of 
absolute  numbers  of  individuals  was  greatest  for  herbs,  trees,  and  climbers, 
which  also  started  out  with  the  largest  number  of  individuals. 

Giant  herbs,  bamboos  and  palms  showed  the  smallest  percentage 
decline  of  species  by  both  1993  and  1996.  Herbs  and  epiphytes  had  the  largest 
percentage  decline  in  species  by  1993,  and  epiphytes,  trees  and  climbers  lost 
the  largest  percentage  of  species  between  1993-1996.  The  epiphyte  growth 
form  was  extinct  by  1996. 

For  comparative  purposes,  plant  growth  forms  in  Figures  3-8  and  3-9 
were  lumped  the  into  5  types  that  are  most  commonly  used  in  other  studies: 
Herb  =  herb  +  giant  herb  +  graminoid;  tree  =  tree  +  arborescent  palm  +  woody 
graminoid.   Though  the  number  of  species  with  each  growth  form  declined 
with  time,  the  relative  proportions  of  species  within  each  growth  form 
changed  only  slightly,  with  the  greatest  change  seen  in  the  increase  in  the 
relative  proportion  of  tree  species. 
Plant  reproduction 

Most  of  the  introduced  pollinators  died  by  the  second  year  of  closure. 
However,  the  following  plant  species  set  fruit  at  least  once  during  the  first  2 
years:    Basella  alba,  Bixa  orellana,  Pachira  aauatica,  Carica  papaya,  Carina 


100  1000  10000 

Number  of  individuals 

Figure  3-5.  Relationship  of  number  of  species  to 
number  of  individuals  over  time  in  Biosphere  2 
rainforest  on  linear  (upper  figure)  and  semilog  (lower 
figure)  scale.  Open  triangles  are  plants  from  1991 
planting,  squares  from  1993  planting,  and  circles 
combine  plants  from  both.  Survey  data  are  from  1991 
1993,  and  1996. 
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Figure  3-6.  Number  of  individuals  within  growth  forms 
in  Biosphere  2  rainforest  in  1991, 1993,  and  1996.  Data 
include  only  plants  from  first  planting.  H=Herb,  T=Tree, 
S=Shrub,  C=Climber/  G=Giant  herb,  P= Arborescent  palm, 
E=Epiphyte,  R=Graminoid,  A=Bamboo. 
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Figure  3-7.  Number  of  species  within  growth  forms  in 
Biosphere  2  rainforest  in  1991, 1993,  and  1996.  Data 
include  plants  from  first  planting.  H=Herb,  T=Tree, 
S=Shrub,  C=Climber,  G=Giant  herb,  P=Arborescent  palm, 
E=Epiphyte,  R=Graminoid,  A=Bamboo. 


Figure  3-8.  Change  in  number  of  species  within  growth  forms 
from  1991  planting  in  Biosphere  2  rainforest  through  time. 
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Figure  3-9.  Change  in  distribution  of  species  within  growth 
forms  from  1991  planting  in  Biosphere  2  rainforest  through  time. 


edulis,  Coffea  arabica,  Ficus  buxifolia,  Moringa  oleifera,  Leucaena  spp., 
Ricinus  communis,  Ipomoea  sp.,  Passiflora  sp.  and  Ceiba  pentandra.    Viability 
of  the  seeds  is  not  known.  Seedlings  of  the  following  species  were  observed: 
Coffea  arabica,  Pachira  aquatica,  Carina  edulis,  and  Leucaena  spp.    In 
subsequent  years,  additional  species  have  set  fruit.   However,  some  plants 
that  reproduced  before  the  1991  closure  failed  to  reproduce  between  1991-1993. 
Examples  are  Clitoria  racemosa,  Averrhoa  carambola,  and  Myrciaria 
cauliflora.    Pachira  aquatica   trees  had  not  reproduced  before  their  planting  in 
the  Biosphere. 

Spatial  Distribution  nf  Plan*: 

The  Biosphere  2  rainforest  plants  were  planted  without  a  formal  plan 
for  their  spatial  distribution.   A  comparison  of  the  distribution  of  the  actual 
plants  with  the  Poisson  distribution  is  shown  in  Figure  3-10.  The  probability 
that  the  actual  distribution  of  plants  was  only  a  change  departure  from  the 
Poisson  distribution  was  less  that  .001.   Thus  the  actual  distribution  is  not 
random.   The  variance  to  mean  ration  of  the  actual  distribution  was  0.355, 
suggesting  regularity  in  the  horizontal  distribution  of  plants  (Whittaker  1975). 

Cumulative  Species  per  Cumulativp  Individuals  Tnnn^H 

Plant  species  in  the  Biosphere  2  rainforest  in  1991  accumulated  as  a 
function  of  the  number  of  individuals  counted,  shown  in  Figure  3-11.   Data 
are  plotted  on  linear,  semi-logarithmic  and  double  logarithmic  scales.    Counts 
made  in  the  lowland  and  east  ginger  belt  habitats  of  the  rainforest  in  1998  are 
graphed  In  Figure  3-12.  The  graphs  show  species /individual  counts  both  with 
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Number  of  individuals  in  2m  x  2m  plots 

Figure  3-10.  Spatial  distribution  of  plants  in  Biosphere  2 
rainforest  in  1991  (solid  squares)  compared  to  Poisson 
distribution  (open  triangles). 
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Figure  3-11.  Cumulative  plant  species  as  a  function  of 
cumulative  individuals  in  Biosphere  2  rainforest  1991 
Data  for  all  habitats  except  cliff  faces  and  surface' 
aquatic  systems  are  plotted  on  linear  (a),  semi-log  (b) 
and  log-log  (c)  scales.  5 
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Figure  3-12.  Cumulative  plant  species  as  a  function  of 
individuals  in  Biosphere  2  rainforest,  April  1998.  Open 
squares  have  Scindapsus  and  Syngonium   subtracted. 


and  without  the  2  most  common  species  that  formed  a  continuous  ground 
cover  in  some  areas,  Scindapsus  aureus  and  Syngonium   podophyllum. 
Together  they  accounted  for  over  half  of  the  individual  plants  in  1998. 
Without  the  2  species,  the  slope  of  the  graph  is  steeper  since  more  species 
accumulate  with  fewer  individuals.  Data  for  the  2  years  are  graphed  together 
in  Figure  3-13. 

Leaf-Area  Index 

The  leaf-area  index  averaged  for  30  observations  in  April,  1998,  was  6.07 
cm2/ cm2,  Table  3-5.  The  range  of  leaf  interceptions  per  point  was  0  -  16. 
Notes  from  field  observations  made  in  November,  1993  (Odum  et  al.)  show  a 
leaf-area  index  of  4.9  averaged  for  10  observations,  with  a  range  of  leaf 
interceptions  from  2-7.   The  mean  number  of  interceptions  of  branches  was 
0.9. 

Number  Of  feedlinPS  and  Green  and  YpTIqw  Fallen  LjagfiS  p»r  m2 

In  1996,  a  single  Leucaena  seedling  was  found  in  one  of  the  30,  0.54  m2 
circular  plots,  Table  3-6. 

Green  and  yellow  fallen  leaves  counted  per  circular  plot  were 
distributed  as  shown  in  Table  3-7.   The  average  number  of  green  leaves  per 
plot  was  1.8  with  a  range  from  0-8  leaves,  and  the  average  number  of  yellow 
leaves  per  plot  was  1.7,  with  a  range  of  0  -  6  leaves. 
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Figure  3-13.  Cumulative  species  as  a  function  of  cumulative 
individuals  in  Biosphere  2  rainforest  in  1991  (solid  circles) 
1998a  (all  plants,  open  squares),  and  1998b  (Syngonium  and 
Sandapsus  deleted,  crosses). 
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Table  3-5.  Leaf  and  branch  interceptions  per  observation 
through  vertical  points  in  Biosphere  2  rainforest.  Data 
for  1993  from  Odum  et  al.  (1993). 

NoTof 
Date      observations  No.  of  interceptions  per  observation 
Leaves 
Nov.  93a  10  2,5,7,4,5,6,6,5,5,4 

Apr.  98  30  7,8,3,4,4,7,3,4,10,7,4,16,1,7,5, 

13,7,4,7,0,7,7,10,5,4,1,3,4,14,6 
Branches 

Nov.  93a  10  0,0,0,1,3,2,1,0,0,2 

a  Data  from  field  notes  by  Odum  et  al.  (1993),  on  file  at 
Biosphere  2. 


Table  3-6.  Distribution  of  seedlings  found  in  30,  0.54  m2 
circular  plots  in  Biosphere  2  rainforest  understory  in  1998. 

No.  of  No.  of 

seedlings  plots 

0  29 

1  1 
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Percent  of  Hnlps  in  T  p^pc 

Table  3-8  shows  the  mean  values  for  percent  holes  in  tree  leaves. 
Visual  estimates  were  made  for  20  trees  throughout  the  rainforest  in  April 
1998  and  for  7  trees  in  Nov.  1993.  The  range  of  leaf  holes  per  leaf  in  1993  was 
0  -  3.3%,  and  in  1998  the  range  was  0.1%  -  1.7%. 

Litterfall  and  r)prnmposition  Ratps 

Litterfall  and  decomposition  rates  for  the  Biosphere  2  rainforest  were 
reported  by  Nelson  (1999)  for  measurements  made  from  Oct.  1992  to  Sept. 
1993.  Litter  was  collected  monthly  from  40,  0.25  m2  baskets  distributed  among 
the  habitats  in  the  rainforest.  Yearly  dry  weight  of  litterfall  averaged  over  all 
habitats  of  the  rainforest  was  723  +  158  g  m2yr  \  The  total  yearly  litterfall  for 
the  rainforest  biome  weighted  by  the  size  of  each  habitat  in  which  it  was 
collected  was  1078  kg  yr1.  Average  monthly  litterfall  values  ranged  from 
about  40  g  m"2  in  Dec.  1992  to  about  95  g  m"2  in  May  1993.  The  high  value  for 
May  was  attributed  to  tree  fall  and  branch  breakage,  but  could  also  include  leaf 
and  stem  drop  from  rising  temperatures  in  the  higher  levels  of  the  rainforest, 
killing  the  new  growth  of  the  taller  trees  and  resulting  in  increased  litterfall. 
Decomposition  rates  were  measured  for  leaves  of  three  species  in  the 
rainforest,  and  the  average  rate  of  disappearance  was  95%  per  year.  Leaf 
samples  from  two  species,  Leucaena  sp.  and  Cecropia  schreberiana, 
disappeared  completely  within  one  year. 


Table  3-7.  Distribution  of  green  and  yellow  leaves  in  0.54  m2 
circular  plots  on  lowland  rainforest  ground,  April  1998. 


No.  of 

plots 

No.  of  leaves 

Green 

Yellow 

0 

9 

12 

1 

7 

3 

2 

5 

7 

3 

5 

3 

4 

2 

1 

5 

1 

2 

6 

0 

2 

7 

0 

0 

8 

1 

0 

Table  3-8.  Percent  of  holes  in  leaves  in  Biosphere  2  rainforest  trees 
Data  for  1993  from  Odum  et  al.  (1993). 


No.  of  Ivs. 

Percent  holes  in  leaves  per  tree 

Date 

observed /tree 

(mean  of  all  leaves  observed  per  tree) 

Apr.  98  a 

10 

1.6, 1.6,  0.1, 1.3,  0.4,  0.1, 1.7, 1.0,  0.9, 
0.4,  0.4,  0.6, 1.2,  0.2, 1.6,  0.1, 1.4 

21 

0.9 

11 

0.3 

7 

1.1 

Nov.  93b 

12 

1.3,  0.2 

10 

0,0.5 

9 

3.3 

8 

2.6 

6 

1.8 

&^,  ^x,  x  ^  ct,_.   1I1C  uuupoint  or  rne 

ranges  was  used  to  calculate  the  mean  value  per  tree  where 
<1  =  0.5, 1-5  =  2.5,  etc. 
'  Data  from  field  notes  by  Odum  et  al.  (1993)  on  file  at  Biosphere  2. 
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Carbon  Uptake  and  Respiration 

Engel  and  Odum  (1999)  calculated  community  respiration  and  average 
daily  C02  absorption  rates  for  the  entire  Biosphere  using  diurnal  curve 
analyses  of  oxygen  and  carbon  dioxide  measured  in  Jan.  -  Feb.  1995.  They 
found  an  average  community  gross  production  rate  of  23  g  02  m"2  d\ 
community  respiration  rate  of  25  g  02  m 2  d1,  and  average  C02  absorption  rate 
of  0.2  g  C02  m"2  h1. 

Net  ecosystem  exchange  of  C02,  total  soil  respiration  and  canopy 
assimilation  were  reported  by  Lin  et  al.  (1999)  for  the  rainforest  biome  for 
June-July  1996.  During  this  period  the  rainforest  was  separated  with  a  secure 
curtain  from  the  rest  of  the  Biosphere.   Estimates  were  made  from 
measurements  of  diurnal  variation  of  C02  in  the  atmosphere  when 
atmospheric  C02  levels  were  kept  at  about  450  ppm.   Results  are  summarized 
in  Table  3-9  in  units  given  in  the  original  publication.   Average  total  canopy 
assimilation  was  18.8  g  C02  m 2  d1  (day  =  13.5  hours),  or,  stated  as  production, 

(18.8 g C02 m"2 d])(32 g 02  /44 g C02)  =  13.7 g 02 m 2 d1. 
Total  respiration  rates  were  25.1  g  C02  m 2  d1,  or  18.2  g  02  m 2  d  \  and  net 
ecosystem  exchange  per  day  was  6.3  g  C02  m 2  d1 ,  or  -4.5  g  02  m"2  d1.  These 
production  and  respiration  rate  estimates  are  somewhat  lower  than  those 
made  by  Odum  and  Engel,  even  though  measurements  were  made  when 
more  light  would  be  entering  the  Biosphere  and  higher  metabolic  rates  would 
be  possible. 
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Cutting  and  Consumers 

An  estimated  22  hours  per  month  was  used  for  pruning  in  the 
Biosphere  2  rainforest,  and  an  approximation  of  10  kg  per  hour  was  cut.  Rate 
of  cutting  for  the  first  2  year  study  period  was  220  kg  per  month.  A  small 
fraction  of  this  was  used  as  fodder  for  the  domestic  animals,  and  the 
remainder  was  air-dried  and  stored  in  the  basement  area. 

Diversity  Index 

Results  of  calculations  of  the  Shannon- Wiener  index  is  given  in 
Table  3-10,   showing  decreasing  plant  diversity  through  time  in  the  Biosphere 
2  rainforest  for  the  plants  from  the  first  planting. 

Biomass  Estimates 

Aboveground  plant  biomass  estimates  for  three  dates  are  listed  on 
Table  3-11.  For  the  1991  estimate,  the  standard  error  was  estimated  to  be  25% 
of  the  biomass.   The  second  estimate  shows  a  161%  increase  over  the 
first,  in  8  months.  From  July  1991  to  July  1993  there  was  a  97%  increase. 
Estimated  increase  in  biomass  from  Nov.  1990  to  July  1993  was  414%. 

Comparisons  Between  Biosphere  2  and  Tahonuco  Forest  in  Puerto  Rim 

Data  were  assembled  from  secondary  and  plantation  forests  within  the 
tabonuco  forest  type  in  Puerto  Rico  for  comparison  with  trends  in  the 
developing  structure,  diversity  and  processes  of  the  Biosphere  2  forest.   The 
secondary  forests  are  located  in  El  Verde  forest  and  Bisley  watersheds,  and  are 
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Table  3-10.  Shannon- Wiener  index  of  diversity  in 
Biosphere  2  rainforest  for  1991, 1993,  and  1996. 

Location Date Index* 

Biosphere  2  1991  539 

1993  4.91 

1996 4.64 

*  Natural  log  (In)  used  in  calculations. 


Table  3-11.  Aboveground  plant  biomass  estimates  in 
Biosphere  2  rainforest  for  1990, 1991,  and  1993. 
From  Haberstock  (1991)  and  Bierner  (1993). 

Date Biomass  (kg)    Increase 

Nov.  1990"  809 

July  1991a  2112  161% 

Julyl993b  ■    4161 97% 

aData  from  Haberstock  (1991). 
bData  from  Bierner  (1993). 
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successional  forests  that  are  recovering  from  land  uses  and  hurricanes, 
developing  without  human  management.    The  plantation  sites,  often  paired 
with  the  secondary  forest  sites  with  regard  to  age  and  soil  type,  are  single- 
species  tree  plantations  in  which  additional  plants  from  nearby  forests  have 
become  established. 

In  addition,  a  count  of  cumulative  species  per  cumulative  individuals 
was  made  for  a  landslide  site  near  the  El  Verde  Field  Station  in  the  tabonuco 
forest  type  in  Puerto  Rico  in  May  1999.   The  vegetation  on  the  site  has  been 
developing  over  11  years  and  is  described  by  Myster  and  Walker  (1997)  as  early 
successional,   continuing  to  increase  in  plant  diversity  and  biomass.   In 
contrast,  the  Biosphere's  rainforest  system  is  9  years  old  with  declining 
diversity  and  increasing  biomass. 

Physical  Environmpnl- 

Table  3-12  shows  values  for  environmental  variables  that  distinguish 
the  Arizona  and  Puerto  Rico  rainforests.   The  Biosphere  2  rainforest  is 
approximately  12°  further  north  than  the  tabonuco  forest  in  Puerto  Rico,  and 
has  a  greater  range  of  daylength  through  the  year.  The  light  reaching  the 
Biosphere  2  rainforest  is  68%  that  of  the  tabonuco  forest  due  to  shading  and 
reflection  by  the  outside  structure.   Vapor  pressure  deficit  for  the  Biosphere 
from  Oct.  '97  -  May  '98  was  6  -  9  millibars.    Vapor  pressure  deficit  for  El  Verde 
was  less  (1.8  -  6  millibars).  There  was  a  greater  range  of  temperatures,  both 
diurnally  and  seasonally  in  Biosphere  2.   During  the  first  2  years  of  material 
closure,  the  C02  ranged  between  500  and  4200  ppm.  Thermal  stratification  of 


Table  3-12.  Environmental  variables  at  El  Verde  and  Biosphere  2  forests. 
From  Odum  (1970),  Romer  (1985),  Ahrain  et  al.  (1998),  and  Cuevas  et  al.  (1991). 
Dashes  mean  no  data  are  available. 


Variable  and  units 

EI  Verde 

Biosphere  2 

Insolation,  kcal  rrf2  day" 

Above  canopy 

Annual  avg.a 

3830" 

2610c 

Near  forest  floor 

Avg. 

134h 

- 

Daylength,  hrs. 

Avg. 

Range 

12-13 

Vapor  pressure  deficit 

millibars 

Avg. 

5.5d 

- 

Range 

1.8  -  6e 

6-9' 

Temperature,  C° 

Above  forest 

April  -  Dec.  avg. 

22.6g 

- 

Jan.  -  Mar.  avg 

20.6g 

- 

Diurnal  range 

21-25g 

23-42f 

Near  forest  floor 

Avg. 

22.9 

27-29f 

Diurnal  range 

20.7  -  25.2  * 

25-35f 

b  Odum  (1970),  p.  1-217. 

c  Romer  (1985),  p.  32,  reports  253  W  m"2  for  Tucson,  AZ. 
approximately  half  of  ambient  light  entered  Biosphere  2  structure 
0.5(253  W  m"2)(2.30  E-4  kcal  secO^OO  sec  hf  x)(24  hrs  )= 
2612  kcal  m"2  d1 

d  Odum  (1970),  H-41. 

e  Odum  (1970),  B-415. 

1  Ahrain  et  al.  (1998),  monthly  mean  values,  Oct.  97-May  98 

8  Odum  (1970),  B-347. 

h  Odum  (1970),  1-218. 

'  Cuevas  et  al.  (1991). 
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the  air  in  the  Biosphere  due  to  low  air  flow  created  temperatures  at  the 
canopy  top  that  were  too  high  for  plant  growth,  limiting  the  height  of  trees. 

Soil 

The  Biosphere  2  soils  were  alkaline  (pH  range  from  7.4  -  7.98);  those  in 
tabonuco  forests  acidic  (pH  range  from  4.2  -  5.7),  Table  3-13.  Values  for  iron 
were  lower  in  Biosphere  2,  but  other  cations  were  much  higher  in  all  profile 
depths.  Total  percent  nitrogen  is  similar  in  the  two  systems  for  the  upper  10 
cm  of  the  soil  profile,  and  higher  in  Biosphere  2  for  in  the  deeper  parts  of  the 
profile.   The  range  of  both  total  nitrogen  and  organic  matter  for  Bisley 
samples  is  greater  than  that  found  in  the  Biosphere.  Soil  bulk  density  values 
for  both  systems  are  shown  in  Table  3-14.   The  highest  value  from  all  sites 
(1.32  g  cm3  at  20-40  cm)  are  in  the  Biosphere  2  lowland  rainforest  habitat;  the 
lowest  values  from  Bisley  watershed  (0.63+0.04  g  cm"3  at  0-10  cm)  and  the 
North  Cut  Center  (0.66  g  cm"3  at  0-15  cm). 

Plants 

Species  and  individual.; 

Thirty-nine  species,  or  12.3%  of  the  flora  planted  in  Biosphere  2  in  1991 
were  species  that  occur  in  El  Verde,  Table  3-15.  Many  of  those  plants  were 
collected  in  Puerto  Rico.  By  1996,  21  species,  or  54%  of  the  original  Puerto 
Rican  species  were  still  alive  compared  to  43%  of  the  total  species  in  the 
rainforest.   About  13%  (241)  of  the  original  individuals  were  El  Verde  species, 
and  34%  of  the  individuals  were  still  alive  in  1996  compared  to  28%  of  the 
total  original  individuals. 
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Table  3-14.  Soil  bulk  density  in  Biosphere  2  rainforest  and 
Puerto  Rico  tabonuco  forest.  From  Scott  (1999),  Edmisten  (1970), 
and  Silver  et  al.  (1994). 


Depth 

Bulk  density 

Location 

cm 

gem'3 

Biosphere  2a 

Lowland  rainforest 

0-10 

1.10 

10-20 

1.26 

20-40 

1.32 

40-80 

1.18 

Ginger  belt 

0-10 

0.94 

10-20 

1.11 

20-40 

1.03 

40-80 

1.10 

Puerto  Rico 

South  Control  Centerb 

0-15 

1.05 

15-30 

1.02 

30-61 

1.00 

Radiation  Centerb 

0-15 

1.18 

15-30 

1.14 

30-61 

1.11 

North  Cut  Center6 

0-15 

0.66 

15-30 

0.92 

30-61 

1.03 

Bisley  Watershed0 

0-10 

0.63+0.04 

10-35 

0.98+0.03 

35-60 

1.13+0.06 

a  Data  from  Scott  (1999).  n=3,  samples  from  Dec.  1993. 
b  Data  from  Edmisten  (1970).  Composite  sample  from 

each  site. 
c  Data  from  Silver  et  al.  (1994).  Mean  value  and 

standard  error  of  87  measurements. 


109 


Table  3-15.  Plant  species  occurring  in  El  Verde  that  were  originally  planted  in 
Biosphere  2  forest,  with  an  inventory  of  individuals  counted  per  species  in 
1991  and  1996. 


Growth  form    Species3 

Trees  Cecropia  schreberiana  Miq. 

Ceiba  pen tandra  (L. )  Gaertn. 

Cordia  alliodora  (Ruiz  &  Pa  v.)  Oken 

Mangifera  indica  L. 

Psidium  guajava  L. 

Spathodea  campanulata  P.  Beauv.. 

Syzygium  jambos   (L.)  Alston 

Tabebuia  heterophylla  (DC.)  Britton 

Arboreal  palms  and  palm-like  plants 

Cnemidaria  horrida  (L.)  C.  PresI 

Cocos  nucifera  L. 

Cyathea  arborea  (L.)Sm. 

Prestoea  montana  (Graham)  G.  Nicholson 

Shrubs  Coffea  arabica  L. 

Hibiscus  rosa-sinensis  L. 
Theobroma  cacao  L. 

Climbers  Dioscorea  alata  L. 

Marcgravia  rectiflora  Triana  &  Planch. 
Philodendron  angustatum  Schott 

Giant  herbs       Alpinia  purpurata  (Vieill.)  ex  K.  Schum. 

Hedychium  coronarium   Koenig 

Heliconia  caribaea  Lam. 
Musasp. 


Number 

of  plants 

1991 

1996 

11 

5 

6 

5 

2 

0 

5 

0 

8 

2 

1 

1 

7 

2 

3 

0 

1 

0 

2 

2 

3 

0 

1 

0 

18 

13 

4 

3 

11 

0 

3 

1 

5 

1 

7 

1 

4 

4 

19 

8 

8 

2 

31 

24 
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Table  3-15  -  continued. 


Growth  form    Species3 


Herbs  Asclepias  curassavica    L. 

Bacopa  monnieri  (L.)  Wettst. 

Blechnum  occidentale  L. 

Colocasia  esculenta  (L.)  Schott 

Eryngium  foetidum   L. 

Etlingera  elatior  (Jack)  R.M.  Sm. 

Ludwigia  octovalvis   Qacq.)  P.H.  Raven 

Nephrolepis  exaltata  (L.)  Schott 

Poly-podium  aureum  L. 

Polytaenium  feei  (W.  Schaffn.  ex  Fee)  Maxon 

Graminoids      Paspalum  plicatulum  Michx. 

Epiphytes         Guzmania  berteroniana  (Schult.  &  Schult.)  Mez 
Guzmania  monostachia  (L.)  Rusby  ex  Mez 
Polypodium  crassifolium  L. 
Rhipsalis  baccifera  (J.S.  Muell.)  Steam 


Number  oi 

:  plants 

1991 

1996 

4 

1 

1 

0 

1 

0 

1 

1 

1 

0 

4 

1 

1 

0 

32 

0 

1 

0 

5 

0 

a  Plant  list  from  Biosphere  2  cross-checked  with  El  Verde  checklist  of 
Charlotte  M.  Taylor  in  Lawrence  (1996). 
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Growth  form  spertra 

Graphs  of  the  spectra  for  the  2  systems  are  shown  in  Figure  3-14. 

The  growth  form  spectra  of  the  2  rainforests  show  a  much  smaller  proportion 

of  herbs  in  the  tabonuco  forest  than  in  Biosphere  2,  and  a  large  percentage  of 

epiphyte  species. 

Cumulative  species  as  function  of  cumulative  individuals 

The  sites,  size  class  of  plants  measured,  dates  of  measurements,  and 

results  of  number  of  species  per  number  of  individuals  counted  are 

summarized  in  Table  3-16. 

Cumulative  species  encountered  per  1000  individuals  counted  are 

shown  in  Figure  3-15  for  the  landslide  site  sampled  at  El  Verde.  The  data  are 

graphed  alongside  data  for  the  Biosphere  2  rainforest  in  1998,  Figure  3-16.  The 
graphs  show  that  species  are  accumulated  at  about  the  same  rate  until  about 
500  individuals  were  counted;  then  the  Biosphere  2  accumulated  species 
more  quickly  than  the  landslide. 

Data  from  the  Bisley  grid  and  the  Radiation  Center  were  graphed 
similarly,  using  databases  provided  by  El  Verde  LTER  staff.  The  3  Puerto  Rico 
sites  are  graphed  together  on  Figure  3-17.  The  Bisley  and  landslide  graphs  are 
nearly  identical,  though  the  Bisley  data  only  account  for  woody  species  greater 
than  or  equal  to  2.5  cm  at  diameter  breast  height.   The  Radiation  Center 
accumulated  the  least  species  for  1000  individuals  counted.   The  Radiation 
Center  data  set  included  all  plants,  including  sprouts  and  seedlings,  and 
according  to  Taylor  et  al.  (1996)  the  1988  counts  were  made  at  a  time  when  a 
lot  of  seeds  were  sprouting. 


B2-'93  B2-'96 

Location 


Figure  3-14.  Growth  form  spectra  of  Biosphere  2  rainforest 
and  El  Verde  tabonuco  forest.  B2  =  Biosphere  2.  Data  for  El 
Verde  from  Smith  (1970). 
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Figure  3-15.  Cumulative  plant  species  as  a  function  of  cumulative 
individuals  >0.5  m  tall  counted  on  an  11-year-old  landslide  at  El 
Verde  site. 
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Figure  3-16.  Cumulative  plant  species  as  a  function 
of  cumulative  individuals  counted  in  Biosphere  2 
(solid  circles)  and  on  an  11-year-old  landslide  at 
El  Verde  (open  squares). 
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Figure  3-17.  Cumulative  plant  species  as  a  function  of 
cumulative  individuals  counted  on  El  Verde  landslide 
(open  circles),  Radiation  Center  (open  triangles),  and 
Bisley  grid  (crosses). 
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In  Figure  3-18  all  sites  are  graphed.  By  1998,  the  Biosphere  2  rainforest 
had  approached  the  values  of  the  landslide  and  Bisley  rainforest. 
Leaf  area  indirps 

Table  3-17  shows  leaf  area  indices  for  a  number  of  sites  in  the  tabonuco 
forest,  and  for  2  dates  in  the  Biosphere  2  rainforest.   Though  measurements 
for  the  Biosphere  are  sparse,  they  indicate  an  increase  in  leaf  area  through 
time,  approaching  or  exceeding  values  for  the  tabonuco  system.   Data  from 
the  tabonuco  forest  also  show  an  increase  after  disturbance  for  the  Radiation 
Center  and  North  Cut  Center. 
Seedling  dpnsify 

Seedling  density  was  higher  in  the  tabonuco  than  the  Biosphere  2 
rainforest,  Table  3-18.   The  Biosphere  data  reflect  the  low  rate  of  reproduction 
by  seed  for  most  of  the  species. 
Percent  of  holes  in  Ip^vp^ 

Table  3-19  shows  the  extent  of  leaf  surface  herbivory,  measured  by  the 
absence  of  leaf  tissue.   The  percent  of  holes  in  leaves  at  El  Verde  ranged  from 
over  30%  to  less  than  5%,  depending  on  the  species  of  tree  (Mercado  1970). 
The  average  number  of  leaf  holes  was  7%  (Odum  1970b),  compared  to  an 
average  in  Biosphere  2  of  less  than  1  in  1998  and  less  than  2  in  1993. 
Litterfall  and  dpromposition  rate* 

The  mean  litterfall  rate  for  the  Biosphere  2  rainforest  was  slightly 
lower  than  the  mean  range  of  the  tabonuco  forest,  shown  on  Table  3-20. 
However,  the  Biosphere  lowland  rainforest  data  include  large  woody 
branches  from  a  treefall,  which  skews  this  value  upward  from  the  mean 
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Figure  3-18.  Cumulative  plant  species  as  a  function  of  cumulative 
individuals  counted  on  El  Verde  landslide,  Radiation  Center,  Bisley  grid 
and  Biosphere  2  rainforest  for  1991  and  1998.  a)  linear  scale;  b)  semi-loe ' 
scale;  c)  log-log  scale. 
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Table  3-17.  Leaf-area  indices  of  Puerto  Rico  tabonuco  forest  and 
Biosphere  2  rainforest.  From  Odum  (1970),  Jordan  (1969), 
Odum  and  Pigeon  (1970),  and  Odum  et  al.  (1993). 


Leaf  area/ 

Location 

Method 

ground  area 

Puerto  Rico  tabonuco  forest 

EI  Verde3 

Direct  harvest 

6.4 

Ravine 

Subset,  n=2 

2.1 

Slope 

Subset,  n=4 

6.4 

Ridge 

Subset,  n=4 

8.7 

El  Verdeb 

Plumb  line 

6.6 

Slope 

Subset,  n=16 

6.68 

Slope 

Subset,  n=16 

5.60 

Ridge 

Subset,  n=16 

8.60 

El  Verdec 

Plumb  line 

5.24 

Radiation  Center 

1  yr.  after  irradiation 

3.16 

3  yrs.  after  irradiation 

5.51 

North  Cut  Center 

1  yr.  after  cutting 

4.03 

3  yrs.  After  cutting 

4.45 

Biosphere  ?  rainfnrpst 

Entire  rainforestd 

Extended  pole,  1993 

4.9 

Lowland  rainforest 

Extended  pole,  1998 

6.1 

a  Odum  (1970). 

b  Jordan  (1969). 

c  Odum  and  Pigeon  (1970).  1-175. 

d  Odum  et  al.  (1993). 
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Table  3-18.  Seedling  density  in  Biosphere  2 
and  tabonuco  rainforest.  Tabonuco  forest 
data  from  Odum  (1970). 


Location 

Seedlings  per  m2 

Biosphere  2,  Apr.  98 

Puerto  Ricoa 
Control  Center 
Cut  Center 
Radiation  Ctr,  Oct.  64 
Radiation  Ctr,  Dec.  64 

<0.1 

5.12 
11.2 
2.25 
3.95 

1  Odum  (1970),  D-139,  for  20-1  m  2  plots 
in  each  center. 


Table  3-19.  Percent  of  holes  in  leaves  in  Biosphere  2  rainforest 

and  tabonuco  forest  trees.  Tabonuco  forest  data  from  Odum  (1970b). 


No.  of  Ivs. 

Percent  holes  in  leaves  per  tree 

Date 

observed /tree   (mean  of  all  leaves  observed  per  tree) 

Biosphere  2 

Apr.  98  a 

10 

1.6, 1.6,  0.1, 1.3,  0.4,  0.1, 1.7, 1.0,  0.9, 
0.4,  0.4,  0.6, 1.2,  0.2, 1.6,  0.1, 1.4 

21 

0.9 

11 

0.3 

7 

1.1 

Nov.  93b 

12 

1.3,  0.2 

10 

0,0.5 

9 

3.3 

8 

2.6 

6 

1.8 

El  Verde 

(No.  of  baskets) 

1964c 

49 

6.34  ±0.42 

101 

5.8+0.30 

97 

5.61  ±0.28 

243 

5.85  ±0.22 

o — '     — '  •*■   '-'/  >->■«—     -nn-  iiuuuuiiu  (jl  Lilt; 

ranges  was  used  to  calculate  the  mean  value  per  tree  where 
b  <1  =  0.5, 1-5  =  2.5,  etc. 
Data  from  field  notes  by  Odum  et  al.  (1993)  on  file  at  Biosphere  2 
Data  from  Odum  (1970b). 


yearly  value.  The  average  disappearance  rate  of  leaf  material  for  the 
Biosphere  was  higher  than  that  measured  for  tabonuco  forest,  Table  3-21. 
Nutrient  content  of  leaves 

Nutrient  content  of  leaves  for  the  two  forests  are  given  in  Table  3-22.   Values 
overlap  for  all  elements  except  phosphorus,  which  was  an  order  of 
magnitude  higher  in  the  Biosphere  rainforest  than  the  tabonuco.    Both 
potassium  and  nitrogen  in  the  leaves  of  the  Biosphere  were  closer  to  the 
maximum  value  of  the  ranges  given  for  tabonuco. 
Carbon  exchange 

Soil  respiration  and  canopy  assimilation  were  both  lower  in  the 
Biosphere  2  rainforest  tabonuco  forest,  shown  in  Table  3-23.   Carbon 
assimilation  in  the  Biosphere  (measured  in  summer)  was  56  to  76  percent 
that  of  the  tabonuco  forest,  and  soil  respiration  was  40%  that  of  the  tabonuco 
forest. 
Aboveground  biomass  density 

The  aboveground  biomass  density  values  for  different  sites  in  the 
tabonuco  rainforest  vary  between  8  and  24.6  kg  m"2,  Table  3-24.  The  Biosphere 
2  rainforest  was  4.2  kg  m2  in  July,  1993,  and  beyond  that  value  later. 

Overview  Simulation  Models 

Two  different  overview  simulation  models  were  used  to  consider 
processes  affecting  diversity  and  production  in  these  ecological  systems.  The 
first  model  simulates  production,  consumption  and  biomass,  the  resource 
basis  for  diversity.   The  second  model  relates  plant  diversity  to  resources, 
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Table  3-20.  Litterfall  rates  in  Biosphere  2  rainforest  habitats  and 
tabonuco  forests.  From  Nelson  (1999),  Wiegert  (1970), 
Cuevas  et  al.  (1991),  Lugo  (1992),  and  Lodge  (1991). 

Location  Rate,  gm'z  yr"1 


Biosphere  2  a 

Bamboo  belt 

578  +  81 

Terraces 

466  ±  70 

Lowland 

1317  ±283 

W.  Ginger  belt 

492  ±  46 

Varzea 

641  +89 

All  habitats,  mean 

723  ±158 

Puerto  Rico  b 

Radiation  Center 

623 

So.  Control  Ctr. 

571 

Secondary  forest c 

820  ±  100 

Mature  secondary  forest d 

970 

Bisley  and  El  Verde,  mean6 

856 

a  Nelson  (1999).  Data  from  collections  made  monthly 
Oct.  '92  -  Sept.  93,  mean  and  standard  error. 

b  Data  from  Wiegert  (1970),  H-92,  for  litter  collections 
made  2/64-12/64  extrapolated  for  one  year. 

c  Cuevas  etal.  (1991). 

d  Lugo  (1992). 

e  Lodge  (1991). 


Table  3-21 .  Decomposition  rates  of  litter  in  Biosphere  2  rainforest  and 
El  Verde  forest  in  Puerto  Rico.  From  Nelson  (1999)  and  Wiegert  (1970). 


Rateofdisappearanrp 

Location 

Material 

gg'V"1 

Annual,  g  m  "2 

Biosphere  2  a 

Leucaenasp. 

1.00 

Cecropia  scftreberiana 

1.00 

Musa 

0.90 

Average 

1024 

Puerto  Rico  b 

Radiation  Center 

0.850 

705 

So.  Control  Ctr. 

0.519 

265 

'  Data  from  Nelson  (1999). 
b  Data  from  Wiegert  (1970)  p.  H-95. 
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Table  3-23.  Carbon  exchange  in  Biosphere  2  (summer)  and  tabonuco 
rainforests.  Rg=  Soil  respiration,  Ac=  Canopy  assimilation. 
From  Lin  et  al.  (1999)  and  Odum  (1970b). 


Location 

Date  &  time 

Total  Rs 
gC  m^hr"1 

Total  Ac 
gC  m^hr'1 

Biosphere  2 
rainforest3 

EI  Verdeb 

Giant  cylinder 
Ten  prisms 

June-July  1996 
Day  (13.5  hours) 
Night  (10.5  hours) 

1963 

-0.28 
-0.267 

-0.68 

0.38 

0.68 
0.53 

a  Lin  et  al.  (1999). 
b  Odum  (1970b). 

Table  3-24.  Estimates  of  aboveground  biomass  density  in 
Biosphere  2  and  tabonuco  forests.  From  Haberstock  (1991), 
Bierner  (1993),  Ovington  and  Odum  (1970),  Scatena  et  al  (1993) 
and  Lugo  (1992). 


Location                                Biomass  density,  kg  m  "2 

Biosphere  2 

1991 a 

1.0 

1993  b 

4.2 

Puerto  Rico 

Primary  forests: 

Radiation  center0 

24.6 

South  Control  Center  c 

14.5 

North  Cut c 

19.8 

Bisleyd 

22.6 

Secondary  forests: 

Mature  e 

8.0 

Mahogany  plantation  e 

12.5 

a  Haberstock  (1991). 

b  Bierner  (1993). 

c  Ovington  and  Odum  (1970). 

d  Scatena  et  al.  (1993). 

e  Lugo  (1992). 
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mechanisms  of  interaction  and  management,  where  diversity  is  measured  as 
the  total  number  of  species  (species  richness). 

Production  and  Biomass  Minimodel  Description 

Figure  3-19  is  a  model  of  relationships  believed  important  to 
metabolism,  a  systems  hypothesis.  Included  are  the  parts  and  processes  of  the 
Biosphere  2  rainforest  that  affected  production  and  biomass.   The  outer 
window  of  the  systems  diagram  shows  the  material  isolation  of  Biosphere  2 
from  the  earth  as  it  was  operated  during  1991-1993.  The  inner  window 
defines  the  rainforest  biome  as  a  separate  system  of  the  Biosphere  2 
mesocosm,  with  air,  humans,  consumers  and  biomass  moving  across  the 
biome  boundaries.  Labels  on  the  flow  pathways  (Kl,  etc.)  are  the  coefficients 
used  to  calculate  the  rates  of  flow  between  storages. 

Main  pathways  in  the  model  are  described  as  follows:  Approximately 
half  of  the  energy  from  the  sun  reaching  Biosphere  2  is  reflected  from  the 
external  glass  or  absorbed  by  the  stainless  steel  spaceframe  structure, 
represented  by  an  arrow  from  the  sunlight  pathway  deflected  from  the  outer 
window.    Most  of  the  sunlight  entering  the  glass  is  used  by  the  production 
process,  which  has  been  split  between  two  producer  groups  of  different 
physiological  types,  weedy  plants  (generally  eutrophic  vines  and  herbaceous 
perennials)  and  plants  of  a  more  mature  forest.   The  flow  pathway  to  the 
weedy  plants  is  labeled  Kl  and  the  flow  K2  to  mature  plants.  The  sunlight 
which  is  not  used  by  the  production  process  (R)  is  re-radiated  out  of  the 
rainforest. 
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The  2  producer  units  (bullet-shaped  symbols)  receive  flows  of  sunlight, 
atmospheric  C02  (K24  and  K25),  and  feedback  from  the  living  biomass  (K3 
and  K5).  The  3  flows  interact  with  producers  in  the  process  of  photosynthesis. 
The  pathway  K4  is  the  gross  productivity  of  weedy  biomass;  K6  the  gross 
productivity  of  mature  biomass. 

The  respiratory  functions  from  plant  biomass  (K8  and  K16)  show 
biomass  conversion  to  carbon  dioxide  by  microconsumers  (boxes  along  the 
pathway),  which  flows  (K9  and  K17)  into  the  atmospheric  C02  storage,  N. 
Pathways  K12  and  K14  leading  from  the  plant  biomass  storage  tanks  represent 
natural  litterfall  feeding  into  the  soil  organic  matter  storage,  and  K10  and  K7 
from  the  weedy  biomass  storage  are  for  the  litterfall  created  by  humans 
through  pruning,  passing  into  either  the  soil  organic  matter  storage  or  stored 
as  dried  biomass  in  the  basement  of  the  rainforest  biome.    Though  dried,  the 
stored  organic  matter  still  had  a  slow  respiration  rate,  hence  the  pathway  K26 
from  stored  biomass  (B)  to  the  atmospheric  C02  storage  (N). 

The  soil  organic  matter  is  consumed  by  detritivores  such  as 
earthworms  and  arthropods  on  one  pathway  (K27)  where  a  storage  of 
macroconsumer  biomass  (C)  accumulates.   The  feedback  from 
macroconsumer  storage  to  the  production  process  (K31)  works  to  increase 
macroconsumer  production.  Some  of  the  organic  biomass  consumed  is 
returned  to  the  atmosphere  along  the  respiration  pathway  (K30  to  K28).   On 
another  pathway  (K22)  soil  organic  matter  is  decomposed  by  microconsumers 
(represented  by  a  rectangle  along  the  pathway)  and  returned  to  atmospheric 
C02  along  another  respiration  pathway  (K23). 
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The  calibration  values  shown  in  Figure  3-20  were  estimated  by 
extrapolating  storages  for  the  steady  state  condition  that  could  be  expected  in 
the  Biosphere  2  rainforest,  projected  to  occur  after  20  years  of  similar 
operation  (year  2011).  Many  of  the  pathway  coefficients  were  calculated  using 
rates  measured  inside  the  Biosphere  2  rainforest.   An  EXCEL  spreadsheet, 
Table  C-l,  shows  the  initial  calibration  values  for  storages  and  rates  and 
details  the  values  used  in  the  calibration. 

A  minimodel  (B2METAB)  was  written  using  QUICKBASIC,  Appendix 
C-2,  for  simulation  of  the  production  and  biomass  model  described  above. 
Baseline  simulation  runs  are  presented  in  the  text,  and  other  runs  where 
conditions  are  varied  are  in  Appendix  C. 

Results  of  Simulating  the  Production  and  Biomass  Minimal 

Baseline  simulations  of  the  minimodel  for  10  and  100  years  were  set 
with  initial  storage  values  that  were  close  to  the  Biosphere  2  rainforest  values 
at  the  start  of  the  first  closure  in  September  1991.    The  10-year  simulation 
(Figure  3-21)  shows  atmospheric  C02,  soil  organic  matter  and  consumer 
storages  declining  at  first  with  a  rapid  increase  in  plant  biomass.    Initially,  the 
net  ecosystem  production  (Pr  =  Pg/Rs)  was  greater  than  1,  which  reflects  the 
initial  fast  growth  of  primary  producers.  Over  100  years,  shown  in 
Figure  3-22,  the  soil  organic  matter  increases,  and  as  it  reaches  its  maximum 
storage  value  the  system  is  able  to  support  more  macroconsumers,  which 
increase  in  a  pulse.  This  pattern  is  repeated,  showing  the  switch  to  a 
heterotrophic  system  (Pr<l)  when  consumers  pulse,  C02  increases,  and  the 
soil  organic  matter  decreases.   The  amplitude  of  the  pulses  decreases  over 
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Figure  3-21.  Ten-year  simulation  of  production  and 
biomass  minimodel,  B2METAB  in  Figure  3-19.  Final  values 
for  plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass. 
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Figure  3-22.  One-hundred-year  baseline  simulation  of  production 
and  biomass  in  the  minimodel  B2METAB  in  Figure  3-19.  Final 
values  for  plant  biomassare  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass. 
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time  as  macroconsumers  assimilate  plant  biomass  and  build  a  storage.   The 
system  sustains  a  stable  population  of  macroconsumers  over  time  (centuries) 
where  the  mean  value  of  Pr  is  1.   All  storages  move  towards  a  stable  value 
with  the  exception  of  stored  biomass,  which  continues  to  increase  due  to  a 
continuous  harvest. 

The  following  alterations  to  the  baseline  simulation,  such  as  increased 
light  or  reduced  pruning,  demonstrate  the  effects  of  changes  on  the 
developing  ecosystem.  Figures  3-23  and  3-24  compare  the  effects  of  each 
alteration  of  the  system  on  biomass  or  production. 
Increased  sunlipht  (I  =  3170) 

When  sunlight  entering  the  Biosphere  was  increased  by  20%  biomass 
storages  (B,  M,  W,  C,  S)  built  more  quickly  than  at  the  lower  light  levels 
(Figure  3-23),   amplitudes  of  the  consumer  pulses  were  higher,  and  stable 
values  of  net  production,  soil  organic  matter,  and  macroconsumer  storage 
were  approached  more  quickly  than  with  less  light  (Figure  3-24).    The  mean 
value  of  carbon  dioxide  in  the  air  was  lower,  and  the  sum  of  plant  biomass 
storages  after  100  years  was  greater  than  with  lower  light. 
No  weedv  biomass  (W=n) 

When  the  minimodel  is  started  and  run  without  weedy  biomass,  the 
mature  plant  species  biomass  increases  faster  and  has  a  higher  steady-state 
value  than  the  baseline  simulation  value.    Total  plant  biomass  accumulation 
is  less  since  there  are  neither  weedy  plants  nor  stored  weedy  biomass. 
Consumer  and  soil  organic  matter  storages  grow  more  slowly,  and  consumer 
pulses  are  smaller.   Though  the  effect  of  no  weedy  species  is  to  increase 
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Figure  3-23.  Comparison  of  simulations  of  biomass  by  the  model 
in  Figure  3-19  showing  effects  of  light,  elimination  of  weedy 
biomass,  limited  pruning,  and  altered  airflow  on  mature  and  weedy 
plant  biomass.  Values  for  biomass  at  100  years  are  given  on 
the  right  side  of  the  graphs  in  grams  of  dry  biomass  per  square  meter. 
Scales  are  the  same  for  each  box.  M=Mature  plants,  W=Weedy  plants. 
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Figure  3-24.  Comparison  of  effects  of  light,  elimination  of  weedy 
biomass  limited  pruning,  and  altered  airflow  on  production-respiration 
ratio,  Pg/Rp,  from  simulations  of  the  model  B2METAB  in  Figure  3-19. 
Scales  are  the  same  for  each  box,  with  the  center  line  =  1. 


mature  species  biomass,  total  ecosystem  development  is  slowed  and  total 

storages  are  less  after  100  years. 

All  pruned  biomass  added  to  soil  (K18*H*W  added  to  PS.  B  removed) 

With  diversion  of  the  pruned  biomass  from  dry  storage  to  the  soil 
(additional  lirterfall),  consumer  storages  build  more  rapidly  and  are  larger. 
Biomass  storages  and  COz  are  virtually  unchanged. 
No  removal  of  wepdv  biomass  fH  =  n) 

Pruning  of  weedy  biomass  during  the  start-up  of  Biosphere  2  was 
thought  to  be  important  for  the  survival  of  some  of  the  key  mature  species  by 
releasing  them  from  light  competition  with  the  weedy  species.   The 
simulation  of  the  system  with  no  pruning  shows  the  increased  growth  of 
weedy  species  and  eventual  total  decline  of  mature  species.  A  larger  plant 
biomass  storage  occurs  with  only  weedy  species  than  with  only  mature  species 
(Figure  3-23),  since  the  weedy  plant  type  has  the  higher  ratio  of  Pg/Rp. 
However,  the  runs  with  both  weedy  and  mature  plants  with  weedy  plant 
pruning  show  the  largest  combined  plant  biomass  storage  over  time.   The 
rapid  increase  of  plant  biomass  when  no  pruning  occurs  promotes  much 
more  rapid  storage  and  cycling  of  soil  organic  matter  and  consumer  pulsing 
than  the  baseline  run. 

Whether  weedy  plant  biomass  could  reach  this  amount  with  climbers 
or  herbaceous  perennial  growth  forms  depends  on  the  ability  of  vines  to  use 
the  superstructure  of  Biosphere  2  to  overgrow  other  species.  To  achieve  the 
biomass  storage  without  well-developed  vertical  structure  is  unlikely.  Inside 
Biosphere  2  the  spaceframe  structure  was  used  by  climbers  for  support,  and  to 
increase  the  volume  in  which  they  were  able  to  produce  biomass. 
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Reduction  of  human  effort  in  removing  weedv  biomass  (H  =  .75) 

A  reduction  of  human  effort  in  removing  weedy  species  was  simulated 
to  see  whether  the  amount  of  pruning  would  change  total  plant  biomass 
production  -  the  sum  of  live  plant  biomass  and  dried,  stored  biomass  -  and 
lower  the  atmospheric  C02  concentration.   Figure  3-23  shows  that  combined 
plant  biomass  storages  would  be  increased  over  the  baseline  values  and  C02 
reduced  when  H  is  reduced.  Weedy  biomass  increases  and  stored  biomass 
increases,  though  mature  biomass  decreases,  with  lower  COz  compared  with 
the  baseline  simulations  where  H=l.   However,  below  a  certain  level  of 
human  intervention,  mature  species  will  be  out-competed  by  the  weeds. 
No  airflow  between  hiomes  (Ta  =  0  and  K21  =  ftt 

The  atmosphere  of  the  Biosphere  2  rainforest  connected  to  the 
atmosphere  of  the  other  biomes  within  the  structure.    Without  airflow 
between  the  different  areas,  atmospheric  C02  level  increased,  plant  growth 
increased,  soil  organic  matter  increased,  and  consumers  increased.   The  Pr 
ratio  of  the  system  was  greater  than  one  until  the  pulse  of  consumers 
occurred,  when  it  switched  to  a  heterotrophic  system  (Figure  3-24).  After 
consumers  used  up  the  organic  matter  resource,  they  crashed  and  the  Pr  ratio 
of  the  system  returned  to  the  pre-pulse  condition. 

Production  and  Diversity  Minimodel  Description 

That  diversity  and  structure  of  an  ecosystem  compete  for  energy  is 
especially  evident  during  the  early  stages  of  ecosystem  development.  In  a 
system  with  limited  resources,  amount  and  quality  of  energy  may  determine 
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the  extent  to  which  diversity,  structure  and  their  connections  can  be 
supported. 

A  model  SPDIV  was  developed  to  represent  the  main  factors 
supporting  biodiversity  in  the  rainforest  study  areas  (Figure  3-25).   The 
number  of  species  at  any  time  is  the  balance  between  species  introduced  and 
those  becoming  extinct.   Species  are  in  two  categories,  those  just  planted  (N) 
but  not  yet  interacting  with  the  surrounding  ecosystem  and  those  which  have 
become  established  (E).  The  latter  species  cause  additional  productivity  in 
proportion  to  their  number  by  drawing  energy  from  the  main  plant  biomass 
in  proportion  to  the  interactions  of  species  (E2). 

The  main  pathways  in  the  model  are  as  follows:   Sunlight  and  species 
are  the  two  outside  sources  that  pass  from  the  outside  to  inside  the  boundary 
of  the  Biosphere  structure.   About  half  of  the  ambient  sunlight  (I)  enters  the 
rainforest,  most  of  it  used  for  the  production  process  with  a  small  re-radiated 
portion  labeled  R.  The  production  process  is  shown  with  2  pathways.  The 
flow  of  sunlight  labeled  Kl  interacts  with  the  diversity  of  established  species 
(an  interaction  of  flow  Kl  and  force  from  E),  and  the  flow  labeled  K2  interacts 
with  a  flow  which  is  fed  back  from  biomass,  K4,  in  the  process  of 
photosynthesis.   The  result  of  those  interactions  is  the  gross  productivity  of 
the  system,  pathways  K3  and  K5.  The  2  gross  productivity  pathways  flow  into 
a  common  biomass  tank,  with  the  respiratory  costs  and  consumer  losses 
associated  with  maintaining  the  system  shown  as  a  drain  from  biomass 
labeled  Kll. 

Major  plantings  of  new  species  to  the  Biosphere  2  rainforest  occurred 
twice,  one  with  the  initial  planting  of  about  300  species  during  1990-1991,  and 
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another  with  about  92  species  (50  of  them  were  new  to  the  system)  after  the 
first  closure  period,  1993-1994.  The  pathway  of  species  introductions  is  labeled 
Ja.  The  value  of  Ja  was  0  in  the  first  two  years  after  the  initial  planting.  After 
planting,   species  were  considered  transient  since  they  lacked  connections  and 
interactions  which  normally  develop  with  surrounding  organisms  and 
environment  as  an  ecosystem  matures.   In  proportion  to  organic  habitat  (B) 
some  of  the  planted  species  developed  those  connections  and  interactions  and 
became  established.  The  transfer  of  species  is  shown  as  pathway  K8  which 
interacts  with  a  force  from  biomass  storage  and  enters  the  established  species 
tank  as  flow  K9.  The  number  of  species  becoming  established  depends  on  the 
biomass  (representing  suitable  habitat  and  soils).   Some  species  became 
extinct,  as  shown  with  the  pathway  K10  draining  the  transient  species  storage. 
When  a  transient  species  (N)  becomes  established  (E),  it  has  additional 
opportunities  for  extinction  and  maintenance.   As  part  of  the  species 
assemblage,  it  consumes  energy  through  competitive  interactions  with  other 
species,   through  the  work  of  strategies  that  prevent  competition,  or  through 
other  interactions  with  the  environment.    Some  established  species  become 
extinct,  shown  as  linear  pathway  K13  and  quadratic  K14  that  represents  effect 
of  competition  from  the  established  species  tank.   The  energy  used  to 
maintain  the  species  information  flows  from  biomass,  K12,  which  interacts 
with  established  species.   Given  the  time  and  opportunity  for  evolution,  K6  is 
the  pathway  along  which  new  species  would  be  added  (K6  >  K7). 

In  microcosms  that  have  been  sealed  for  long  periods  without  species 
being  added,  life  appears  to  persist.  Examples  were  not  found  where  last 
species  or  several  species  went  extinct.  Apparently  some  species  are  capable  of 
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operating  an  ecosystem  with  low  diversity.   The  model  has  a  provision  so 
that  extinctions  do  not  continue  once  the  species  number  is  less  than  the 
threshold  C.  This  is  diagrammed  as  a  switch  X.  The  switch  is  set  to  1  when 
the  number  of  established  species  is  greater  than  C,  and  to  zero  when  the 
number  of  established  species  is  less  than  C,  preventing  further  extinctions 
while  still  using  energy  from  biomass  storage  for  maintenance  of  extant 
species. 

The  calibration  values  shown  in  Figure  3-26  were  made  from 
estimations  of  biomass  from  measurements  made  in  Biosphere  2, 
productivity  rates  measured  for  the  El  Verde  rainforest  (adjusted  at  one-third 
of  that  value  for  the  system  inside  Biosphere  2),  and  extinction  rates 
experienced  in  the  Biosphere  2  rainforest  after  the  first  5  years  of  operation. 
Details  of  the  calculations  are  presented  as  notes  to  the  EXCEL  spreadsheet 
used  to  make  the  calculations,  Table  D-l. 

A  minimodel  was  written  using  QUICKBASIC,  Table  D-2,  for 
simulation  of  the  production  and  diversity  model  described  above. 

Results  of  Simulating  the  Production  and  Divprsitv  MinimoHpl 

A  baseline  simulation  of  the  minimodel  was  run  for  both  10  and  100 
years,  starting  with  values  close  the  1991  values  of  Biosphere  2  (biomass  = 
1000  grams  dry  weight  per  square  meter,  transient  species  =  300,  established 
species  =  60  ).  Figure  3-27  shows  that  within  the  first  5  years,  the  transient 
species  have  all  either  become  extinct  or  become  established  species  using 
energy  from  the  system  to  support  the  additional  information.   As  the  total 
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Figure  3-27.  Baseline  simulation  of  diversity  and  production 
minimodel,  SPDIV,  in  Figure  3-25  for  Biosphere  2  rainforest  for 
10  years.  Final  values  for  plant  biomass  and  species  are  on  the 
right  side  of  their  graphs.    g/m2  =  grams  of  dry  biomass  per 
square  meter. 
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number  of  species  in  the  system  (transient  plus  established)  declines  during 
the  first  year,  the  biomass  increases.  As  the  number  of  species  in  the  system 
reaches  a  relatively  stable  level  after  about  3  years,  the  biomass  also 
approaches  a  stable  value.  When  the  established  species  begin  to  decline 
around  year  5  due  to  interactions  with  other  species  and  factors  other  than  the 
amount  of  energy  available,  biomass  increases,  using  the  energy  once  used  for 
species  support  for  biomass  production.  The  100  year  simulation, 
Figure  3-28,  shows  the  continual  increase  of  biomass  as  diversity  declines. 
When  there  are  only  10  species  remaining,  the  switch  in  the  model  prevents 
further  decline  and  biomass  also  remains  stable,  dividing  energy  between 
information  support  and  biomass  production. 
Continuous  addition  of  species  (5,  or  20  per  year) 

The  simulation  results  show  the  effects  of  continuous  additions  of 
species  to  the  system  over  10  and  100  years.  The  10  year  simulation  shows  5 
and  20  species  per  year  added  to  the  system  (Ja=5/365  and  Ja=20/365),  resulting 
in  higher  diversity  after  10  years  of  both  transient  and  established  species  and 
lower  biomass  as  the  diversity  increases.   The  100  year  simulation  with 
additions  of  20  species  per  year  shows  higher  final  diversity  than  the  baseline 
-  25  transient  species  rather  than  none,  150  established  species  rather  than  10; 
and  lower  biomass,  4290  g/m2  rather  than  26000. 
One-time  addition  of  many  sppHpc; 

A  simulation  was  run  with  the  addition  of  50  new  species  after  the  first 
2  years  of  operating  Biosphere  2.  After  10  years,  the  effect  is  to  increase  the 
number  of  established  species  to  112  above  the  baseline  value  of  96,  with 
accompanying  decrease  of  biomass  from  5364  in  the  baseline  run  to  4920  with 
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Figure  3-28.    Baseline  simulation  of  diversity  and  production 
minimodel  SPDIV  in  Figure  3-25  for  Biosphere  2  rainforest 
for  100  years.  Final  values  for  plant  biomass  and  number  of 
species  are  on  the  right  side  of  their  graphs.  g/m2  =  grams  of 
dry  biomass  per  square  meter. 
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species  additions.  Over  100  years,  the  addition  of  50  species  in  the  second  year 
had  no  apparent  effect  on  biomass  or  number  of  species  in  the  system. 

Emergy  Evaluation  of  Rainforest  in  Biosphere  7 

The  design  and  construction  phase  of  Biosphere  2  lasted  from  1985 
until  the  first  closure  in  September  1991.   The  costs  accumulated  during  that 
phase  were  assembled  so  that  the  total  could  be  annualized  over  the  planned 
lifespan  of  the  Biosphere.  Figure  3-29  shows  a  summary  diagram  of  inputs  to 
the  Biosphere  2  rainforest  used  in  evaluating  the  emergy  of  system  start-up. 
The  start-up  costs  include  costs  of  construction,  intellectual  capital  invested  in 
the  design  of  structures  and  ecosystems,  the  stored  emergy  of  soil  and  biotic 
collections  made  from  both  captive  and  wild  sources,  and  the  direct  input  of 
sun  to  plants  and  of  wind  for  the  evaporative  cooling  tower  process  used  after 
the  mechanical  systems  were  in  operation.   Table  3-25  shows  the  values  for 
each  item,  with  notes  explaining  the  origin  of  the  values.   The  emergy  value 
of  design  and  construction,  which  included  all  goods  and  services  of  building 
and  designing  the  structure,  dwarfs  that  of  the  collection  of  ecosystem 
components  during  the  start-up  phase. 

The  emergy  output  of  the  Biosphere  2  system  is  information  and 
tourism.   Information  is  in  the  form  of  published  scientific  results,  media 
coverage,  web  page  usage,  and  the  undergraduate  program  through  Columbia 
University  that  has  been  in  place  since  1997.     Tourism  includes  the  visitation 
by  a  certain  number  of  people  yearly  to  see  the  Biosphere  and  learn  about  how 
the  earth  works. 
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Table  3-25.    Accumulated  emergy  inputs  to  Biosphere  2  rainforest  for  start-up 
of  the  system  prior  to  material  closure  in  1991. 


Note     Item 

Data  &  Units 

Solar 
Transformity 

Solar  EMergy 
E13sej 

EM$va 
lue 

Environmental  Sources 

1  Sun 

2  Wind 

5.66  E12  J 
1.75  E14  J 

lsej/J 
1.5E3sej/J 

0.566 
26,300 

3.45 
160,365 

Ecosystem  Components 

3  Plants  at  closure 

4  Plant  collection 

4.18  E10  J 
6.0  E5$ 

1.63E7sej/J 
1.64E12sej/$ 

68,134 
98,400 

415,451 
600,000 

5         Soil,  mineral 
fraction 

4.78  E9g 

1.0E9sej/g 

47,762 

291,231 

6         Soil,  organic 
fraction 

1.32  E12  J 

7.4E4sej/J 

9,768 

59,561 

Design,  construction  and  < 
7        Design  and 
construction 

aperations 

22.5  E6  $ 

1.64  E12 
sej/$ 

3,690,000 

22,500,000 

8         Electricity 

5.4  E13J 

2.0  E5 

1,080,000 

6,585,366 

Rainforest  is  approx.  15%  (1900  m2/12,766  m2)  of  the  total  surface  area  and  17%  (34,690 
m  /204,000  m3)  of  the  total  volume  of  Biosphere  2.  Transformity  values  from  Odum  (1996). 

Environmental  Sources 

1  Sun 

Average  outside  insolation  for  Tucson  is  5200  kcal/m2/day  (Romer  1985). 
Approximately  50%  of  the  outside  light  enters  the  Biosphere  and  approximately 
50%  of  the  sun  was  intercepted  by  plant  biomass.  The  rainforest  biome  is  1900  m2. 
Planting  began  about  1.5  years  before  the  1991  closure. 

(5200  kcal/m2/day)(.5)(.5))(1900  m2)(1.5  years)(365  days/yr)(4184  J/kcal)  = 
5.66  E12  J 

2  Wind 

Wind  contributed  3.37  E14  J/yr  of  kinetic  energy  for  evaporative  water  cooling 
external  to  the  Biosphere  (SBV  data  cited  in  Engel  1994).  Cooling  began  in  Sept  1989 
Wind  energy  assumed  to  have  been  contributed  in  proportion  to  volume.    Solar 
transformity  for  wind  from  Odum  (1996). 

(3.37  E14  J/yr)(2  yrs)(.26)  =  1.75  E14  J 
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Table  3-25-continued. 


Ecosystem  Components 

3  Rainforest  plants 

Biomass  at  closure  was  approximately  2500  kg  dry  weight  (Bierner  (1994)  estimate  for 
July,  1991). 

(2500  kg*1000  g/kg)(4  kcal/gm)(4184  J/kcal)  =  4.18  E10  J 

4  Plant  collection 

Emergy /money  ratio  for  1986-1991.  The  cost  of  collections,  including  labor,  transportation, 
and  permits,  was  approximately  $600,000.  Average  $/sej  ratio  for  the  years  1987-1991  is 
1.64  E12  sej/$  (Odum  1996). 

($6.0  E5)(1.64  E12  sej/$)  =  9.84  E17  sej 

5  Soil,  organic  fraction. 

Transformity  of  topsoil  organic  matter  =  7.4  E4  sej /J  (Odum  1996).  Average  organic  matter 
content  of  topsoil  is  3%  (Scott  1999).  Total  amount  of  topsoil  in  rainforest  is  1766  cubic 
meters.  Avg.  bulk  density  of  topsoil  =  1.1  g/cm3. 

(.03)(1766  E6  cm3)(l.l  g/cm3)(5.4  kcal/g)(4184  J/kcal)  =  1.32  E12  J 

6  Soil,  mineral  fraction 

Solar  transformity  for  world  sedimentary  cycle  is  1.0  E9  sej/g  (Odum  1996).  Bulk  density 
for  subsoil  is  1.43  g/cm3  and  for  topsoil  =  1.1  g/cm3  (Scott  1999).  Volume  of  subsoil  is  3340 
cubic  meters  and  for  topsoil  is  1766  cubic  meters  (Scarborough  1994).  Mineral  fraction  of 
topsoil  is  97%. 

(3340  E6  cm3)(1.43  g/cm3)  +  (0.97)(1.1  g/cm3)(1766  E6  cm3)  =  4.78  E9  g 

Design,  construction  and  operations 

8       Overall  design,  construction  and  operation  prior  to  1991  closure 

Total  cost  for  Biosphere  2  of  design,  construction  and  operation  prior  to  1991  closure  was 
$150,000,000  (SBV,  personal  communication).  The  rainforest  surface  area  is  approx. 
15%  of  the  total  Biosphere  2  area.  Average  $/sej  ratio  for  the  years  1987-1991  is 
1.64  El 2  sej/ $  (Odum  1996). 

($150  E6)(.15)(1.64  E12  sej/$)  =  3.69  E19  sej 

10     Electricity 

Electrical  consumption  for  Biosphere  2  is  approximately  5  E6  kWh/yr.  The  energy  center 
was  supporting  the  Biosphere  for  3  years  prior  to  closure  in  1991. 

(5  E6  kWh/yr)(3.6E6  J/kWh)(3  yrs)  =  5.4  E13  J 


A  comparison  of  emergy  in  the  developing  rainforest  in  Biosphere  2, 
including  only  start-up  costs  incurred  before  September  1991,  and  a  25  year  old 
natural  tabonuco  forest  are  shown  in  Table  3-26. 
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Table  3-26.  Comparison  of  emergy  in  developing  rainforest  in  Biosphere  2  i 
1991  and  25  year  old  natural  tabonuco  forest. 

Note        Location  Emergy/ha         Em$a/ha  $  cost/ha 

1  Biosphere  2  2.64  E20  230,000,000  (147,400,000) 

2  Tabonuco  forest  1.15  E17  100,152 


4050 


aEmergy/(1.15  E12  1995  sej/$) 
Notes 

1  See  Table  3-27. 

2  Odum  et  al.  (1999). 


CHAPTER  4 
DISCUSSION 

Understanding  the  unusual  process  of  ecological  succession  with 

declining  species  in  Biosphere  2  was  studied  considering  the  species  becoming 

established,  species-individual  relationships,  models  of  mechanisms,  and 

comparisons  with  early  succession  in  the  rainforest  at  El  Verde,  Puerto  Rico. 

Succession  with  Declining  Diversity 

The  self-organization  of  a  rainforest  ecosystem  in  Biosphere  2  was  the 
combination  of  human  efforts  to  plant  a  rainforest  and  patterns  of  self- 
organization  by  the  developing  ecosystem.   The  process  resembled  previous 
efforts  to  establish  small  aquatic  ecosystems  in  microcosms  by  seeding  with 
many  more  species  than  could  be  supported,  letting  the  system  develop  a 
compatible  set  by  processes  of  reinforcement,  selection,  and  extinction. 

Multiple  seeding  of  more  species  than  could  be  supported  provided  a 
way  to  find  the  carrying  capacity  for  diversity  and  its  limits.  The  biota  of  the 
Biosphere  2  rainforest  was  assembled  at  the  start  in  functional  groups.  In  the 
case  of  the  plants,  species  were  chosen  for  the  following  attributes:  early 
shade  producers,  food  web  species  support,  potential  foods  or  medicines, 
evolutionary  interest,  shade-tolerant  canopy  trees,  and  others.   In  the  case  of 
animals,  functional  groups  included  the  following:   food  web  complexity, 
detritivores,  pollinators,  and  prey  items.   Opportunities  for  a  unique 
ecosystem  existed  as  a  result  of  the  assembly  of  plants  and  animals  from 
different  regions  of  the  world,  and  the  packing  of  species  into  an  area  that  in 
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nature  would  not  normally  support  such  high  diversity.   In  the  Biosphere  2 
rainforest,  the  rate  of  species  decline  was  expected  to  decrease  as  food  webs 
became  organized  and  vertical  forest  structure  matured.  For  example,  species 
that  could  not  survive  the  open  canopy  of  the  first  year  of  plant  succession 
may  have  survived  during  later  succession,  when  a  more  continuous  canopy 
and  more  consistent  humidity  were  established.   The  vascular  epiphytes  were 
an  example. 

Effect  of  First  Arrivals 

First  arrivals,  or  colonizing  species,  in  natural  system  may  have  either 
positive,  negative,  or  neutral  net  effect  on  later  species  colonization  (Connell 
et  al.  1987).  Walker  et  al.  (1996)  noted  that  first  arrival  on  landslides  may  slow 
the  rate  of  plant  succession,  largely  due  to  effects  of  shading.  This  was  noted 
especially  for  a  landslide-colonizing  fern  that  formed  dense  thickets  in  high- 
light areas. 

Comparison  of  survival  data  over  the  two  planting  dates  in  the 
Biosphere  2  rainforest  suggests  that  individuals  planted  earlier  (1990-1991) 
had  a  survival  advantage  over  those  planted  later  (1993-1994).   Between  1993 
and  1996,  75%  of  the  individuals  from  the  second  group  (total=339)  died, 
whereas  only  39%  of  the  first  group  individuals  (total=872)  died;  and  about 
60%  of  the  species  from  the  second  group  (total=92)  and  30%  of  the  species 
from  the  first  group  (total=192)  went  extinct  over  the  same  period.   An 
exception  was  epiphytes,  where  the  second  group  had  the  higher  survival 
rate.   Mueller-Dombois  and  Ellenberg  (1974)  have  noted  that  vascular 
epiphytes  are  usually  the  last  plant  growth  form  to  establish  in  a  tropical 
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rainforest,  presumably  due  to  their  need  for  structure  providing  moist 
microsites  (Kohyama  1997). 

Limited  Access  to  Gpnps 

Most  systems  do  lose  species,  but  they  are  constantly  replaced  by  other 
species.   Though  the  number  of  species  may  appear  constant,  the  composition 
of  species  may  change. 

With  succession  in  a  closed  system,  the  species  will  not  be  replaced  by 
species  from  outside  the  present  system.  Unless  there  are  species  with 
unusual  plasticity  of  genetic  change,  the  organization  of  the  species  that  are  in 
place  will  not  be  improved  by  addition  of  genetic  information  above  the 
microorganism  scale.   Selection  by  competition  will  eliminate  species  but 
there  will  be  no  new  species  unless  planned  and  implemented  by  managers. 
Without  seeding,  plant  diversity  in  the  Biosphere  2  rainforest  can  either 
decrease  or  stay  the  same. 

Figure  3-5  showed  the  decline  in  species.  The  3  survey  results  over  a  5 
year  period  showed  plant  species  decreasing  through  time  as  a  log  function  of 
cumulative  individual  plant  loss.   A  straight  line  in  a  semi-log  graph  means 
that  the  rate  of  loss  was  proportional  to  the  number  remaining,  a  linear  rate 
of  extinction.   This  suggests  that  the  species  losses  at  first  did  not  involve 
many  interactions  or  organization. 

Declining  Diversity  and  a  gggcjgs  Uatsaii 

As  expected  and  as  shown  in  Figures  3-  5  and  Tables  A-l  and  A-2,  the 
diversity  in  the  Biosphere  2  rainforest  declined  from  the  initial  condition  in 
which  an  excess  of  plant  species  was  introduced.  Many  of  the  species  were  not 
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adapted  to  the  unique  conditions  of  soil  (amended  desert  grassland  soil), 
atmosphere  (elevated  COz  and  lowered  Oz) ,  and  rainwater  composition 
(varying  salinity  and  nutrients).     Also,  after  time,  the  Biosphere  supported 
few  animals  to  aid  pollination  and  seeding.   Although  no  example  of  an 
isolated  terrestrial  mesocosm  was  found  to  suggest  the  species  diversity  that 
could  prevail,  studies  of  sealed  aquatic  microcosms  showed  that  a  few  species 
were  supported  for  years  in  small  volumes  (Folsome  and  Hanson  1986). 
Some  models  of  species  would  predict  no  species  surviving  without  some 
introductions,  but  the  microcosm  studies  suggest  that  there  are  species 
adapted  to  survive  indefinitely  in  low  diversity  situations.   Figure  3-5  shows 
the  rate  of  extinction  is  slowing  down.   The  central  question  is  when  would 
the  species  decline  stop  and  how  many  species  would  be  sustained. 

Imagine  if,  as  in  recent  management  of  the  Biosphere  2  rainforest,  no 
species  were  added.  Where  would  extinction  stop?  Hubbell  (1979)  described 
single  species  dominance  when  there  is  no  addition  of  species.   Experience 
with  sealed  microcosms  over  long  periods  of  time  does  not  reveal  cases 
where  life  disappeared. 

After  initial  organization  to  maximize  energy  intake  from  its  limited 
gene  pool,  the  Biosphere  2  system  can  go  no  further  in  addition  of  plant 
species  without  additional  seeding.   Evolution  among  the  vascular  plants 
may  not  occur  at  all  over  the  100-year  life  of  the  Biosphere.   However, 
evolution  could  take  place  in  microorganisms,  affecting  the  organization  and 
feedbacks  of  the  overall  ecosystem.    For  example,  two  previously  undescribed 
species  were  described  from  the  Biosphere,  an  amoeba  in  the  ocean  and  a 
nematode  in  the  agriculture  area.    It  is  possible  that  these  species  existed 
elsewhere  before  their  discovery  in  Biosphere  2,  but  were  able  to  achieve 
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noticeable  densities  inside  the  Biosphere  in  the  unusual  conditions  there  that 
had  not  been  seen  before. 

For  the  plants,  however,  further  self-organization  required  building 
forest  structure  with  fewer  species  and  limited  reproduction.   Growth  in  the 
Biosphere  2  rainforest  may  build  similar  age  structure  to  El  Verde,  thus  the 
number  of  individuals  within  similar  areas  may  be  similar.   Trees  in  the 
Biosphere  2  rainforest  will  continue  to  close  the  canopy,  shading  the  plants 
beneath.   The  decline  in  species  may  be  the  result  of  both  senescence  and 
failure  to  reproduce. 

Figure  4-1  compares  the  patterns  for  species  and  individuals  at  the  start 
in  1991  with  the  condition  in  1998  (Figures  3-3  and  3-4).  By  1998  the  number 
of  individuals  from  the  original  plantings  of  1890  and  339  species  had 
decreased,  but  new  individuals  had  developed  mostly  from  vegetative 
propagation  rather  than  seeds.   Defining  an  individual  as  a  new  stem  or  stem 
cluster,  there  were  about  90  species  per  1000  individuals  based  on  the  1998 
counts.  As  Figure  4-1  shows,  the  self-organizational  process  had  a  typical 
species-individual  curve  made  up  of  original  plants  and  the  asexually 
developed  plants. 

An  analysis  follows  of  the  species  remaining  to  estimate  the  possible 
long  term  species  diversity  that  could  be  sustained  under  the  1998  conditions 
where  the  animals  needed  for  normal  propagation  were  absent. 

Extrapolation  pf  Species  Composition  Arr0rdinfT  to  RpprnH^ti™ 

An  examination  of  species  growing  in  the  rainforest  in  1998  showed 
some  species  reproducing,  some  not  reproducing,  and  many  spreading 
clonally.  Based  on  the  survey  in  Table  B-2,  the  list  of  plant  species  was 
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Figure  4-1.  Extrapolation  of  number  of  plants  in  El  Verde 
rainforest  to  a  possible  diversity  plateau  (solid  triangle)  for 
the  Biosphere  2  rainforest.  Measured  data  from  Biosphere  2 
are  shown  with  solid  circles  (1991)  and  open  squares  (1998). 
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modified  to  predict  species  that  might  be  part  of  the  eventual  composition  if 
there  is  little  further  seeding  from  outside  and  all  the  non-reproducers  and 
short-lived  species  have  lived  and  died. 

Some  of  the  species  expected  to  live  for  100  years  are  the  trees  that  form 
the  canopy  layer.  Those  that  do  not  senesce  earlier  will  be  part  of  the  later 
flora.  Cecropia,  for  example,  is  part  of  the  current  canopy  but  will  eventually 
either  be  overshaded  by  other  species  or,  since  it  is  a  short-lived  species,  it  will 
die.   Additionally,  several  Cecropia  treefalls  have  occurred  during  the  first  3 
years  of  forest  development,  but  the  fallen  trees  sprouted.   If  sprouting  occurs, 
Cecropia  will  persist  as  part  of  the  flora.   Other  species,  such  as  Heliconia  and 
Strelitzia,  may  flower  yet  not  set  seed,  propagating  by  extending  underground 
organs  which  sprout  culms  a  distance  from  the  original  plant.   The 
extrapolated  composition,  Table  4-1,  shows  a  potential  species  pool  of  30 
species  for  the  0.19  ha  area.  See  the  species-individual  curve  (Figure  4-1)  for 
an  extrapolation  of  the  number  of  individual  plants. 

Simulation  of  Species  Derlinp 

A  simulation  of  the  diversity  decline  was  made  with  the  model  in 
Figure  3-25  with  a  one-time  addition  of  50  species  after  2  years  of 
development.  The  model  was  run  for  10  years  (Figure  D-3).  After  10  years, 
diversity  was  increased  over  the  baseline  simulation  with  no  additions.    The 
addition  had  no  effect  on  diversity  at  100  years.    A  simulation  with  yearly 
additions  of  20  species,  Figure  C-2,  resulted  in  175  species  at  100  years 
compared  to  the  minimum  species  number  of  30  established  when  no 
additions  were  made.   Extrapolations  made  from  the  actual  measurements 
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Table  4-1.  Extrapolation  of  species  composition  of  Biosphere  2  rainforest  to 
100  years  with  30  species.  H=Herb,  T=Tree,  P=Palm,  G=Giant  herb, 
A=Bamboo,  C=CIimber,  S=Shrub. 


Family 

Species 

Growth 
form 

Reproductive  mode 
in  Biosphere  2 

Amaryllidaceae    tucnaris  grandiflora 

H 

Glonal 

Araceae 

Colocasia  sp. 
Dieffenbachia  sp. 
Monstera  deliciosa 
Scindapsus  aureus 
Zamioculcas  zamiifolia 

H 
H 
C 
C 
H 

Clonal 
Clonal 
Clonal 
Clonal 

Seed 

Bignoniaceae 

Spathodea  campanulata 

T 

Long-lived 

Bombacaceae 

Ceiba  pentandra 
Pachira  aquatica 

T 
T 

Will  survive  if  long-lived 
Seed,  needs  high  light  to 

Cecropiaceae 

Cecropia  schreberiana 

T 

reproduce 
Sprouting 

Cyclanthaceae 

Carludovica  palmata 

G 

Seed 

Euphorbiaceae 

Hum  crepitans 

T 

Long-lived 

Heliconiaceae 

Heliconia  sp. 

G 

Clonal 

Malvaceae 

Hibiscus  rosa-sinensis 

S 

Clonal 

Moraceae 

Ficus  pumila 

c 

Long-lived 

Musaceae 

Musa  textilis 
Musa  paradisiaca 

G 
G 

Clonal 
Clonal 

Phytolaccaceae 

Phytolacca  dioica 

T 

Long-lived 

Poaceae 

Bambusa  multiplex 
Bambusa  tuldoides 

A 
A 

Clonal 
Clonal 

Rubiaceae 

Coffea  arabica 

S 

Seed 

Solanaceae 

Brugmansia  suaveolens 

S 

Clonal,  sprouts 

Strelitziaceae 

Phenakospermum  guyanense 
Strelitzia  reginae 

G 
G 

Clonal 
Clonal 
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Table  4-1-continued. 


Family 


Vitaceae 

Zamiaceae 

Zingiberaceae 


Species 


Cissus  sia/oides 

Zamiafurfuracea 

Alpinia  purpurata 
Costus  sp. 

Hedychium  coronarium 
Zingiber  spectabile 


Growth      Reproductive  mode 
form        in  Biosphere  2 


Clonal 

Long-lived 

Clonal 
Clonal 
Clonal 
Clonal 


also  showed  that  the  addition  of  50  species  after  2  years  will  have  little  effect 
on  the  diversity  of  the  rainforest  at  100  years. 

The  turn-over  times  of  species  within  geographic  areas  varies  with 
internal  and  external  processes.     The  extinction  rate  declined  over  time 
(Figure  3-5),  from  about  60  species  per  year  from  1991-1993  to  about  30  species 
per  year  from  1993-1996.  A  total  of  366  species  were  added  during  that  period, 
averaging  73  per  year.  Species  seeding  on  a  continuous  basis  may  be  required 
to  maintain  the  diversity  of  the  system  above  the  approximately  30  species 
that  are  expected  to  survive.  Most  of  the  additions  may  be  considered  the 
"transient  species"  in  the  model  (Figure  3-25). 

Mechanisms  Affecting  Diversity  in  Biosphere  7 

In  the  introduction  many  mechanisms  affecting  diversity  were  cited 
from  published  literature.   The  following  seemed  to  be  operating  in 
Biosphere  2. 

Excess  Resources  and  Diversity 

When  nutrients  and /or  light  are  in  excess,  diversity  is  often  reduced  by 
competitive  overgrowth  by  a  few  species.  This  concept  was  applied  to  the 
ecosystem  in  Silver  Springs  by  Yount  (1956).   For  example,  Kent  (1996)  found 
species  diversity  lower  in  fertilized  pond  plots  in  Puerto  Rico  over  his  clear 
water  controls,  and  vegetation  came  in  more  quickly  to  the  fertilized  plots. 
Odum  and  Odum  (in  press)  used  a  model  NUTRSPEC,  and  Odum  (1999)  used 
a  more  complex  version  PIONINFO  to  show  the  mechanisms. 

In  1991-1995  in  Biosphere  2  one  of  the  most  obvious  resource  excesses 
was  C02.  Growth  of  'weedy'  species,  and  in  particular  vines  was  observed, 
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helping  sequester  carbon  dioxide  rapidly.  These  species  were  aided  by  the 
large  amount  of  structure  for  climbing  provided  by  the  stainless  steel 
spaceframe.  The  result  was  a  dense  weedy  canopy  above  the  upper  tree 
canopy  of  the  rainforest.  This  diversity-reducing  factor  was  at  least  partly 
counteracted  by  the  frequent  pruning  of  the  weeds  by  the  author  and  others. 

Effect  of  Trimming  Weedy  Growth 

A  major  effort  was  made  to  trim  out  rapidly  growing  vines  and  other 
successional  weeds  in  Biosphere  2  during  the  first  two  years  after  seeding  and 
thereafter.    During  the  first  two  years,  most  of  the  pruned  biomass  was 
removed  from  the  rainforest  and  stored  and  air  dried  in  the  underground 
stainless  steel  and  concrete  passageways.  A  fraction  of  the  fresh  biomass  was 
fed  to  the  domestic  livestock.   Thus,  human  effort  was  directed  to  truncate 
early  successional  processes  and  to  accelerate  the  development  of  the  climax 
state.  Pruning  was  analagous  to  a  large  herbivore  or  a  large  pulse  of 
herbivorous  insects.   The  total  number  of  person  hours  pruning  was  about  22 
hours  per  month  during  the  first  2-year  closure. 

Alternatively,  a  different  climax  state  may  have  developed  had  early 
successional  processes  not  been  arrested.    Removal  of  the  early  successional 
biomass  may  have  altered  the  relative  frequencies  of  species  in  the  rainforest. 
Those  species  that  were  adapted  to  conditions  of  higher  light  may  have  been 
favored  over  shade-adapted  species.  The  result  may  have  been  loss  of  species 
that  were  adapted  to  shading.  An  extreme  possibility  is  that  over  a  longer 
period,  when  the  light-adapted  species  formed  a  continuous  canopy,  no 
species  would  remain  to  grow  in  the  shaded  conditions. 
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Additionally,  removal  of  biomass  that  may  have  contributed  to  the 
litter  accumulation  could  have  an  effect  on  nutrient  cycles,  nutrient 
availability  and  productivity  of  the  rainforest  as  has  been  suggested  for  the 
rainforest  of  Puerto  Rico  (Lodge  et  al.  1991). 

An  effect  of  removing  the  two  most  abundant  species  from  the  data 
and  replotting  the  graph,  Figures  3-12,  is  a  higher  species  diversity  number  of 
individuals  counted.  The  species  diversity  graphs  during  this  were  steeper 
than  they  might  have  been  without  the  pruning,  partly  reversing  the  effect  of 
high  nutrient  overgrowths. 

The  production  and  biomass  model  in  Figure  3-19  was  run  with  and 
without  pruning  (Figure  3-23)  to  see  whether  possible  effects  of  pruning  were 
indicated.     The  simulation  results  of  the  minimodel  suggest  that  without 
pruning,  diversity  would  have  been  reduced  through  loss  of  mature  species. 
Over  100  years,  no  mature  species  remained  and  the  simulated  system  was 
composed  of  only  weedy  species. 

Pulsing  Disruption 

Many  ecologists  offer  the  theory  that  disruption  increases  diversity  by 
preventing  competitive  exclusion  overgrowths  between  species.    Where 
disruptions  create  open  patches  of  successional  vegetation  in  formerly  closed 
canopy  forest,  diversity  is  augmented  by  the  successional  colonizers. 
Hurricanes  and  landslides  at  El  Verde  create  patches  in  various  seasons. 
Disturbances  in  Biosphere  2  include  those  by  pruning,  trampling  by  visitors, 
and  the  shift  to  lower  C02  regime  with  throughflow  of  air  in  1995.   Some 
shocks  obviously  reduce  diversity  such  as  transplant  shock  and  excessive  heat 
on  the  tops  of  trees  in  summer. 
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A  more  recent  management  strategy  (starting  in  1997)  for  the  Biosphere 
2  rainforest  has  been  to  trim  trees  that  are  touching  the  glass  boundary  and  to 
deposit  all  of  the  biomass  on  the  forest  floor.  Pruning  is  scheduled  cyclically, 
and  creates  a  pulse  of  organic  matter  added  to  certain  areas.  An  interesting 
question  would  be  what  effect  the  pulsed  pruning  and  pulsed  biomass 
deposition  has  on  forest  productivity,  and  to  what  extent  this  would  vary 
from  normal  leaf  and  branch  fall  and  litter  accumulation. 

Consumers 

Several  hundred  insect  species  were  either  intentionally  introduced, 
accidentally  introduced,  or  became  established  in  Biosphere  2  before  the 
structure  was  completed.   By  1997,  the  most  conspicuous  insects  in  Biosphere 
2  were  ants,  homopterans,  katydids,  and  cockroaches  (Wetterer  et  al.  in  press) 
The  ecologically  dominant  species  is  the  ant  Paratrechina  longicornis,  a 
species  which  has  spread  from  the  Old  World  throughout  disturbed  areas  in 
the  tropics  and  subtropics.  It  is  the  same  tramp  ant  species  that  has  taken  over 
nests  normally  occupied  by  native  species  in  the  Dry  Tortugas  in  the  Florida 
Keys  (Wetterer  et  al.  in  press). 

The  lack  of  herbivorous  insects  is  shown  by  the  low  percentage  of  leaf 
holes  in  Biosphere  2  compared  to  the  tabonuco  forest  (Table  3-21).  Canopy 
arthropods  may  have  a  significant  effect  on  primary  productivity,  nutrient 
cycling  and  hydrology. 

Graphical  Representation  of  Cumulative  SppHPs-Individnal  fnnntc 

To  help  understand  patterns  of  diversity,  species-individual  counts 
were  represented  on  several  coordinates   with  linear  scales,  with  semi- 
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logarithmic  plot,  and  double  logarithmic  plot.   See  the  example  of  the  study 
sites  in  Figures  3-18a-c.  Consider  the  kind  of  graphs  that  some  hypotheses 
might  predict. 

First,  suppose  there  were  many  species  available  for  seeding  in  a  large 
area  so  that  each  individual  added  was  a  different  species.  In  that  case  the 
species  added  would  increase  in  a  straight  line  on  the  linear  plot,  but  curve 
sharply  up  on  the  semi-log  plot. 

Next,  suppose  there  was  a  great  excess  of  individuals  of  one  or  two 
species,  a  situation  sometimes  occurring  where  nutrients  are  in  excess  and 
there  is  competitive  exclusion  or  the  seeding  from  a  low  diversity  source.   In 
that  case  the  graph  will  be  shifted  to  the  right,  and  will  increase  very  slowly. 

Third,  suppose  that  there  are  interactions  among  species  so  that  the 
species  are  being  organized  in  systems  relationships  with  resource  use  of  the 
organization  increasing  as  the  square  of  the  number  of  kinds.   In  that  case  the 
number  of  individuals  required  to  support  each  additional  species  increases 
as  the  square  of  the  species.  Lines  are  added  to  the  graphs  to  show  the 
situation  where  S  is  proportional  to  N2  or  conversely  that  N  is  proportional 
to  the  square  root  of  species  S.  Hierarchical  relationships  are  suggested  when 
there  are  many  of  a  few  species  and  few  of  many  species. 

Fourth,  loss  of  pollinators  in  a  system  would  show  a  loss  of  certain 
species  that  can  no  longer  reproduce.   The  number  of  individuals  of  species 
that  do  not  require  pollinators  or  of  clonal  species  would  increase,  reducing 
the  number  of  species  while  not  necessarily  reducing  the  number  of 
individuals  where  an  individual  was  defined  as  a  separate  stem. 

Another  characteristic  curve-shape  reflects  extremes  of  available  energy 
such  as  in  both  shallow  tropical  seas  (high  energy  inputs  from  currents, 
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sunlight,  etc.)  and  deep  ocean  bottoms  (low  energy  inputs),  both  with  many 
species  and  individuals  per  1000  count,  but  in  the  former  they  are 
concentrated  in  less  area  compared  to  the  latter,  where  there  is  a  larger 
distance  between  individuals. 

The  shapes  of  curves  of  cumulative  species  versus  cumulative 
individuals  observed  in  Biosphere  2  and  in  the  colonization  in  the  rainforest 
were  compared  with  the  theoretical  patterns.  Initial  patterns  of  seeding  in 
Biosphere  2  and  in  natural  seeding  of  landslide  openings  within  the  complex 
forest  turned  up  on  semi-log  plot  suggesting  less  organizational  interactions. 
Slopes  of  species  curves  from  the  natural  forest  were  not  so  steep  or  declining 
with  increased  individuals-suggesting  more  hierarchy  or  organization. 

If  when  species  are  first  added,  there  is  one  new  species  for  each 
individual,  the  line  is  straight  on  the  linear  plot  and  sharply  curving  upward 
on  semi-log.   As  more  individuals  are  seeded,  duplications  occur  either  on 
purpose  as  in  Biosphere  2  or  out  of  probabilities  of  seeding  access  as  in  the 
landslides.  The  curve  bends  to  the  right  on  linear  plot  and  develops  as  a 
straight  zoned  beyond  the  first  few  individuals  on  a  semi-log  plot. 

The  effect  of  many  individuals  of  a  few  species  was  demonstrated  by 
removing  two  abundant  species  and  replotting  (Figure  3-12).  The  effect  is  to 
shift  the  curve  to  the  left.  Adding  an  excess  of  species  that  dilute  the  others 
shifts  the  straight  zone  of  the  graph  to  the  right.  In  the  Radiation  Center  23 
years  after  radiation  (before  hurricanes)  there  was  a  predominance  of  a  few 
species  with  many  seedlings  making  the  linear  curve  flat  and  the  semilog 
irregular. 

The  graphs  for  Biosphere  2  were  compared  with  the  El  Verde  and 
Bisley  graphs  used  as  a  reference  of  normal  organization  of  diversity.   In  the 
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semi-log  graph,  the  Biosphere  2  curve  turns  up  at  the  end  whereas  the 
natural  forest  line  is  straighter.   This  curvature  can  be  interpreted  as  an 
indicator  that  the  species  are  less  organized.  See  Table  4-2  for  comparison  of 
the  slopes  on  the  double  log  plot. 

Effect  of  Spatial  Characteristics  of  Plants  in  BinsphprP  7 

With  the  exception  of  Leucaena  spp.  trees,  the  plants  at  the  time  of 
placement  in  the  Biosphere  2  rainforest  were  between  approximately  0.5 
meters  and  3  meters  tall.   A  comparison  of  the  horizontal  distribution  of 
plants  in  their  initial  configuration  suggests  a  regular  distribution  of  plants 
(Figure  3-10).  This  departure  from  successional  pathways  that  are  found  in 
nature  -  beginning  with  a  seed  pool  rather  than  uneven  size  classes  of  plants 
and  regular  rather  than  clumped  distributions  ~  could  have  an  effect  on 
survival  of  species  as  related  to  shading.   In  particular,  species  that  were  not 
shade-tolerant  may  have  been  planted  in  the  shade  of  other  species,  without 
an  opportunity  to  overgrow  other  species.   Resulting  composition  of  the 
rainforest  may  have  been  shaped  by  such  differences. 

Com  pari  sons  with  El  VernV 

With  the  passage  of  time,  the  rainforest  in  Biosphere  2  developed 
properties  similar  to  those  at  El  Verde.   There  were  similarities  between  early 
self-organization  and  the  successional  stages  at  El  Verde. 
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CQmparispns  Between  Rainforest  Structure  in  Biosphere  ?  and  El  Verde 

A  visual  survey  of  the  Biosphere  2  and  El  Verde  rainforests  showed 
differences  and  similarities.   In  the  tabonuco  forest,  a  deep  organic  mat  of  fine 
roots  forms  above  the  heavy  clay  soil.  In  Biosphere  2  the  roots  were  deeper  in 
the  ground  with  no  mat  at  the  surface.   Previous  authors  accounted  the  mat 
of  small  roots  at  El  Verde  as  an  adaptation  to  low  transpiration  rates  and  low 
soil  nutrients  (Odum  1970).   Saturation  deficit  aiding  transpiration  was  6  to  9 
millibars  in  Biosphere  2  compared  to  1.8  to  6  mb  at  El  Verde.   Nutrients  were 
higher  in  Biosphere  2,  especially  phosphorus.   Root  adaptations  were 
consistent  with  the  theory.   Wind  velocity  was  lower  in  the  Biosphere  2 
rainforest  breeze.    Another  striking  difference  was  the  greater  density  of 
certain  growth  forms,  and  in  particular  of  epiphytes,  both  vascular  and  non- 
vascular, in  the  tabonuco  forest  (Figure  3-14).   Most  of  the  vascular  epiphytes 
had  perished  from  the  Biosphere  2  rainforest  by  Sept.,  1993.  There  were  no 
mosses,  liverworts,  and  filamentous  algae  species  covering  leaves,  trunks  and 
rocks. 

Soil  StrnrhirP 

The  soils  in  the  Biosphere  2  rainforest  had  a  coarser  texture  than  the 
clays  and  organics  at  El  Verde.  Scott  (1999)  found  some  differentiation  of  soil 
profiles  in  his  1993  measurements.   Earthworms  were  processing  clays  and 
organic  matter.  However  the  pH  was  still  above  7  below  the  litter,  compared 
to  the  normally  acid  conditions  in  El  Verde.   Whereas  soil  was  generating 
carbon  dioxide  from  high  levels  of  organic  matter,  some  was  being  absorbed 
by  the  basic  substances. 
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In  a  study  by  Gonzalez  et  al.  (1999),  the  tabonuco  forest  was  compared 
with  wet  areas  dominated  by  palms  and  Heliconia.  Litterfall  in  the  two  areas 
studied  was  similar.  The  tabonuco  had  lower  pH,  with  a  mean  value  of  4.8 
(4.5  -  6.5),  and  twice  as  many  earthworms.   The  palm  area  had  more  soil 
water  and  higher  pH  (5.1-5.9).    The  earthworm  Pontoscolex  corethrurus  was 
dominant  in  both,  with  20  to  160  grams/m2  and  30  to  400  individuals /m2. 

Absence  of  Most  Animals  and    Simplification  of  Food  Wph 

Most  of  the  animals  of  typical  rainforests  either  died  or  were  never 
seeded  in  Biosphere  2.  The  hundreds  of  insect  species  at  EI  Verde  were  absent 
from  Biosphere  2,  yet  decomposition  was  slightly  faster  in  Biosphere  2 
(Nelson  1999).   Litter  fall  developed  values  like  those  at  El  Verde  and  organic 
differentiation  of  the  soil  was  accelerated  (Scott  1999).   Earthworm 
populations  were  present  and  the  numerous  tiny  ants,  Paratrechina 
longicornis  (Latreille),  became  the  dominant  ant  species. 

The  effect  of  the  lack  of  animals  in  the  Biosphere  2  rainforest  can  be 
related  to  the  possible  effects  of  the  loss  of  animal  diversity  in  increasingly 
fragmented  ecosystems. 

Plant  Reproduction 

Minimal  plant  reproduction  by  seed  has  occurred  in  the  Biosphere  2 
rainforest  since  the  1991  closure.  Since  most  potential  pollinators  perished 
during  the  first  2-year  closure,  a  lack  of  pollen  dispersal  vectors  along  with 
low  light  and  other  environmental  factors  would  prevent  both  pollination  of 
flowers  or  flowering  of  plants.  Most  of  the  species  that  reproduce  sexually 
were  located  in  sunnier  locations,  such  as  the  ginger  belt.   Species  that  fruited 
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at  first  stopped  producing  after  they  were  shaded  later.  The  most  conspicuous 
characteristic  of  the  forest  floor  in  Biosphere  2  was  the  many  stems  of  plants 
emerging  from  asexual  growth  from  runners  of  several  species. 

In  contrast,  most  of  the  species  in  the  El  Verde  rainforest  reproduce  by 
seed.  In  disturbed  areas,  clonal  taxa  such  as  Heliconia  may  be  part  of  the  early 
successional  flora,  but  they  are  not  a  dominant  species  of  the  forest.  The 
energy  involved  in  normal  reproduction  (formation  of  fruits)  in  the  El  Verde 
forest  was  0.075  kcal  m"2  d"1  (Odum  1970).   By  favoring  vegetative  reproduction 
in  Biosphere  2,  this  energy  was  saved  for  growth,  but  genetic  diversity  was 
lost  increasing  risk  of  epidemic  mortalities  later.  The  species  that  flowered 
without  setting  seed,  either  due  to  lack  of  pollinators  or  light  quantity  or 
quality,  were  using  energy  for  reproduction  that  will  not  be  manifest,  taking 
the  energy  from  growth. 

Species  and  Individuals 

At  El  Verde,  in  mature  forest  there  were  7063  plants  per  hectare 
(average  of  pre-radiation  center,  south  control  center  (Rushing  1970),  and 
counts  by  R.  F.  Smith  (1970)).  For  the  equivalent  area  of  Biosphere  2  (0.19  ha) 
there  were  1342  plants. 

The  species-individual  count  in  the  successional  development  of  bare 
areas  of  landslides  after  11  years  (Figure  3-16)  was  60  species  per  1000 
individuals  compared  to  about  80  species  in  Biosphere  2  in  1998  after  8  years. 
Extrapolating  the  extinction  rate,  diversity  in  Biosphere  2  will  be  similar  to 
that  in  succession  at  El  Verde  in  1999.   At  EI  Verde,  diversity  may  plateau  at  a 
lower  level,  while  developing  more  vegetative  structure  among  those  species 
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Table  4-2.  Comparison  of  diversity 
indices:  species/ 1000  individuals 
and  slope  of  log-log  graphs  (K). 


Species/ 

Site  name 

1000 

K 

Biosphere  2 

1991 

273 

0.621 

1998-subset 

115 

0.486 

1998-all 

96 

0.413 

Bisley  grid 

63 

0.364 

Landslide 

60 

0.339 

Radiation  Centei 

41 

0.276 

K=slope  of  the  graph  of  log  of  number 
of  species  as  a  function  of  log  of 
number  of  individuals. 
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that  prevail.  Figure  4-1  shows  a  possible  future  with  1342  plants  and  30 
species. 

Table  4-2  compares  Biosphere  2  with  tabonuco  forests  and  introduces 
another  diversity  index,  where  the  slope  of  the  log-log  curve  of  species/ 1000 
individuals  is  measured.   The  index  may  indicate  the  amount  of  organization 
that  exists  in  the  ecosystems  measured.  Margalef  (1997)  used  this  double 
logarithmic  index  to  contrast  low  diversity  plankton  blooms  with  mature 
high  diversity  populations. 

General  Implications 

Comparison  of  Biosphere  2  with  RoIp  and  Trends  of  Global  Diversity 

Biosphere  2  is  useful  to  help  understand  the  global  biosphere.  Earth  is 
also  a  closed  system  whereby  seeding  of  species  does  not  occur  from  outside 
the  system.   On  a  much  larger  scale,  much  like  the  Biosphere,   the  process  of 
diversity  on  earth  includes  only  extinction  and  evolution.    Because  evolution 
occurs  over  a  relatively  long  time  scale,  it  probably  will  not  occur  in 
Biosphere  2  for  macroorganisms,  but  could  feasibly  occur  for  microorganisms. 
There  is  controversy  now  in  determining  if  diversity  on  earth  is  downsizing 
as  in  Biosphere  2.   Species  extinction  may  be  occurring  faster  than  species 
evolution.  Adding  1  species  per  hectare  per  year  in  rainforest  systems  at 
EI  Verde,  Puerto  Rico  (Scatena,  personal  communication)  may  be  enough  to 
replace  species  and  maintain  the  current  diversity.    The  diversity  minimodels 
can  indicate  the  number  needed  if  the  extinction  rate  is  known. 
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Declining  Diversity  and  Downsizing 

The  decrease  of  diversity  in  Biosphere  2  may  have  parallels  with  a 
civilization  which  is  downsizing.   Many  specialized  jobs  (analogues  to 
species)  only  have  few  or  no  species  who  can  do  them,  as  a  more  generalist 
way  of  making  a  living  prevails.   In  a  civilization  which  is  downsizing  -  the 
many  specialized  jobs  now  only  have  one  person  who  can  do  them,  as  others 
have  had  to  find  another  more  generalist  way  of  making  a  living.   The  ethnic 
diversity  in  colonizing  systems  and  the  costs  of  organization  may  be  like  a 
forest. 

Comparing  Succession  of  Diversity  in  Biosphere  2  with  Global  Cultural 
Change 

On  a  different  scale  there  may  be  analogy  between  diversity  in 
ecosystems  and  in  human  cultural  change.   Diversity  includes  the  number  of 
species  and  patterns  of  their  distribution.  Higher  levels  of  organization  might 
be  expected  to  increase  ethnic  diversity  of  occupations,  cultures  and  social 
differences. 

The  Balkan  war  of  1999  started  with  separation  of  cultures  with 
destructive  competition  (ethnic  cleansing).   At  first  this  was  a  loss  of  diversity. 
The  patterns  that  developed  through  interactions  with  the  larger,  developed 
system  appeared  to  be  diminishing  competition  and  increasing  diversity  in 
1999. 
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Carrying  Capacity  of  Systems  for  Diversity 

There  is  growing  evidence  that  there  is  a  carrying  capacity  for  diversity 
for  a  given  system  at  a  given  time.    For  instance,  Ashton  (1977)  states:   "We 
now  have  growing  evidence  .  .  .  that  there  is  a  maximum  number  of  species 
that  a  forest  can  accommodate,  that  this  varies  with  site  conditions  and  that  it 
has  already  been  approached  in  west  Malesian  forests."    Diversity  indices  on 
an  area  basis  measure  the  carrying  capacity. 

Carrying  capacity  is  the  diversity  that  is  eventually  stabilized  in 
Biosphere  2.   According  to  some  models  such  as  the  one  in  Figure  3-28,  the 
diversity  achieved  coming  down  should  be  similar  to  that  achieved  by 
gradually  adding  species  -  typical  of  natural  succession,  such  as  the  landslide 
at  EI  Verde.   Even  now  the  diversities  are  similar. 


APPENDIX  A 
BIOSPHERE  2  RAINFOREST  MAPS 


Figure  A-l.  Map  cells  showing  the  location  of  individual  plants  in  the  Biosphere 
2  rainforest.  See  Table  A-l  for  corresponding  list  of  species  and  survey  of 
individuals,  a)  Key  to  the  location  of  individual  map  cells  within  the  rainforest- 
s' ii  tCnl2;udlCel\  3;  6)  Ce"  4;  f)  CdI  5;  S}  Cdl  6' h> Cel1  7> »  CeI1  8;  J)  Cell 
9;  k)  Cell  10;  1)  Cell  11;  m)  Cell  12;  n)  Cell  13;  o)  Cell  14;  p)  Cell  15;  q)  Cell  16-  r) 
Cell  17;  s)  Cell  18;  t)  Cell  19;  u)  Cell  20.  Surveying  and  original  graphics  by ' 
Teague  and  Co.,  Tucson,  AZ.  Approximate  scale:  1  cm  =  0.6  m 
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Figure  A-l.  Map  cells  showing  the  location  of  individual  plants  in  the 
Biosphere  2  rainforest.  See  Table  A-l  for  corresponding  list  of  species 
and  survey  of  individuals,  a)  Key  to  the  location  of  individual  map  cells 
T^n  ?t  !"^0reSJ;  b)  CeI1 1;  c)  CeI1 2-  d)  Ce»  3;  *)  Cell  4;  f)  Cell  5; 
8  ?    J ih)  ft",,7''  0  CdI  8;  j)  CeI1  9;  k)  Cel1 10'  »  CeI1  ";  m)  Cell  12; 
n  ?    J?  °J  Ce"  M;  P)  CeI1 15;  ^  CeI1 16'  r) CeI1 17;  s)  Cell  18;  t)  Cell  19- 
u)  Cell  20.  Surveying  and  original  graphics  by  Teague  and  Co., 
lucson,AZ.  Approximate  scale:  1  cm  =  0.6  m. 
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b)  Cell  1 
Figure  A-l  -  continued. 


c)  Cell  2 
Fieure  A-l  -  continupH 


d)  Cell  3 
Figure  A-l  -  continued. 


*3^ 


e)  Cell  4 
Figure  A-l  —  continued. 
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f )  Cell  5 
Fieure  A-l  --  continued. 


i®r~ 


*"&%" 


g)  Cell  6 
Fieure  A-l  -  continued. 
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h)  Cell  7 
Figure  A-1  -contimipH 


i)  Cell  8 
Fieure  A-l  -  continnpH 
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j)  Cell  9 
Fieure  A-l  -  contirmpH 


k)  Cell  10 
FieureA-1-conHnnpH 


I)  Cell  11 
Figure  A-1-continiiPrl 
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m)  Cell  12 
Figure  A-1-continiipH 
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n)  Cell  13 
Figure  A-1  -  conhmiPrl 


o)  Cell  14 
Figure  A-l  -  conrimipH 
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p)  Cell  15 
Fieure  A-l  -  contimiPH 
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q)  Cell  16 
Fieure  A-l  -  continued 
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r)  Cell  17 
Fieure  A-l  -  contirmpH 


s)  Cell  18 
Figure  A-l  -  continnpH 


t)  Cell  19 
Fieure  A-l  -  continnpd 


u)  Cell  20 
Figure  A-l  -  continued. 


APPENDIX  B 
BIOSPHERE  2  RAINFOREST  PLANT  LISTS 

Table  B-l.      Biosphere  2  rainforest  plant  survey  list.  Planting-   l=pre-1991- 
2=Oct.-Dec.  1993.  GF=Growth  form:  T=Tree,  S=Shrub,  C=Climber,  G=Gia'nt 
Herb,  P=ArboreaI  palm,  R=Graminoid,  A=Woody  graminoid,  such  as 
bamboo,  H=Herb,  E=Epiphyte.  A=Alive,  D=Dead.  Blank  cells  indicate 
information  available. 

Table  B-2.      Species  from  the  first  planting  of  the  Biosphere  2  rainforest,  with 
inventories  from  1991,  1993,  and  1996.  T=Tree,  S=Shrub,  C=CIimber,  G=Giant 
Herb,  P=ArboreaI  palm,  R^Graminoid,  A=Woody  graminoid,  such  as 
bamboo,  H=Herb,  E=Epiphyte. 

Table  B-3.      Species  from  the  second  planting  of  the  Biosphere  2  rainforest 
with  inventories  from  1993  and  1996,  or  self-propagated.  T=Tree,  S=Shrub' 
C=CIimber,  G=Giant  Herb,  P=Arboreal  palm,  R=Graminoid,  A=Woody 
graminoid,  such  as  bamboo,  H=Herb,  E=Epiphyte. 
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APPENDIX  C 
PRODUCTION  AND  BIOMASS  MINIMODEL 

Table  C-l.  EXCEL  spreadsheet  used  to  calculate  coefficients  for  Biosphere 
rainforest  production  and  biomass  minimodel  under  predicted  steady-state 
conditions. 

Table  C-2.  Program  in  QUICKBASIC  for  the  simulation  of  metabolism  of 
the  rainforest  in  Biosphere  2,  model  B2METAB  in  Figure  3-19. 

Figure  C-l.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with  light  increased  20% 
over  baseline.   Final  values  for  plant  biomass  are  on  the  right  side  of  their 
graphs.  g/m2=grams  per  square  meter,  dry  weight  for  biomass. 

Figure  C-2.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with  no  weedy  biomass. 
Final^  values  for  plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m2=grams  per  square  meter,  dry  weight  for  biomass. 

Figure  C-3.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  where  all  pruned  biomass  is 
put  onto  the  soil.   Final  values  for  plant  biomass  are  on  the  right  side  of  their 
graphs,  g/m  =grams  per  square  meter,  dry  weight  for  biomass. 

Figure  C-4.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB  with  no  pruning  by  humans 
Filial  values  for  plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m  =grams  per  square  meter,  dry  weight  for  biomass. 

Figure  C-5.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  human  effort  is  reduced  to 
.75  of  the  baseline    Final  values  for  plant  biomass  are  on  the  right  side  of 
their  graphs,  g/m  =grams  per  square  meter,  dry  weight  for  biomass. 

Figure  C-6.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  human  effort  is  reduced  to 
.85  of  the  baseline    Final  values  for  plant  biomass  are  on  the  right  side  of 
their  graphs,  g/m  =grams  per  square  meter,  dry  weight  for  biomass 


Figure  C-7.     One-hundred-year  simulation  of  Biosphere  2  rainforest 
production  and  biomass  minimodel,  B2METAB,  with  airflow  cut  off  from  the 
rest  of  the  Biosphere.  g/m2=grams  per  square  meter,  dry  weight  for  biomass. 
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Table  C-2.  Program  in  BASIC  for  the  simulation  of  metabolism  of  the 
rainforest  in  Biosphere  2,  Model  B2METAB.bas  in  Figure  3-19. 

Line   Code 

10  REM  B2METAB  - 18  APR  1999 

11  REM  Autocatalytic  consumers  C 

12  REM  Set  for  100  years'  run 

13  REM  started  with  high  soil  organics,  low  plants 
15  LINE  (0,5)  -  (400,90X3 

20  LINE  (0,110)  -  (400,190)„B 

25  LINE  (0,210)  -  (400,290)„B 

30  LINE  (0,310)  -  (400,350)„B 

35  LINE  (0,370M400,410)„B 

40  LINE  (0,330)-(400, 330),3 

44  REM  Outside  Sources 

45  Im  =  2600:REM  mean  insolation  about  which  sine  wave  varies 

46  Ir  =  500:REM  annual  range  of  insolation 

47  H  =  l 

48  Ja  - 1 

50  REM  Initial  storages 

55  W  =  100:  REM  Calibrated  at  1000 

60  M  -  900:REM  Calibrated  at  9000 

70  N  =  18:  REM  Calibrated  at  18 

80  S  =  28400!:REM  Calibrated  at  28400 

90  B  -  0:REM:  calibrated  at  23000 

95  C  =  .1:REM:  calibrated  at  1 

100  REM  COEFFICIENTS 

110  Kl  =  .00005 

120  K2  =  .  00005 

130  K3  =  6.41E-08 

140  K4  =  2.564E-07 

150  K5  =  9.853E-08 

160  K6  =  2.564E-07 

170  K7  =  . 00015 

180  K8  =  .0003 

190  K9  =  .0000732 

200  K10  =  .00015 

210  Kll -.00015 

220  K12  =  .0003 

230  K13  =  .0003 

240  K14  =  .0002778 

250  K15  =  . 0002778 

260  K16  =  .0004611 

270  K17  =  .0001122 

280  K18  =  .  00015 


Table  C-2  -  continued. 
Line   Code 
290    K19  =  l 
300    K20  =  l 


310  K21  =  .055556 

320  K22  =  6.866E-05 

325  K23  =  1.675E-05 

330  K24  =  4.08E-08 

335  K25  =  3.341E-08 

340  K26  =  6.522E-06 

345  K27  =  3.521E-05:REM  K27*C*S  =  1.  (C=l  S=28400) 

350  K28  =  . 00244  ' 

355  K29  =  1.591E-06 

360  K30  =  .01:REM  If  K30*C  =  .01  where  C  =  1 

375  ft\  QQ3-521E"07:REM  3-5E"7  sin<*  K27*C*S  =  .01  (C=1,S=28400) 

400    REM  Scaling  Factors 

410  TO  =  .01:  REM  When  TO  =  .1,  3650  days  =  10  years 

420  DT  =  10  y 

430  M0  =  .007 

440  W0 -.007 

450  NO  =  1.75 

460  SO  =  .002 

470  B0  =  .002 

480  CO  =  6 

490  10  =  .0225 

495  Pr0  =  40 

500    REM  Equations 

504  ll  =  K%IZZ  +KX^yiEEM  T°tal  Gross  Eduction  Rate 

505  RFM  fJi  ?  W  +  f5  R;N  M  +  K8*W  +K16*M  +K22*S+K27*S*C 
i>U5  REM  Total  Respiratory  Consumption 

506  PR  =  Pg/Rp  r 
508  I  =  Im  +  Ir*SIN  (.017*t) 

510    R  =  1/  (1+K1*N*W+K2*N*M) 

520    DW=  K4*N*R*W-K3*N*R*W-K7*W*H-K10*W*H-K12*WK8*W 

L3n    Sif =  ^N^KS'N'R'M-KM'M-Kl^M  W 

565o  Sb:^^^Tk15*m-^s*c-^ 

570  DC  =  K31*S*C-K30*C 

580  W  =  W+DW*DT 

590  M  =  M+DM*DT 

600  N  =  N+DN*DT 

610  S  =  S+DS*DT 


Table  C-2  -  continued. 


Line    Code 


620  B  =  B+DB*DT 

630  C  =  C+DC*DT 

700  REM  Plotting 

710  PSET  (t*T0,  90-1*10) 

720  PSET  (t*T0,190-W*W0) 

730  PSET  (t*T0,190-M*M0) 

740  PSET  (t*T0,90-N*N0) 

750  PSET  (t*T0,290-S*S0) 

760  PSET  (t*T0,290-B*B0) 

770  PSET  (t*T0,410-C*C0) 

780  PSET  (t*T0,  330-(PR-l)*PrO) 

800  t  =  t  +  DT 

805  IF  t*T0<400  GOTO  500 
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Figure  C-l.    One-hundred-year  simulation  of  Biosphere  2 
rainforest  production  and  biomass  minimodel,  B2METAB, 
with  light  increased  20%  over  baseline.  Final  values  for 
plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass. 
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Figure  C-2.  One-hundred-year  simulation  of  Biosphere  2 
rainforest  minimodel,  B2METAB,  with  no  weedy  biomass 
Final  values  for  plant  biomassare  on  the  right  side  of  their 
graphs.  g/m2  =  grams  per  square  meter,  dry  weight  for 
biomass. 
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Figure  C-3.  One-hundred-year  simulation  of  Biosphere  2 
rainforest  production  and  biomass  minimodel,  B2METAB, 
where  all  pruned  biomass  is  put  onto  the  soil.  Final  values  for 
plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass. 
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Figure  C-4.  One-hundred-year  simulation  of  Biosphere  2 
production  and  biomass  minimodel,  B2METAB  with  with 
no  pruning  by  humans.  Final  values  for  plant  biomass  are 
on  the  nght  side  of  their  graphs.  g/m2  =  grams  per  square 
meter,  dry  weight  for  biomass. 
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Figure  C-5.  One-hundred-year  simulation  of  Biosphere  2 
rainforest  production  and  biomass  minimodel,  B2METAB 
where  human  effort  is  reduced  to  .75  of  baseline.  Final  values 
for  plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass. 
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Figure  C-6.  One-hundred-year  simulation  of  Biosphere  2 
rainforest  minimodel,  B2METAB,  where  human  effort  in 
pruning  is  reduced  to  .85  of  baseline.  Final  values  for 
plant  biomass  are  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass 
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Figure  C-7.    One-hundred-year  simulation  of  Biosphere  2 
rainforest  production  and  biomass  minimodel  with  airflow 
cut  off  from  the  rest  of  the  Biosphere.  Final  values  for  plant 
biomass  are  on  the  right  side  of  their  graphs. 
g/m2  =  grams  per  square  meter,  dry  weight  for  biomass. 


APPENDIX  D 
PRODUCTION  AND  DIVERSITY  MINIMODEL 

Table  D-l.      EXCEL  spreadsheet  used  to  calculate  coefficients  for  Biosphere  2 
rainforest  production  and  diversity  minimodel  SPDIV. 

Table  D-2.      Program  in  QUICKBASIC  for  the  simulation  of  production  and 
diversity  in  the  Biosphere  2  rainforest,  model  SPDIV  in  Figure  3-28. 

Figure  D-l.  Simulation  of  Biosphere  2  rainforest  with  continuous  additions  of  5 
and  20  species  per  year  over  10  years.  Final  values  for  plant  biomass  and  species 
meter"  ^  °f  ^^  ^^  ^/m2=^ams  of  drY  biomass  per  square 

Figure  D-2.    Simulation  of  Biosphere  2  rainforest  with  continuous  additions  of  5 
and  20  species  per  year  over  100  years.  Final  values  for  plant  biomass  and 
species  are  on  the  right  side  of  their  graphs,  g/  m2=grams  of  dry  biomass  per 
square  meter.  r 

Figure  D-3.     Simulation  of  Biosphere  2  rainforest  when  50  species  are  added 
after  the  second  year.  Final  velues  for  plant  biomass  and  species  are  on  the  right 
side  of  their  graphs,  g/  m2=grams  of  dry  biomass  per  square  meter 
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Table  D-2.  Program  in  BASIC  for  the  simulation  of  production  and  diversity 
in  the  Biosphere  2  rainforest,  Model  SPDIV  in  Figure  3-25. 


Line  Code  " 

10  REM  SPU1V  -  2  May  1999 

1 1  REM  Species  diversity  decline  in  Biosphere  2 

12  REM  Set  for  25  years' run 

13  REM  started  with  high  transient  species,  low  established  plants 
15  LINE  (0,5)  -  (360,80)„B 

15  LINE  (0,100)  -  (360,180)„B 

15  LINE  (0,200)  -  (360,300)„B 

44  REM  Outside  Sources 

45  Im  =  2600!:REM  mean  insolation  about  which  sine  wave  varies 

46  Ir  =  500.-REM  annual  range  of  insolation 

47  S  =  l 

48  Ja  =  0 

50  REM  Initial  Storages 

55  B  =  1000  :  REM  Starting  biomass 

60  N  -  300:  REM  Starting  transient  species  diversity 

70  E  =  10:  REM  Established  species  at  start 

on  x  -  30'  REM  Canying  caPacity  for  species  in  system 

100  REM  Coefficients 

110  Kl  =  .0346 

120  K2  =  .  000692 

130  K3  =  .  000256 

140  K4  =  7.69E-07 

150  K5  =  3.08E-06 

160  K6  =  2.78E-07 

170  K7  =  2.78E-07 

180  K8  =  2.67E-07:REM  calibrated  for  N  =  300 

190  K9  =  2.67E-07:REM  same  as  K8 

200  K10  =  . 00105 

210  Kll  =  .000855 

220  K12  =  4.04E-08 

230  K13  =  2.85E-11 

240  K14  =  .  0000171 

400  REM  Scaling  Factors 

4^o  St -m  REM  When  T°  =  X  365°  days  =  10  years 

430  B0  =  .01 

440  NO  =  .2 

450  E0  =  .2 

460  10  =  .02 


328 


Table  D-2  —  continued. 
Line    Code 


5U0  REM  Equations  ' " 

508  I  =  Im  +  Ir*SIN  (.017*T) 

510  R  =  I/(1+K1*E*+K2*B) 

520  DB=  K3*E*R+K5*B*R-K4*B*R-K11*B-K12*E*E*B 

530  DN  =  Ja-K10*N-K8*N*B 

540  DE=  K9*N*B-K14*X*E-K13*X*B*E*E:REMK6  =  0 

551  IFE>CTHENX=1 

552  IFE<CTHENX=0 
580  B  -  B  +DB*DT 
590  N  =  N+DN*DT 

595  IFN<0THENN  =  0 

600  E  =  E+DE*DT 

700  REM  Plotting 

710  PSET  (T*T0,  80-1*10) 

720  PSET  (T*T0,180-B*B0) 

730  PSET  (T*T0,300-N*N0) 

740  PSET  (T*T0,300-E*E0) 

800  T  -  T  +  DT 

805  IFT*T0<360  GOTO  500 
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Figure  D-l    Simulation  of  Biosphere  2  rainforest  with  continuous 
additions  of  5  and  20  species  per  year  over  10  years.  Final 
values  for  plant  biomass  and  species  are  on  the  right  sides  of 
their  graphs.  g/m2  =  grams  of  dry  biomass  per  square  meter 
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Figure  D-2.  Simulation  of  Biosphere  2  rainforest  with  continuous 
additions  of  5  and  20  species  per  year  over  100  years.  Final 
values  for  plant  biomass  and  species  are  on  the  right  sides  of 
their  graphs,  g/  m2  =  grams  of  dry  biomass  per  square  meter 
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Figure  D-3.  Simulation  of  Biosphere  2  rainforest  when 
50  species  are  added  after  the  second  year.  Final 
values  for  plant  biomass  and  species  are  on  the  right  side 
of  their  graphs.    g/m2  =  grams  of  dry  biomass  per  square 
meter. 
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